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XI. The Origin of Fatigue Fracture in Copper 


By N. THompson, N. Wapswortu* and N. Lovatt 
H. H. Wills Physical Laboratory, University of Bristol 


[Received December 12, 1955] 


SUMMARY 


The paper describes metallographic observations of the development of 
fatigue cracks in copper. The specimens (polycrystal and single crystals) 
were tested in push-pull at 1000 ¢.p.s. It was shown that the fatigue crack 
started in a slip band inside a single grain. An electropolishing technique 
showed changes in the character of the slip band after 5° of the life had 
expired. Cracks more than one grain long were present when half the 
life had expired. It was confirmed that the fatigue life could be increased 
by preventing access of oxygen to the surface. The bearing of these 
observations on the problem of the origin of a fatigue crack is discussed. 


§ 1. INTRODUCTION 


THe fundamental problem of the nature of fatigue failure in metals may 
conveniently be divided into two parts—although the division between the 
two may prove to be somewhat nebulous. Firstly, we must ask by what 
process an alternating stress can give rise to a singularity in the medium 
which can be called a crack ; and, secondly, we must investigate the laws 
governing the growth of such a crack. The final stages of crack propaga- 
tion are undoubtedly very rapid and are also accelerated and complicated. 
by the accompanying macroscopic re-distribution of stress. This phase is 
of little fundamental interest. There is a growing body of evidence, how- 
ever, that the rate of growth of a small crack may in fact be very small. 
A recent paper by Head (1953) has dealt with this question. To obtain 
further evidence on the matter, and to make the most useful kind of 
observation on the origin of the crack, it is clearly necessary to detect the 
incipient fatigue crack as early as possible. Ideally, one would like to 
observe the point at which the crack will form, before it has in fact 
appeared. It is possible to restrict the area of search by deliberately 
introducing local stress-raisers (e.g. notches or grooves) but in so doing the 
conditions are necessarily made non-representative. 

In the experiments to be described in this paper, it was possible to single 
out a very few places on the surface of a specimen from which fatigue 
cracks would grow, when less than half of the fatigue life had expired ; 


and of the few sites thus chosen, one would almost certainly give rise to the 
SL a eee ee eee 
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eventual fatigue failure. In some early, exploratory work, use was made of 
artificial stress-raisers to localize the phenomena, but in almost all the 
results to be described, this was not done. It is not claimed that the 
methods or results have any validity outside the particular material on 
which observations have so far been carried out. Nor have the conclusions 
any status more than that of a working hypothesis. But it has been 
thought worth while to publish the results of such experiments and 
speculations as have so far been made, in case they may stimulate other 
workers to confirm or refute them. 


§ 2. APPARATUS AND TECHNIQUES 

When embarking on experiments on the mechanism of fatigue, it was 
considered essential that the stresses be macroscopically uniform over the 
gauge length, so as to avoid at least one of the factors which complicate 
the interpretation of results. For this reason a direct stress push-pull 
machine was favoured. Such a one was designed, built and tested by one 
of us (N. L.) and has been used for all the subsequent work. It is similar 
to the Amsler ‘ Vibrophore ’ but has sufficient points of difference to merit 
a brief description. The specimen is, typically, of the form shown at S 
in fig. 1. It is supported vertically from the top, and carries at its lower 
end a massive coil holder (M, fig. 1) and a formerless coil, C. The coil lies 
in the annular gap G of a permanent magnet, PM, and provides the driving 
force to excite the oscillations. In order to obtain sufficiently large 
forces from this arrangement, it is run at the resonant frequency deter- 
mined by the inertia of the coil holder and associated parts, and the 
compliance of the specimen itself. This frequency is of the order of 1000 
cycles per second : by altering the mass of the coil holder and the dimen- 
sions of the specimen, it is possible to increase or decrease this by almost a 
factor of two. 

In view of the sharpness of the resonance peak, such a system must be 
self-excited, rather than driven by an independent oscillator. The 
controlling signal is obtained from a piezo-electric crystal (quartz) incor- 
porated in the top grip (SM, fig. 1). Part of the same signal is fed into a 
valve voltmeter, and is used to measure, continuously, the stress in the 
specimen. The amplitude of the oscillations is fixed by a limiting device, 
which automatically adjusts the gain of the system so that the output 
signal from the quartz crystal has a pre-determined value. The specimen 
is thus driven at constant stress throughout the test, after having taken 
one or two seconds to stabilize itself at the pre-set level on first switching 
on. When fatigue failure is imminent, the power needed to maintain the 
oscillations rises sharply: this increase in demand is made to operate a 
cut-out, which switches off the machine. This usually happens when the 
crack has spread along 4-4 of the circumference. The endurance is 
calculated from the known frequency of the oscillations and the duration 
of the test, as measured by a mains-driven electric clock. 

The degree of axiality of loading is believed to be good. 'The driving 
coil and its holder are axially symmetrical and are laterally unconstrained 
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in any way. In addition, the shape of the specimen lends itself to accurate 
aligning, since, during its production, all four conical bearing surfaces, and 
the gauge length itself, can be turned between centres without removin 
the work from the lathe. Finally, the method of driving the system at is 
resonant frequency for longitudinal vibrations—which will be aimeno 
from the resonances corresponding to other modes—tends to improve 
the degree of axiality. Tests using resistance strain gauges have failed 
to show any asymmetry greater than the errors of the method of test 
and the position of the origin of fatigue fracture has shown a sratifying 
randomness in azimuth. . 


Fig. 1 
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Diagram of fatigue machine. 


The fracture always occurs within the gauge length and usually near the 
centre of it. This is taken as indicating that the radius of curvature of the 
‘ fillets ’ is so large that no important stress-concentration effects are to be 
found there, and as justifying the use of the nominal stress at the centre of 
the gauge length as the actual maximum stress. Before this can be 
given in absolute units, however, the quartz-crystal and valve-voltmeter 
assembly must be calibrated. The most satisfactory method yet devised 
of doing this has been to observe the peak-to-peak amplitude of the 
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excursion of the coil holder when the machine is running, using a micro- 
scope with an eyepiece scale. A knowledge of the mass of the moving 
parts, and the frequency of the oscillation, then gives the peak force 
immediately. The absolute accuracy of this method is probably not 
greater than 2-3%,, but the relative accuracy of stress readings should be 
two or three times better than this. Observations with resistance strain- 
gauges have given confirmatory results. 

All the observations to be described in this paper have been made on 
commercially pure copper. As the work has extended over a number 
of years, two or three different batches of materials have been involved, 
some O.F.H.C. quality, and some H.C. No significant difference in 
behaviour has been apparent between the different batches, although it 
has been found that the fatigue life for a given stress depends a good deal 
on the previous thermal history of the material. The normal procedure 
for the production of a specimen was as follows: The starting material 
was usually 3 in. diameter rolled bar, as supplied by the makers. The 
specimens were machined to size, the finishing cuts along the gauge 
length being not more than 0-001 in. They were then polished mechani- 
cally down to 000 emery followed by ‘ Brasso’ and ‘Silvo’. With 
each new grade, the polishing was done first circumferentially and then 
longitudinally. Specimens were then electropolished in syrupy phosphoric 
acid, using a potential of 2 volts. Polishing for 45 seconds removed. 
about 2 microns from the surface. Finally they were annealed im vacuo 
for 1 hour at 600°c and furnace cooled: a second light electro-polish 
produced a surface suitable for microphotography. The resulting 
grain diameter was of the order of 0-1 mm, with the usual numerous 
annealing twins. 


§ 3. EXPERIMENTAL OBSERVATIONS 
(a) General 


In the course of the work some data have been accumulated on the 
endurance of the material under the conditions of test. Some of this is 
shown in fig. 2. Curve A shows the relation between fatigue life and 
stress, obtained when the specimens were annealed at 400°c instead of 
the usual 600°c. This was done in an attempt to compare our results, 
obtained at 1000 cycles per second, with data obtained at the N.P.L. 
using a frequency of 33 cycles per second. It was later discovered that 
the fatigue properties of this material depend critically on the annealing 
temperature, when this is in the neighbourhood of 400°c. Thus if any 
speed effect were present, it would readily be masked by changes due to 
small variations in heat treatment. 

Subsequent specimens were annealed at 600°C at which temperature 
the conditions are not so critical. Curves B and C both refer to material 
which had had this treatment, but in different furnaces. One was 
continuously evacuated, and in the other, the specimens were treated 
in sealed tubes containing argon. In the former, the work was cooled 
in the furnace: in the latter the tubes were removed and cooled in air. 
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Jt is thought that the difference between curves B and C is most probably 
due to these differences in heat treatment, although this cannot be 
considered as established with certainty. 

In later paragraphs, whenever comparisons of fatigue life are made, 
they always refer to specimens which were prepared under identical 
conditions, so that this uncertainty is not present. It will be noticed 
from fig. 2 that when conditions are carefully controlled, the amount of 
scatter on the endurance curve is quite small. A number of factors 
probably contribute to this result. The material itself is reasonably 
‘pure and homogeneous: the surface was carefully prepared in a repro- 
ducible manner, and in the gauge length there is a homogeneous, uniaxial 
tensile stress which will nowhere exceed its nominal value by any significant 
amount. The consequence is that a ratio of 1-5: 1 between the fatigue 
lives of two individual specimens tested at the same stress can probably 
be considered to be significant. 


Fig. 2 


DYNES cM? 


x EX-N.PL MATERIAL 400°C ANNEAL 
© EX1.C.1, MATERIAL 500°C OR 600°C ANNEAL RAPID COOL 
e EX-1.C4 ate “ " SLOWCOOL 


10 107 108 10? 
CYCLES 
Endurance curves for material used. 


(b) Microscopic Observations 
When tested at stresses giving a fatigue life of 10°10’ cycles, the 
specimens showed slip markings on the surface at a very early stage. 
These became more numerous and some of them became more intense 
as the test proceeded. - It was found that the energy dissipated per cycle 
in the fatigue specimen increased gradually during most of the test 
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after an initial rapid decrease, and that this increase was correlated with 
the increase in the number of slip bands. Microscopic examination of 
the eventual fatigue crack showed that for most of its length it was 
surrounded by a zone of dense slip bands. The general direction of the 
crack was always normal to the tensile stress. Near its centre there 
was always at least one straight section, about one grain diameter long, 
free from the dense slip. The regions above and below this central line 
often had a density of slip markings somewhat less than the average for 
the specimen as a whole. The centre portion of the crack was parallel 
to any slip lines appearing in the same grain as itself. Figure 3 (Plate 3) 
shows the centre portion of a typical crack, and illustrates these features. 

These observations suggest that the crack originated in a slip-band 
at a fairly early stage of the test. Its presence would give rise to a 
stress concentration at its tip, so causing abnormally dense slip there : 
the crack would advance through this region. In addition, the presence 
of the crack would reduce the stress on either side of itself and so produce 
the observed patches with less than the normal density of slip markings. 
The suggestion that the fatigue crack starts in a slip band is by no means. 
new ; but the present evidence suggests that it must have been present, 
with the mechanical properties of a crack, at a fairly early stage of the 
test. 

The high density of slip markings made detection of a crack in its 
early stages very difficult : in order to restrict the field of search, a small 
shallow groove was filed in the side of a specimen to act as a stress 
concentrator. The surface of this groove was polished mechanically and 
electrolytically in the same way as the rest of the specimen before testing. 
The centre portion of the groove was photographed at intervals throughout 
a fatigue test, which lasted for 3-2 10° cycles. A line following much 
of the path of the eventual crack was clearly visible after 2-15 x 105 
cycles (67°, of the life) and some part of it showed after 6-5 104 cycles 
(20°%,). The place where the crack eventually formed was marked by 
slip-bands as early as 1-5 104 eycles (5°). 

In an attempt to distinguish between a slip-band and the early stages 
of a crack, it was decided to electro-polish the specimen at intervals 
throughout a fatigue test. Preliminary trials established a standard 
‘dose ’ of electro-polish which was sufficient to remove almost all of the 
slip markings ; about 2 microns was removed from the surface by this 
treatment. Another grooved specimen was prepared and the previous 
experiment repeated, the sequence of operations being fatigue, photo- 
graph, electro-polish, photograph, etc. It was found that although 
almost all the slip-bands were removed by the standard electro-polishing 
treatment, some remained, and even appeared to be brought into greater 
prominence by the process. A second dose of electro-polish immediately 
following the first usually left these persistent bands. still present, 
indicating that they were not just abnormally large steps on the surface, 
but were in some way qualitatively different from the remainder, On 
repeating the sequence of operations, it was found that some of these 
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persistent bands increased in intensity, and after continued application 
of the fatigue stress, a few of them eventually developed into cracks 
which spread into neighbouring grains. 

Following this preliminary trial, all subsequent observations were made 
on specimens without any groove, so as to avoid the complication of the 
uncertain stress concentration factor. The use of the electro-polishing 
technique meant that places likely to be of interest were not obscured 
by a welter of normal slip bands, and were less likely to be overlooked. 

After fatiguing for a known time, chosen to be a small fraction of the 
estimated fatigue life, a specimen was electro-polished. Examination of 
the whole surface revealed a number of typical ‘ persistent slip bands ’ 
and the more prominent of these were photographed before the specimen 
was replaced in the fatigue machine. This cycle of processes was 
repeated at intervals until the specimen finally broke after 3-65 10® 
cycles. 

The characteristic persistent slip bands were observed the first time 
the test was interrupted, which was after about 4° of the fatigue life : 
these became more marked as the test progressed, but most of those 
which had been kept under observation from the early stages did not 
spread beyond the grain in which they started. However, figs. 4-7 
(Plate 3) show one such field of view in which progressive changes were 
observed. Figures 4 and 5 show it, respectively, before and after electro- 
polishing, when 74°, of the fatigue life had elapsed. A few persistent 
slip bands are visible in fig. 5, but unfortunately the point which sub- 
sequently proved most interesting is just off the edge of the negative. 
Figures 6 and 7 show the same field after electro-polishing when 42°% 
and 77%, respectively of the life had elapsed. The marks are more 
numerous in fig. 6, and by fig. 7 one of them has developed into a crack, 
extending into neighbouring grains. 

The crack which eventually led to the fatigue failure in this specimen 
was first photographed after 42°, of the life (fig. 8, Plate 4). At that 
stage it was the only one of the many similar markings observed which 
had already shown signs of spreading into neighbouring grains—a 
behaviour similar to that shown at a later stage by the other crack 
already illustrated. Figures 9 and 10 (Plate 4) show this same crack 
at 77%, and almost 100% of the life.* Taken together, the two series 
of photographs present a fairly complete picture of the stages in the 
development of a fatigue crack under the conditions of these observations. 
Persistent slip bands were visible after 74°, of the life ; one had developed 
into a crack involving neighbouring grains after 42% of the life, and the 
most serious crack was already many grains long after 77°% of the life. 
Similar slow rates of growth were shown for other cracks both on this 
specimen and on other specimens. It has also been verified that the 
general course of events is not affected by the number of times that the 


i rr ell esse Eres 
* It should be noted that figs. 9 and 10 are both mosaics, and that the 
specimen was not electro-polished at the conclusion of the test, before fig. 10 was 


taken. 
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specimen is polished during a test; nor is the fatigue life significantly 
altered by these interruptions. The fatigue process, 1m short, is not 
being modified by the steps taken to observe it. 


(c) Tensile Tests 
It is still far from obvious at what stage a ‘ persistent slip band © may 
- properly be called a ‘crack’. In an attempt to throw more light on this 
problem, a partially fatigued specimen was chosen which showed a 
number of ‘ cracks ’ ranging from 3 to 12 grains long, and many * persistent 
slip bands’ extending across a single grain. The specimen was then 
given a 5°, slow tensile extension, and re-examined. The cracks more 
than one grain long had all opened into large cavities, but none had 
shown any tendency to extend lengthways. Persistent slip bands 
confined to one grain had not opened up by any detectable amount, but 
there is a suggestion of some broadening of those that were just beginning 
to invade the neighbouring grains. Figures 11 and 12 (Plate 4) show 
one such borderline case before and after the 5°, strain. This experiment 
would suggest that the grain boundaries play an important role in the 
early history of a fatigue crack. Such a view also receives some support 
from the results already quoted, that while persistent slip bands were 
visible after 4°, of the life of one specimen, they were only beginning to 
invade neighbour grains at about 40°,,—and, in this one case at least— 
it was the crack which first spread in this way which led to the final 
failure. 
(d) Haperiments on Single Crystals 

In view of these results, some interest attaches to observations made 
on single crystals. Four such specimens have been examined so far. 
The crystals were grown from the melt by the usual techniques, carefully 
machined to size and then treated in the same way as the polycrystals 
in the matter of surface preparation. The phenomena observed were 
qualitatively exactly parallel to those shown by polycrystals, and the 
observations are not sufficiently numerous to justify any close comparison - 
of details. Figures 13 to 18 (Plates 5 and 6) show some photographs 
of one of these specimens at various stages of its life.* Apart from 
figs. 14 and 18, the specimen had been electro-polished just prior to the 
photographs being taken. This crystal was orientated so that slip took 
place preferentially on one system, although the stress on a second 
system was only slightly less. (A number of curious patches of slip on 
this second system were observed ; fig. 14 shows one such.) For the 
most part, however, slip was confined to the single system. Some of 
these slip markings persisted after an electro-polish at 18°, of the eventual 
fatigue life (fig. 13). The persistent marks were rather short, typically 
about 50 microns long and arranged in rows along the slip bands. 
Further stressing showed a tendency to increase the intensity and length 


In these and all other cases the tensile axis of the specimen is parallel to the 
long edge of the individual frames of the mosaic. 
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of some of these markings—or it may be that the shorter ones in or near 
the same slip band joined together (figs. 15 and 1636 anda o. 
respectively). After 90% of the life, it was seen that some of the persistent 
slip marks had developed. into cracks which in places did not exactly 
follow the slip bands (fig. 17). This was much more obvious when the 
test was finally stopped after a further 10°/ when numerous large cracks 
were present (fig. 18). For appreciable parts of their length the cracks 
followed the direction of the first or—more rarely—the second slip 
system, but the general direction of the main failure, and of other, shorter 
cracks in the same specimen was the usual one, i.e. normal to the direction 
of the tensile stress. Another single crystal specimen, which was run to 
failure without any intermediate polishing, showed the same general 
form of crack—initially following the active slip plane but later deviating 
from it, to run approximately normal to the tensile stress. 

An interesting feature on fig. 14 is the traces of slip on the second system 
which have appeared at the ends of the persistent slip bands left by the 
previous polish. (Figure 19 (Plate 7) shows an example, at greater 
magnification, from the fourth single crystal specimen which was also 
run to failure without interruption.) Their presence clearly indicates 
that in these regions the local stress is abnormal. Such a stress concen- 
tration might quite well be produced at the tip of a crack ; but similar 
effects would also be expected from, for example, a closely packed group 
-of dislocations. 


(e) Depth of Persistent Slip Bands 

Two further experiments should be mentioned concerning these 
persistent slip markings and their relation to fatigue fracture. The 
first concerns their surface origin. Although they are revealed in the 
first place by a light electro-polish which removes about 2 microns from 
the surface, continued electro-polishing will eventually remove them 
completely. Figures 20 to 22 (Plate 7) show a series of photographs of 
a partially fatigued specimen which illustrate this point. The same 
field of view is shown after electro-polishing to remove about 2, 10 and 
18 microns. A further photograph, taken after 34 microns had been 
removed, showed no trace of the marks. Experience suggests that, in 
the first half of fatigue test on a polycrystalline specimen, most of the 
persistent marks are fairly shallow, and that the deepest of them does 
not penetrate more than about 100 microns. This is the same order of 
magnitude as the grain diameter on the material used. It is noteworthy, 
that in all the observations of this kind which have been made, the 
removal of metal from the surface has never revealed any new markings 
of the same kind ; all of them have their origin on the free surface. 

Nevertheless, even at an early stage of the test the persistent slip bands 
represent a source of weakness, and contain within them the seeds of 
fatigue failure. A striking illustration of this point was provided by the 
following experiment: A fatigue test was stopped after 25% of the 
estimated endurance, and the specimen was heavily electro-polished 
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until all the persistent markings had been removed. This necessitated 
the removal of a layer of metal about 30 microns thick. The test was then 
continued at the same stress as before, for a further, equal time, and the 
polishing repeated. This was done nine times in all, so that the specimen 
had survived 225°, of its normal life. In view of the small amount of 
scatter obtained ‘in our observations, this is quite a significant increase. 
Moreover, even after this treatment, microscopic observation showed no 
difference between this specimen and a piece of virgin material, and there 
appeared to be no reason why the process could not be continued inde- 
finitely, as long as any metal remained. 


({) Bxperiments on Annealing 

Finally. we might mention some experiments on the effect of periodic 
annealing during the course of a fatigue test. Sinclair and Dolan (1951) 
have carried out an extensive investigation of this kind on polycrystalline 
z-brass. They annealed their specimens for 30 minutes, at 400°C, every 
20°, of the expected fatigue life. They found that the endurance of 
such specimens was statistically indistinguishable from that of untreated 
specimens. 

We have made similar observations on three polycrystalline copper 
specimens. ‘The annealing conditions were | hour at 600°c—a treatment 
which completely re-softens the material without producing any noticeable 
grain growth or recrystallization. One specimen was given this treatment 
after 50°, of its expected, unannealed life, and one was given it every 
20%. Both of them broke after a life not significantly different from that 
of normal specimens. 

The third specimen was annealed every 10°, of its expected life. 
After the first run, it was electro-polished and showed the usual persistent 
slip bands. After annealing. followed by a second electro-polish, the 
same markings still persisted. The sequence of fatigue, polish, anneal, 
polish, was repeated five more times. A photograph taken at this stage 
showed persistent markings already spreading into neighbouring grains, 
the general appearance at this stage being intermediate between that 
shown by another, unannealed, specimen photographed at 42°, and 
77%, of its actual life (figs. 8 and 9). There is little doubt, in fact. that 
the events were following their normal course, unaffected by the annealing. 
Unfortunately the specimen had been slightly dented while being removed 
from the fatigue machine for the first anneal; a crack formed at the 
bottom of this dent, and the specimen broke there, mid-way through 
the seventh fatigue period. . 


§ 4. Discussion 
To summarize, it appears that after only a few per cent of the fatigue 
life has elapsed, some of the numerous slip bands undergo some modifica- 
tion of their structure which causes them to react differently to the 
electro-polishing conditions. The change, whatever its nature. originates 
on the free surface of the specimen. At some stage, the modification 
confers upon the slip band the essential property of a crack—a surface 
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across which it is not possible to transmit a tensile stress. A few of the 
slip bands go through this sequence more rapidly than the majority, and 
have developed into obvious fatigue cracks, several grain diameters long, 
before the longest of them is so large that the specimen can be said to 
be broken. The modification to the slip band structure cannot be 
reversed by annealing.* 

It has been noticed that, of these persistent slip bands which develop 
into cracks, an abnormally high proportion lie very close to twin bounda- 
ries. The explanation of this is not clear: but the proximity of the 
twin boundary cannot be an essential feature in the formation of the 
crack since the events follow the same general pattern in single crystals. 


(a) Effect of Oxygen Atmosphere 


The observations are consistent with the hypothesis that the slip 
band changes directly to a series of micro-cracks, but it is not easy to 
conceive of a detailed mechanism for this which requires, as an essential 
feature, that the change must originate on a free surface. As an alter- 
native, therefore, we have explored the possibility that the diffusion of 
oxygen from the atmosphere into the slip band forms an intermediate 
step in the process. The persistent slip bands in their earliest stages 
might then well be regions of abnormally high oxygen concentration. 
This explanation, too, is consistent with all the observations so far 
recorded. If it is correct, one would expect that the endurance of a 
specimen might be increased if access of oxygen to the surface were 
prevented. It was found that just such an effect had been reported a 
number of years ago by Gough and Sopwith (1932, 1935, 1946). They 
showed that the life of a copper specimen was prolonged if, during the 
fatigue test, it was placed in a partial vacuum (10-* mm Hg), or in an 
atmosphere of dry, oxygen-free nitrogen. 

It was thought desirable to repeat the experiments of Gough and 
Sopwith to verify that the effect was obtained with our material, and 
under our conditions of test. The construction of the apparatus did 
not readily lend itself to experiments im vacuo, but an atmosphere of 
nitrogen could readily be provided by surrounding the specimen with a 
light polythene bag, more or less air-tight, and maintaining in it a 
pressure of gas slightly greater than atmospheric. The nitrogen was 
obtained from a cylinder, and was claimed by the suppliers to contain 
not more than | part in 10° of oxygen or water vapour. Six specimens 
from the same batch were tested, all at the same stress, either in nitrogen 


* Some reference should be made at this point to an early paper by Ewing 
and Humfrey (1903) which is in danger of being forgotten. In this they describe 
microscopic observations on flat specimens of Swedish iron, fatigued in rotating 
bending. The sequence of events described is closely similart to that found by 
us for copper: slip lines growing broader, more diffuse and more numerous. 
with cracks eventually appearing in such broadened slip bands. ‘The differences 
in material, in stress system, and in technique help to establish the facts more 
clearly in the case of copper but it would appear that, at least in broad outline, 
the story is the same for iron. 


124 N. Thompson, N. Wadsworth and N. Louat on the 


or in air. The values obtained for the fatigue life are shown in the 


table : 
CR LOA VOOM PRS eS ee ee 


Specimen No. Atmosphere Life 
58 Air Bite a 
66 Nitrogen 20-3 
67 Air 3-4 
68 Nitrogen 10-4* 
69 Air 3:2 
70 Nitrogen 17-5 


* The fatigue crack started at one of a number of very small depressions on 
the surface. These had presumably been overlooked when the specimen was 
originally polished. This was the only occasion on which they wereseen, and their 
presence could well account for the short life of No. 68, compared to Nos. 66 and 70. 


There is little doubt that the life in nitrogen is considerably greater 
than that in air, particularly in view of the good reproducibility normally 
obtained. The average increase in life—of the order of 5-6 times—is 
also in rough quantitative agreement with Gough and Sopwith, who found 
an increase in fatigue ‘limit’ of 13° (calculated on the basis of an 
endurance of 3x10? cycles): the two figures are consistent with the 
slope of the stress-log (cycles) curve. 

With specimen No. 70 the test was interrupted at intervals for the 
usual electro-polishing and surface examination. Figures 24 to 27 
(Plate 8) are a selection of the photographs taken. At about 5% of the 
expected life in air, the surface appearance was not very different from 
that which would have been found if the test had been running in air 
(fig. 24). By 42°, of the expected life, there was little change (fig. 25). 
At 164% the persistent marks had developed further (fig. 26); this 
point proved to represent 30° of the actual life in nitrogen. By 400%, 
there was little further change in the length of most of the persistent 
slip bands that had been kept under observation although they were more 
numerous (fig. 27). However, at one place a crack in a more advanced 
stage was discovered, and this eventually proved fatal. It is not easy 
to speak with confidence about qualitative observations of this kind, 
but the general impression gained from a study of the numerous photo- 
graphs taken was that the absence of oxygen had little effect on the 
first appearance of the persistent slip bands, and that the major effect 
was to slow down the transition from this stage to the point at which 
there could be no doubt that the bands had turned into cracks. Further 
propagation of the crack was slowed down also. 


(b) Role of the Oxygen 


In the light of these results, there is little room for doubt that the 
pressure of oxygen in the atmosphere around a copper specimen has an 
effect on the fatigue life. Further, the effect of the atmospheric oxygen 
seems to be intimately concerned with the transition from slip band to 
crack, and with the early stages of crack development. 
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Tt is well known that an adsorbed mono-layer of oxygen forms very 
rapidly on a fresh copper surface and that, once formed, it is very difficult 
to remove (see, for example, Schlier and Farnsworth 1954). A rough 
calculation suggests that there would be ample opportunity for such a 
layer to form on the freshly exposed surface in the slip band during one 
half cycle of a fatigue test. On the following half cycle, one can envisage 
three possibilities : (i) The oxygen will be ‘scraped off’ and slip will 
take place in the reverse sense on the same plane; this would appear to 
be somewhat improbable. (ii) The slip in the reverse sense will take 
place on some neighbouring plane : this could well be the beginning of a 
process leading to extrusion of metal from the slip bands, as reported 
already by Forsyth (1955). Such an effect has, in fact, been observed 
by us on a single crystal of copper (figs. 23 (a) and (b), Plate 7). The loss 

of this material from the interior of the crystal just below the surface 
would give rise to an effective cavity, the stress field of which would 
represent a potential source of a crack. (iii) Slip in the reverse sense 
will take place to some extent on the same plane, and some of the freshly 
adsorbed oxygen will be drawn in to the body of the crystal in the form 
of dissolved atoms. This oxygen, replenished by fresh supplies drawn 
into the crystal in the same manner each cycle, could then diffuse down 
the slip band, building up an abnormally high concentration of oxygen 
in and near the band and weakening the crystal so much that a crack 
eventually forms there. 

On this latter view the persistent slip bands would be regions of high 
oxygen concentration. It is at once clear that no ordinary diffusion 
process could account for the penetration of the persistent bands to an 
appreciable fraction of a grain diameter (~10-2cm) in times of the 
order of 10? sec. This would require a diffusion coefficient of order 
10-? cm?/sec—which is only reached for copper at about 800°c (Ransley 
1939). However, examples are multiplying of metallurgical changes 
occurring under fatigue conditions which are explicable if one postulates 
abnormally fast diffusion (see, for example, Holden 1954, Hanstock 1954). 
Several possible mechanisms suggest themselves, such as an increased 
vacancy concentration, and high local temperatures due to ‘thermal spikes’ 
(Seitz 1952) ; such efforts will be most pronounced near an active slip band. 

An objection to any mechanism in which diffusion is the rate-deter- 
mining process is that it would tend to make fatigue a very temperature 

_ sensitive phenomenon, and would predict a considerable frequency effect, 
since the production of a given state of the crystal would depend on the 
duration of the test in seconds rather than in cycles. Neither of these 
effects is well marked in fatigue phenomena in general. If, however, 
diffusion were merely capable of equalizing the oxygen concentration 
across a slip plane at the end of each half-cycle the relative movement 
of the two parts of the crystal during slip would result in oxygen being 
drawn down the slip band. In this case the number and extent of the 
cycles of slip would be the rate-determining process, and the above 
objection would be removed. 
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However, any hypothesis in which the oxygen plays an essential role 
will have difficulty in accounting for the fact that fatigue fracture takes 
place in copper at 4°k (Rosenberg, private communication). An alter- 
native interpretation of our results which avoids this trouble would be 
that the persistent slip bands already have the character of minute 
cracks at a very early stage (say 5% of the total life, orless). This 
~ evades the question of how the crack forms in the first place and why it 
forms only on the surface, and assigns to the oxygen the part of expediting 
the growth of the crack in its earliest stages. Certainly the formation 
of an adsorbed layer of oxygen on the faces of a freshly opened crack 
would tend to prevent it closing again on the next half-cycle. If, following 
this line of thought, a persistent slip band in its earliest stages consists 
of many small cracks, each too small to be seen individually, then its 
appearance would not be expected to change appreciably due to any 
opening of the cracks under the influence of a steady tensile stress. 

On either of the two views one can see that a moderate annealing 
treatment would not suffice to undo the damage already done to the 
crystal early in the fatigue test. Thus both hypotheses are in accord 
with the observations which we have made to date, and more experimental 
work is needed before the mechanism can be clearly established. 
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ABSTRACT 


A survey is presented of our present knowledge concerning radiative 
transitions in nuclei of A<20. Over a hundred transitions, chiefly dipole, 
are analysed in terms of the Weisskopf units, and the following chief 
conclusions are drawn : 

(i) E1 transitions have a most probable speed of about 0-032 Weisskopf 
units with a spread in speed of about a factor of seven either way ; 

(ii) The corresponding quantities for M1 transitions are a speed of 
0-15 Weisskopf units and a spread of a factor of 20 either way ; 

(iii) These distributions are consistent with shell model expectations ; 

(iv) If a transition shows I)(ev)/H,3 (Mev)>0-02, then it is very 
probably E1 ; 

(v) There is no clear dependence of transition speed on level spacing 
or other feature ; 

(vi) Some of the higher multipole transitions are inconsistent with shell 
model expectation and give evidence of collective motion. 


§ 1. INTRODUCTION 


Two years ago a survey of radiative transitions in light nuclei (Wilkinson 
1953 a—referred to as I) revealed some 22 dipole transitions, equally 
divided numerically between electric and magnetic in character, which had 
been established with fair to good probability, and whose widths were 
known with accuracies varying from a few per cent to a factor of two or 
three. 

From this scanty material it emerged that El transitions were by no 
means weak as a rule, and indeed possessed strengths of about one- 
hundredth to one-tenth of the rough single-particle shell model estimate 
that has come to be called the Weisskopf unit (Weisskopf 1951—for a 
detailed treatment see Moszkowski 1955) : 

P,»(H1)=0-11A?" #,? ev. 
HE, is the energy of the transition measured in Mev, and we have used a 
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nuclear radius of 1:5 Ax 10-13 em.* The radiative width of a transi- 
tion in these Weisskopf units is then a measure of the associated matrix 
element, and is written as | M |?. It also appeared from I that the spread 
in the values of | M |? was not very great for E1 transitions that did not 
suffer the inhibition of the isotopic spin selection rule, but so few transitions 
were involved that little confidence could be felt on this point. It was 
pointed out by Lane and Radicati (1954) and emphasized by Lane and 
Wilkinson (1955) that, even within the pure shell model, we must expect 
transition strengths relative to the single-particle strength to vary widely, 
ranging from very small values (including zero) to values of greater than 
unity. This wide departure to be expected from the single-particle width 
or the Weisskopf unit is due to five causes : 

(i) In areal transition we are not concerned as a rule with the behaviour 
of one isolated particle, but rather with the behaviour of the often many 
(n) particles needed for a full specification of the states between which the 
transition takes place. Since we do not know which particle actually 
makes the transition, we must multiply the one-particle matrix element 
by n. This leads to a speeding up of the transition, and accounts for the 
fact that some theoretical shell-model transitions have | M |? values of 
greater than unity. We may think of this enhancement effect as due to 
the correlations imposed on the motions of equivalent nucleons by the 
Pauli principle. 

(ii) The requirement that radiative transitions proceed only through 
common parents gives a diminution of the matrix element except in rare 
cases, and can lead to the transitions being completely forbidden. 

(ii) In general the (Racah) coefficients that describe the various vector 
couplings taking place in the transition lead to a diminution, often 
substantial, of the matrix element. 

(iv) The ‘ potential well ’ within which the transition takes place is not 
infinitely massive and uncharged as is assumed to calculate the single- 
particle estimate or Weisskopf unit, but rather has a charge and is of finite 
mass. This means that the ‘ well ’ recoils from the motion of the radiating 


* There is some evidence from the treatment of the Coulomb energy of 
light nuclei, taking into account exchange contributions (see e.g. Cooper and 
Henley 1953, Carlson and Talmi 1954), that a better value for the constant 
in this expression for the radius would be about 1-3. This is in line with 
other familiar evidence in heavy nuclei that the radius constant is about 1-2. 
However, the evidence from the dynamical properties of light nuclei seems to 
favour the larger radius for computations of the present sort. Use of the 
larger radius in computing the Weisskopf units also makes some allowance 
for the fact that, in most of the transitions with which we concern ourselves 
the radiating nucleon is unbound in its initial state, and so gives a contribution 
to the matrix element from outside the nucleus which does not occur in the 
computation of the Weisskopf unit. The extra contribution for the cases of 
interest here usually speeds the transition so the smaller radius with the 
proper wave functions gives about the same result as the Weisskopf unit 
computed for the big radius. 
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particle and, since the well is charged, this constitutes a current opposing 
that of the radiating particle if it is a proton. In a classical calculation 
this results, for E1 radiation, in replacing the charge e of the radiating 
proton by the familiar ‘ effective charge ’ (1—Z/A)e, or, if the ‘ radiating 
particle * is a neutron, by (Z/A)e. In a proper shell model treatment the 
situation is more complicated, but a similar effect obtains, and so we see 
that, even neglecting the many-particle complications (i), (ii), and (111) 
above, E1 transitions of a single particle should not be expected to display 
strengths of one Weisskopf unit, but of only about (N/A)? or (Z/A)2—one- 
quarter of a Weisskopf unit. This point is much less important for 
magnetic transitions, while for electric transitions of higher order, for a 
radiating proton, the expected strength is about (1—Z/A*)? Weisskopf 
units (for a 2*-pole transition), which is very close to unity, and for a 
‘radiating neutron ’ the expected strength is about (Z/A")?, which is very 
close to zero. 


(v) The Weisskopf unit is correct only for transitions ending in a state of 
zero orbital angular momentum and where the initial state has J=k-+-}. 
If this is not so the expected width even in the idealized one-particle case is 

different from the Weisskopf estimate. The difference may be easily a 
factor of three or more in a practical case. This is discussed in more detail 
later. Another difference is that between the radial overlap integral 
evaluated with ‘realistic’ shell model wave-functions and with the 
rectangular wave-funtcion of the Weisskopf unit. This latter difference is 
not great in the cases of present interest. 

To avoid confusion we should state that throughout this survey | M |? in 
connection with an experimental transition means always the actual 
radiative width I’, divided by the Weisskopf estimate I’,,, appropriate 
to that type (multipole and parity change) of transition as written out 
above for El transitions and later for others. We do not, in quoting 
| VM |?, apply any ‘ correction factors ’ for any of these five points just 
enumerated. 

Such considerations show that the observed dispersion of | M |? values is 
qualitatively quite consistent with transitions taking place within the pure 
shell model ; this point we return to later. 

The M1 transitions analysed in I were also of very roughly the strength 
expected from shell model type calculations. Here the appropriate 
Weisskopf unit has the value : 

Ty (M1)=0-021 EB, ev. 

It appeared both that the spread in | M |? values was somewhat greater 
for M1 than for E1 transitions and that these | M |? values were somewhat 
higher for the magnetic than for the electric transitions. We might note 
that if considerations (i), (ii), and (iii) above are similar for M1 and for El 
transitions then, ceteris paribus, we should expect point (iv) to result in 
values of | M [? for El transitions running, in the mean, about a factor 
four less than those for M1 transitions. 
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The spreads in | M |? values recorded in I were so great that it was felt 
that no safe distinction could in general be made between El and M1 
transitions on the basis of radiative width alone. 

The intervening years since the publication of I have seen a great increase 
in the volume of material pertaining to the radiative properties of light 
nuclei, and it is now possible to analyse some 105 transitions instead of the 
22 of 1953. This is a sufficient bulk of material to begin to treat statistic- 
ally, and from this present survey it emerges that the main tentative 
conclusions of I concerning dipole transitions are substantiated ; they may 
now be stated more quantitatively. A new feature is the beginnings of 
information on transitions of higher multipole order, but as yet data are 
extremely sparse. For those transitions of interest in the present survey 
we have the Weisskopf units (for a radiating proton) : 


,,,,(H2)=1-2 x 10-7 A4? B® ev, 
r,,,(E3)=8-7 x 10-14 A? B7 ev. 


yw 

Also since the publication of I there has appeared a number of detailed 
attempts to apply the shell model to calculate in detail various dynamical 
properties of light nuclei including radiative widths of certain special 
transitions (Lane 1953, 1955). These attempts, on nuclei in the 1p-shell, 
have met with an encouraging success and suggest that the shell model may 
indeed be the basis for a first-order description of certain states of light 
nuclei. Again the attempt by Inglis (1953) to account for the level 
positions of lp-shell nuclei in terms of the shell model shows a sufficient 
relationship to experiment to give the impression that any prediction of the 
shell model in this region of the periodic table must be taken seriously. 
The work of Redlich (19&5) and of Elliott and Flowers (1955) on the nuclei 
just beyond 1°O shows kow agreement in respect both of level positions and 
dynamical properties may also be achieved between experiment and shell 
model theory at the beginning of the 1d—2s shell. These results provide 
encouragement for a continued examination of the predictions of the shell 
model about radiative transitions in light nuclei. This examination may 
take the form of a search for detailed shell model descriptions of specific 
transitions as, in particular, Lane has done, or we may make an omnibus 
comparison of the predictions of the shell model as to strengths and dis- 
persion of strengths of radiative transitions with the whole ensemble of 
experimental fact. This latter approach is pursued in this survey. 

Other more recent publications have confirmed the feeling that much 
qualitative if not quantitative truth resides in the shell model’s view of 
radiative transitions as summarized in the Weisskopf units. Thus the 
initial survey of isomeric transitions by Goldhaber and Sunyar (1951), : 
which was in large measure responsible for triggering off the present 
considerable interest in the dynamical predictions of the shell model, has 
been extended (Goldhaber and Sunyar 1955) and confirms the earlier 
impression that the gross structure of isomeric transitions conforms to the 
general pattern of shell model expectation. Detailed comparisons are, 
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however, difficult to make in this field of the heavier nuclei, and in 
particular it can be seen from the work of Goldhaber and Sunyar that in 
certain respects the simple shell model breaks down—for example, there is 
little sign that selection rules specific to the shell model are well obeyed, 
and again there is little difference evident between the strengths of odd- 
proton and odd-neutron electric transitions of high order. To some 
degree such behaviour may be understood in terms of configuration 
interaction, but it seems likely that some more far-reaching revision will 
be needed. At one point the shell model seems to fail rather completely, 
and that is for E2 transitions in heavy nuclei which frequently have values 
of | .M |? so much greater than unity (up to 100 or more) that it appears 
rather artificial to attempt a shell model interpretation, and the collective 
model (Bohr and Mottelson 1953) must be accepted as the more realistic 
approximation. Much interest will centre on the systematics of E2 and 
3 transitions in light nuclei when more data become available. The 
present survey can only make a small beginning in this direction but, as will 
be seen, there is some suggestion that, for transitions of this type, we 
may be already confronted by evidences of collective motion even in the 
Ip shell. Collective motion of the type envisaged by Bohr and Mottelson 
does not produce enhancement of E1 transitions and, since considerable 
equilibrium distortions are not to be expected in the lightest nuclei. 
we should not expect shell model predictions to be very greatly modified 
by collective properties of the ‘ rotational’ type; ‘ vibrational ’ oscilla- 
tions are to be expected. On the other hand we must remember the 
possibility of collective oscillations of the ‘ polarization ’ type (Goldhaber 
and Teller 1948), which may result in enhanced E1 transitions. As yet 
there is little direct evidence for such polarization oscillations: the 
‘giant resonances’ of nuclear photodisintegration for which they were 
postulated have been shown to be susceptible of an alternative inter- 
pretation within the frame-work of the shell model (Wilkinson 1955). 
A further recent important contribution to our knowledge of E1 
transitions comes from work on the gamma-ray spectra following capture of 
slow neutrons (Kinsey and Bartholomew 1954, Kinsey 1955), where it is 
found that, where they are possible, E1 transitions predominate and have 
a radiative width that, on the basis of about a dozen transitions, may be 


written : r~0-06Dr 


Y yw 


where D is the spacing (in Mev) of levels, at the excitation from which the 
gamma-ray is emitted, of spin and parity such that they may combine 
by radiation of the given type with the final state. Such a factor (Blatt 
and Weisskopf 1952) is needed to take into account the configuration 
interaction, many-particle excitation and similar effects that play so 
important a rdle in the close crowding of levels at high excitation. ~The 
very rough expectation of the shell model is that, on account of these 
effects, the numerical factor in the above equation should be as an upper 
jimit somewhere between 0-05 and 0-5 (taking into account the effective 
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charge factor in (iv) above). We see therefore that the results of Survey I 
for the light elements and the work of Kinsey chiefly for heavier elements 
agree that the expectations of the shell model so far as El transitions go 
are fairly well realized. Indeed if we insert for D in Kinsey’s expression a 
value of about unity appropriate for the light nuclei of Survey I, we find 
that the properties of transitions in heavy nuclei extrapolated in this way 
into the region of the light nuclei agree very well with what we find in 
practice for light nuclei—transitions whose strength is of the order of a 
tenth of the Weisskopf unit. It is not yet clear how important configura- 
tion interaction and allied effects will be for the light nuleci with which we 
are here concerned, and what place the spacing of levels will have in deter- 
mining radiative properties ; this point will be returned to later. 

Other interesting recent work concerns the total radiative widths of 
individual slow neutron resonances (Hughes and Harvey 1954). It seems 
likely that such results will be more informative about statistical aspects 
of the nucleus such as the dependence of level density upon various 
parameters than they will about the nature of the radiative process itself. | 

We now present the experimental material which, as before, has been 
arbitrarily restricted to A <20. 


§ 2. EXPERIMENTAL MATERIAL 


The experimental material that is the source of this survey is necessarily 
of very variable quality and reliability. In some cases the type of the 
transition (multipolarity and whether electric or magnetic) is established 
with good probability ; such transitions are granted Status A in the 
following tables. In other cases the evidence for one assignment rather 
than another is direct but only fairly good, or good but only circumstantial ; 
here Status B is given. Where the evidence is not strong enough to give a 
clear preference for one assignment over others the transition is ignored 
even though it may be well-established and of accurately-measured 
strength. Where the strength of a transition is known to a few tens of 
per cent or better a star is attached under Status ; strengths of transitions 
not so marked may be in error by as much as a factor of two or three or 
even more. 

In determining the type of a transition ‘ bootstrap ’ arguments have been 
avoided as much as possible. We have very occasionally, however, in 
special cases allowed ourselves to be guided by the strength of a transition 
in determining whether it is M1 or El. As will be seen later, it appears 
that if a transition shows a value of the quantity I(ev) divided by 
E,? (Mev) of greater than 0-02, then there is about a 10: 1 chance that it is 
E1 rather than M1. This criterion has been admitted to assign El] 
character to the transitions in 14N between the following pairs of states— 
8:62-6:23 Mev; 8-70-5-69 Mev; 9-18-6-44 Mev; 9-49-5-10 mMev—which 
show values of the above quantity of 0-10, 0-18, 0:06, and 0-10 respectively. 
These assignments receive Status B (in any case the existence of the 
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second of these is not established with certainty). A case such as the 
transition between the 8-90 and 5-83 mev states of 14N, which has a value 
of the above quantity of 0-014, cannot be classed by the use of this criterion 
and so is ignored. 

It sometimes arises that while the product of the radiative width and the 
statistical weight factor of a disintegration cross section is well known. the 
spin J of the radiating state has only been determined within limits. 
This means that only (2J+-1) | M [? can be quoted with accuracy rather 
than the more interesting | M |?. However, spins are known in a sufficient 
majority of cases now for us to tabulate | M |? rather than (2.J-+1) eats 
as we had to do in I.* In cases of unknown J we have quoted a value of 
I’, and so of | M |? which is the mean of the values derived from the 
various alternative spin possibilities. This has been done in three cases 
in 1°O and one in “Be where the radiating state is in each case J=1/2 or 
3/2. The error introduced by our procedure is 50° or less, which is not 
worrying in view of the large dispersion in the distribution of | M |? 

values to which these transitions are contributing. In order to enable 

corrections to be made in the event of re-assignments we list in the 
following tables the value or values of the spin J of the initial (emitting) 
state that we have assumed in computing I’, from the experimental 
results. 

In one case, the probable M1 transition from the 18-14 Mev state of ®Be, 
it is not known to which of two states the transition goes, the ground 


* We shall for the greater part of this survey eschew any explicit consideration 
of statistical weight factors as such. The subject is rather a vexed one and it is 
the responsibility of each author to justify his own use of any particular 
statistical weight factor or set of symbols likely to be mistaken for a statistical 
weight factor (or equally to justify his neglect of such factors). In I we stated 
with care that we were quoting (2/+1) | |? because that was what we 
knew about from experiment, but now the experimental situation is much 
better and we can quote | M |?. We believe that this is the most sensible 
quantity to tabulate in our case and that, in a general tabulation such as this, 
application of any statistical weight factor at any stage after extraction of Ty 
from experiment is meaningless. The Weisskopf unit in terms of which we 
tabulate | 17 f is computed for the transition of one particle in a featureless, 
spinless well from a state of J/=k-+4 to one of J’=}; kis the multipole order 
of the transition and the final state is an s-state. The radiating particle is 
responsible for the entire spin of the ‘nucleus’. Consider now another such 
one-particle transition of the same type, again in a spinless well, but between 
states whose spins are not given by the above relationships (say, for example, 
the final state is not an s-state). Then the statistical weight factor giving the 
speed of this second transition in terms of that of the first (that is to say in 
Weisskopf units) is perfectly definite. For electric transitions for example it is 
(Kennedy and Sharp 1954—whom see for the most complete statement of 
statistical weight factors in both electric and magnetic transitions) : 

rae 
(QF EAD, +Y Wyss’ ; He [Cee | 


J, and 1, are the initial and final orbital angular momentum quantum numbers. 
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state or that at 2-9 Mev, or whether it goes to both. It has here been 
assumed that the transition is to the ground state only. Any other 
possibility will introduce a change of less than a factor of two in the value 
of |. M ? deduced from this datum. 

In computing the M1 radiative width of the 8-89 Mev level in 1°B an 
ambiguity arises because of the two alternative values for the competing 
alpha-particle width ; however these values differ only by a factor two ; we 
have used the mean of the two values thereby introducing an error of 
about 40%. 

In certain transitions from states in 1°O formed in 'N (p,y)1°O we have 
assumed that the greater part of the radiative width represents transitions 
to the ground state. This is based on the report that at these resonances 
the ground state transition is ‘ relatively intense ’, although detailed decay 
schemes have not been published. 

A word is necessary about the choice J=1— for the 9-49 mev level of 
14N since this is not one of the possibilities appearing in the references 
from which other properties of the transition (to the 5-10 Mey level) are 
drawn. As mentioned above we believe this to be an El transition. The 
level at 5-10 Mev decays in fairly equal strength to the ground state (1+-) 
and first excited state (0+, 77=1) of 14N. This seems to fix it as 1+ since 
1— would mean that one or other of these two equal-strength branches 
would be forbidden by the isotopic spin rule. This in turn implies 
0; 1, or 2 for the 9-49 Mev state. Zero is eliminated because of the 


This is the factor that was given explicitly by Moszkowski (1953) for dipole 
and quadrupole transitions ; Hughes (1954) has given it for k=1. It in fact 
depends only on J and J’ and perhaps for this reason the impression has 
sometimes been given that it is a meaningful factor to apply to a real nuclear 
transition between states of spin J and J’ irrespective of how these spins are 
made up. However this factor is significant only if we believe there is some 
sense in supposing that the nuclear spin in both initial and final states is due 
to one particle only—to the one that is making the radiating jump. This may 
sometimes be true in low-lying isomeric transitions, but it obviously cannot 
be true for example in any nucleus of even A or for a transition following the 
addition of a slow neutron to a nucleus of non-zero spin. And it obviously 
cannot be true either for the great majority of the light nucleus transitions 
with which we deal in this survey since so many of the states are produced 
by rearrangement within a configuration. As soon as we begin to deal with 
transitions between states whose spins are made by coupling the spins of 
two or more nucleons together, as is the case for nearly all our light nucleus 
transitions and indeed for nearly all nuclear transitions except very special 
cases, the situation is much more complicated and no simple statistical weight 
factor depending only on J and J’ can be specified. Nor does the above 
expression factor out of the complete statistical weight expression except in 
rare and particular cases. We therefore see that it would be in general erroneous 
and confusing to apply to the experimental Py any statistical weight factor 
depending only on the spins of the initial and final states and on the multipolarity 
of the transition. 

The application of a statistical factor depending on /, and l,, the orbits 
between which the radiating particle jumps (which in general have no direct 
relationship with J and J’), is a different matter. It is certainly meaningful 
to express the experimental radiative width in terms of that expected for a 
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observed anisotropy of the gamma-radiation in the reaction 8C (pea tNe 
2 gives a correlation in the wrong sense; 1 gives a correlation in accord 


with experiment for pure d-waves but with the possibility of a large s-wave 
admixture. 


A further approximate procedure has been resorted to in computing 
radiative widths from radiative capture cross sections for thermal neutrons. 
Such cross sections are known for several light elements and in many cases 
the spectrum of the ensuring gamma-radiation has been measured. Since 
we are here certain of the parity of the radiating state we may usually be 
sure of the type of the transition if the spin and parity of the final state are 
known. It is then possible to compute the associated radiative widths 
using methods outlined earlier (Wilkinson 1953 b). These estimates may 
well be in error by a factor of three or so. Where the target nucleus has 
non-zero spin and where both states corresponding to the two initial 
spin possibilities may radiate to the final state, it has been assumed that 
they both do so in computing the radiative width (which has been taken 
to be the same for both possible transitions, but has been listed only once). 
The maximum error introduced here is a factor of three. One case has 
been omitted from these analyses of thermal neutron cross sections, 
namely 7Li (n, y) 8Li, which is known (Thomas 1951) to proceed chiefly 
through extra-nuclear contributions to the El matrix element. 

The only material admitted to this survey concerns direct measurements 
of radiative widths or lifetimes of transitions supposed to be largely pure 


particle making the same basic orbital jump but in a featureless well. To do 
this we should have to modify the Weisskopf unit (which has /,=0) by what 
could be called a statistical weight factor, but which would depend only on /) 
and /,, and not at all on J and J’. It is in terms of such a modified unit that 
Lane and Radicati (1954), for example, express their theoretical radiative 
widths, and it is, of course, essential that in attempted one-to-one corre- 
‘ spondences this modification be made. But as a rule we do not know /, and /, 
and in any case to apply such a factor would to some degree implicitly prejudge 
the issue of whether the shell model had any relevance for the transition. 
Even in cases such as transitions in our Ip-shell nuclei where we may be fairly 
sure in some cases that /,—1, it is still not an easy matter to settle J). Take 
as an example a compound nucleus state formed by adding rather slow protons 
to a target nucleus. By angular distribution measurements we may find that 
these added protons are almost entirely s-wave. We are therefore tempted 
to say /,=0. But it may be that the actual compound nuclear wave-function 
thereby set up consists predominantly of the target nucleus plus d-wave 
protons with very little s-wave; the small s-wave admixture however pre- 
dominates in setting up the state by addition of protons from outside because 
the centrifugal barrier cuts down the d-wave so strongly. In this case we 
should really use J;=2, not J,=0. Such considerations make it very difficult 
to know what basic orbital jump is taking place even if we wish at the outset 
to commit ourselves to a shell model viewpoint. This makes the application 
of a statistical weight factor of the orbital type too hazardous a procedure 
except in especially clear cases. 
These considerations lead to the view that in all but very special cases, and. 
certainly in the present case, we should use simply | (r in terms of 
unmodified Weisskopf units and not apply any statistical weight factor. 
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in multipole character. For a few transitions angular correlation measure- 
ments have suggested small admixtures of higher multipoles and so from 
them an estimate may be made of the partial gamma-ray width appropriate 
to the higher multipole. However such interpretations are usually 
ambiguous and in any case the operation of extracting the mixing coeffici- 
ent is one of notorious delicacy, so this small body of data has been 
systematically ignored. 

As well as actual widths, the following tabulations include limits 
to widths where these are of interest. It is felt that much valuable 
information could be obtained from records of limits to widths and 
branching ratios even though no actual determination may have been 
made. 

In the following tables the transitions have been arranged under the 
various types, within which the sequence is that of increasing mass of the 
radiating nucleus. #, is the excitation energy in Mev of the level from 
which radiation is emitted, and J is its assumed spin ;. #; is the excitation 
in Mev of the level to which the transition takes place, and H,, is the energy 
of the gamma-ray. J’,is the radiative width in ev and | M [? is this width 
divided by the Weisskopf unit I’,,, appropriate to the type of transition. 
The meaning of Status has already been explained. The references have 
been kept as few as possible and the guiding principle has been to quote 
only the best determination or determinations of the radiative width and 
to trust that those interested in the history of a particular gamma-ray will 
be able to work backwards from the quoted reference. An exception is 
made where the best value is believed to come from unpublished work ; 
here a reference to earlier published work is also given. Where the 
quoted width depends on more than one measurement, for example a 
measurement of total gamma-ray width by one source and of branching 
ratios by another or a determination of thermal neutron capture cross 
section by one group and of emitted gamma-ray spectrum by another, 
both measurements are given as reference. In general no reference is 
given for the spins and parities of the states concerned or for the type of 
transition. This information can in most cases be easily found in the 
review article of Ajzenberg and Lauritsen (1955), whose assessment of the 
reliability of particular assignments has by-and-large been followed, 
although it has sometimes been found possible to award Status A even 
though spins have not been determined with certainty. As explained 
earlier, in those few cases where assignment as to type has been made on 
the basis of internal evidence, Status B only has been awarded. The only 
exceptions made to the rule of not giving references for spins and parities 
are when the references are not self-evident from Ajzenberg and Lauritsen. 
A few difficult cases have already been discussed. 

In the cases of the reactions #C (e, e’) ?C, 160 (y, n) 15O and 0 (y, p)eN 
the experimentally-observed transition is ‘upwards’ from the lower to 


the upper state; the quoted I’, is that inferred for the ‘downward’ 
transition. 
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§ 3. Discussion oF Drpote TRANSITIONS 


We will discuss first of all the E1 and M1 transitions and then the higher 
multipoles. A glance at the tables shows that ||? displays wide 
dispersion. This is seen graphically in the following figures. 

Figure 1 shows the ||? values for E1 transitions plotted against the 
atomic weight of the radiating nucleus. The closed circles indicate 
transitions forbidden by the isotopic spin selection rule for El transitions 
(AT=-+1 when 7,=0): Radicati (1952), Gell-Mann and Telegdi (1953) ; 
the open circles represent transitions that do not suffer this inhibition. 


Fig. 1 
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Matrix elements for El transitions. Transitions forbidden by the isotopic 
spin selection rule in self-conjugate nuclei are distinguished. In this 
figure and in figs. 3 and 6 points have been prevented from overlapping 
by small upward or downward displacements by their diameter or less. 
The exact | M |? values may be taken from the tables. 


The effectiveness of this isotopic spin rule has been established in a number 
of cases additional to those displayed in fig. 1 in experiments where only 
comparisons are made between transitions of different type and so which 
are not considered in this survey which concerns itself only with absolute 
measurements (see Wilkinson 1953 ¢ and MacDonald 1955 for comment) ; 
but it is also clear from the limited data of fig. 1 that such a heehee Se 
is very marked. The forbidden transitions of high |[M|? value in 1B 


138 D. H. Wilkinson on the 


Electric Dipole 


ee SS eee ee 


Reaction E; Ey By J Ty | M |? Status Refs. 
®*He(d, y)sLi | 16-80]0 | 16-80] 3 ll 7-3x10-8 | A* (a) 
SLi(p, y)7Be | 6-35 | 0 6-35 | 4,3 0-7 7-0x10-3 | B (b) (ce) 
6Li(p, y)7Be | 6-35 | 0-43] 5-92| 4,4 0-4 4-7x10-3 | B (b) (c) 
*Be(n, y)!Be | 6-81 | 0 6-81| 1 7-0 0-043 A (a) (e) 
*Be(n, y)’Be | 6-81 | 3-37| 3-44] 1,2 0-9 0-043 A (a) (e) 
SLi(a, y)B | 5-11] 0 Bll] 2 <0-012 | <1:8x10-*| B* (f) 
®Li(z, y)B | 5-11] 0-72] 4:39] 2 <0-012 | <2:8x10-+]| B* (f) 
6Li(a, y)B | 5-11|2-15| 2-96] 2 <0-012. | <9-1x10-4| B* (f) 
*Be(p, y)B | 6-89]0-72] 617| 1 0-7 5-9x10-*| B (g) (h) 
*Be(p, y)’B | 6:89] 1-74] 5:15] 1 3-8 0-055 B (g) (bh) 
*Be(p, YB | 6-89}215| 4:74] 1 0-7 0-013 B (g) (h) 
*Be(p, y)°B | 6-89] 5-16] 1-73] 1 0-7 0-27 B (g) (h) 
*Be(p, y)B | 7-48 | 0 7-48| 2 21 0-10 A* (i) 
*Be(p, y)°B | 7-48]0-72| 6-76] 2 <2 <0-013 A (i) 
Lila, y)4B | 9-28] 0 9-28| $ 0-8 1-9x10-3 | B (j) (k) 
7Li(a, y)4B | 9-28] 4-46] 4-82] 3 1-9 0-031 B (j) (is) 
Lila, y)4B | 9-28] 6-81| 2-47]. 3 0-6 0-073 B (j) (ix) 
1B(p, y)2C | 16-57] 4-44] 12-13] 2 15 0-015 A* (1) 
UB(p, y)2C |17-22}0 |17-22| 1 40 0-014 B* (1) 
UB(p, y)2C | 17-22] 4-44] 12-78] 1 20 0-017 B* (1) 
LO(n, yO | 4-95 | 0 4-95] ¥ 1:3 0-018 A (e) (m) 

1 22C(n, y)#C | 4-95] 3-68] 1-27] 4 0-6 0-48 B (e) (m) 
2C(p, y)3N | 2-37 | 0 2-37) 4 0-67 0-083 A* (n) 
2C(p, y)3N | 3-51 | 2-37| 1-14] 3 0-04 0-044 A (0) 
BC(n, y)MC | 8-17] 0 S17) 1 0-5 1:5x10-2 | A (p) 
BO(p, y)4N | 8-06 | 0 8-06 | 1 10-2 0-031 A* (q) (r) 
8Q(p, y)4N | 8-06] 2-31] 5-75| 1 5x 10-2 4-1x10-4 | A (s) (t) 
BC(p, y4N | 8-06]3-95| 4-11| 1 13 0-029 A* (q) (x) 
18C(p, y)4N | 8-62] 6-23] 2-39| 0 1-4 0-16 B* (q) (r) 
18C(p, y)4N | 8-70] 0 8-70 | 0 46 0-11 A* (q) (r) 
BC(p, y4N | 8-70| 5-69] 3-01] 0 5 0-28 B (q) (r) 

BO(D, y4N | 9-18] 0 9-18] 2 10-6 0-022 A* (q) (a) 
18C(p, yMN | 9-18] 6-44] 2-74] 2 1-18 0-091 B* (q) (x) 
18C(p, y4N | 9-49] 5-10| 4:39] 1 8-2 0-15 B (q) (r) 
UN(n, y)8N |10-83]}0 | 10-83] 4,3 4-6 54x 10-3 A (a) (u) 
UN(n, y)§N | 10-83 | 6-33 | 4:50] 42 4-0 0-067 B (a) (u) 
U4C(p, y)8N | 11-43 | 0 12-48 + 2-3 2-3 x 10-3 A* (v) 
UQ(p, y)5N | 11-57|0 | 11-57] 4 28 0-028 Ax (v) 
“N(p, yO | 8-790 | 8-79] 4 0-48 1-1x10-8 | A (w) 
UN(p, yO | 9-77] 0 9:77 | 4.3 92 0-15 A ae 
16Q)* lifetime | 7-12 | 0 712} 1 | >0-058 | >2-3x10-4] A =) 
N(p, yO | 13-050 |13-05| 1 150 0-098 A y) 
BN(p, y)°O |13-05| 6-05] 7-00| 1 < 0-2 <8-4x10-4| A C 
LO(y, n)O | 2033/0 | 20-38] 1 43 7-4x10-3 | B* 2 
16 15 o).KQ = ase (at) 

Oly, n)80O | 20-5810 | 20-58] 1 54 8-9x10-3 | B* 1 
WO(y,n)*O | 20-79 | 0 20-79 1 58 9-3 10-3 B* io) 
’O(y,n)8O |20-93/0 | 20-93] 1 83 0-013 B* (ar) 
w’O(y,n)8O |21-21]0 |21-21] 1 143 0-022 B* ie 
WO(y,n)6O | 21-52}0 |21-52] 1 98 0-015 B* (at) 
WO(y, n)%O | 22-3710 | 22-37] 1 322 0-042 Ax ta) 
14O(y,n)O | 22-54 | 0 22-54 1 480 0-059 Ax ta) 
1O(y,n)®O | 22-76) 0 22°76" ea 549 0-066 A* inf 
W’O(y, n)9O | 23-02 | 0 23:02] 1 501 0-059 ‘Ax (a) 
WO(p, y)7F | 3-86] 0 3-86| 3 0-6 0-015 A fe a 
saith lifetime Ce WEG 0-11 4 6:7 10-7 6-4 10-4 ‘Ax ap ae ) 
19H (p, y)?2Ne | 13-70 | 1-63 | 12-07 2 <2-8 = 1Gise 10-3 ‘AX Ca 
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Magnetic Dipole 


Reaction E; Ey Ey J Ty | M ? Status Refs. 

7Li* lifetime | 0-48 | 0 0-48) 4 | 46x10 2-0 A* (£1) 

| 7Be* lifetime | 0-43 | 0 0-43) 4 3:3 x 10-8 2-0 A* (£2) 

| 7Li(p, y)®Be | 17-63 | 0 17-63| 1 16-7 0-14 A* | (gi) (bh) 
7Li(p, y)®Be | 17-63 | 2-90] 14-73] 1 8:3 0-13 A* | (gt) (ht) 
*Li(p, y)®Be | 17-63] 4-2 | 13-4 | 1 0-45 9:0 x 10-3 B (i2) 
7Li(p, y)®Be | 17-63| 5-4 | 12-2 | 1 0-43 0-011 B (it) 
7Li(p, y)®Be | 17-63] 7-55 | 10-08] 1 0-25 0-012 B (i2) 
7Li(p, y)®Be | 18-14] 0 18-14] 1 53 0-042 B (j*) (Ist) 
SLi(a, y)B | 4:77) 0-72| 4-05] 1 0-15 0-11 B (f) 
Li(x, y)B | 5-16] 0 5:16| 2 0-03 0-010 B (f) 
SLi(a, y)"B | 5-16]0-72| 4-44] 2 0-15 0-082 B (f) 
SLi(a, y)"B | 5-16]2-15| 3-01} 2 0-42 0-74 B (£) 
*Be(p, y)’B | 7-56] 0-72| 6-84] 0 4:8 0-72 AY (i) 
*Be(p, y)!B | 7-56] 2-15| 5-41] 0 1-2 0-36 A* (i) 
*Be(p, y)°B | 8-89] 0-72] 8-17] 2 8 0-75 At an) 
7Li(x, y)4B | 9-19] 0 9-19] $ <0-010 | <6-2x10-*| B (j) (x) 
7Li(a, y)4B | 9:19] 4:46] 4-73] 3 0-4 0-18 B (j) (x) 
"Li(a, y)"@B | 9-19|6-81| 2-38] 3 0-07 0-25 B (j) (x) 
WB(p, yC | 9-701 0 9-70| 3 2-3 0-12 B (m?) (n1) 
UB(p, y)2C | 16-10 | 4-44] 11-66] 2 70 2-1 A (1) 
BC(p, y)E=N | 3-51] 0 3°61] 2 0-67 0-74. A* (n) 
18C(p, y)4N | 8-62] 0 8-62} 0 0-73 0-055 A* (q) (x) 
BC(p, y)#N | 8-62|3-95| 4-67] 0 3:1 1-48 A* (q) (x) 
C(p, y4N | 8-98] 0 8-98] 1 0-17 0-011 A* (q) (x) 
MC(p, y)*N | 10-70] 0 10-70) |, 3 0-12 4:7x10- A* @v) . 
U4C(p, y)=N | 10-81] 0 10-81| ¢ 0-01 3-8 x 10-4 A (v) 
MN(n, y)8N | 10:83 | 7:36| 3-47] 4,3 2:5 2-8 B (d) (u) 
M4N(n, y)=N | 10-83] 8-28 | 2-55] 4,3 0-9 2-6 B (d) (u) 
MN (n, y)5N | 10-83] 9-16} 1-67] 4,3 0-4 4-0 B (d) (u) 
MC(p, y)5N | 11-30 | 0 11-30] 4 0-25 8-4 10-3 A* (v) 
VN(p, y)8O | 8:34] 0 8-34| 4,3 1:3 0-11 A (w) 
UN(p, yO | 8-98} 0 8-98 | 4 0-30 0-020 B (w) 
UN(p, y)*O | 9-04] 0 9-04] 3,3 1-04 0-067 B (w) 
19F(n, y)2°F | 6-60] 0 6-60 | 0,1 0:2 0-033 A (0) (p*) 
LF (p, y)2°Ne | 13-19 | 0 13:19) 1 <0-7 <0-015 A (q?) 
L9E(p, y)2°Ne | 13-19 | 1-63} 11-56] 1 <0-9 <0-028 A (q4) 
19F(p, y)2°Ne | 13-51 | 0 13-51] 1 <0-02 <3-9x10-4| A (q}) 
19 (p, y)2°Ne | 13-51 | 1-63} 11-88] 1 2:3 0-065 A* | (st) (64) 
19F'(p, y)2°Ne | 13-76 | 0 13; 76a ak <2-1 < 0-039 A (q?) (t*) 
19F(p, y)2°Ne | 13-76 | 1-63 | 12-13 | 1 <3-5 < 0-094 A (q+) (t4) 

Electric Quadrupole and Octupole 
es hep ee Se ee 
Reaction E, |Hy| By | J Ty Type | M |? Status] Refs. 

10B* lifetime | 0-72] 0 | 0-72] 1] 1-0x10-% |. H2 2-0 
12C(e,e’)2C | 4:44] 0] 4:44) 2] 1-0x10-? | E2 |] 1-7 
22C(e, e’)12C 9-61|70°} 9-61 | 2 | 20x10 |" B2 } 7x10-* 
UB(p, y)2C | 16-10} 0 | 16-10} 2 3 E2 0-84 
16Q* lifetime | 6:14] 0} 6-14] 3 | 6-0x10-° | E38 7:3 
16Q* lifetime | 6-91] 0] 6-91] 2 | >0-038 E2 | >0-50 
O(y, p)ON | 14-7 | 0 | 14-7 | 2 46 E2 14 
170* lifetime | 0-87] 0 | 0:87] 4 | 2-7x10-° | KH2 1-0 
19F* lifetime | 0-20] 0| 0-20] 3 | 4:6x10-° | E2 2-4 


a ell ee ee ee 


140 D. H. Wilkinson on the 


TABLE REFERENCES 


(a) Buatr, J. M., Huvrz, N. M., and van Parrer, D. M., 1954, Phys. Rev., 96, 1023. 

(b) Erpman, K. L., Warren, J. B., James, D. B., and ALEXANDER, T., 1954, Phys. 
Rev., 96, 858. 

(c)  Basuxty, 8., and Cartson, R. R., 1955, Phys. Rev., 97, 1245. 

(d) Hueuss, D. J., e¢ al., 1952, AHC Unclassified Report No. 2040. 

(e) BarrHotomew, G. A., and Krysry, B. B., 1953, Canad. J. Phys. 31, 49. 

(f) Jonzs, G. A., and Witxrson, D. H., 1954, Phil. Mag., 45, 703. 

(g) Cartson, R. R., and Netson, E. B., 1955, Phys. Rev., 98, 1310. 

(h) COxeac, A. B., and Wiixrnson, D. H., unpublished work. 

(i) Hornyak, W. F., and Coor, T., 1953, Phys. Rev., 92, 675. 

(j) BENNETT, W. E., Roys, P. A., and Torrent, B. J., 1951, Phys. Rev., 82, 20. 

(ete Jones, G. A., and Witkinson, D. H., 1952, Phys. Rev., 88, 423. 

(1) Beckman, O., Huvus, T., and Zupanéié, C., 1953, Phys. Rev., 91, 606. 

(m) Way, K., Fano, L., Scott, M. R., and THEw, K., 1950, National Bureau of Standards 
Circular 499. 

n) SHaGRAVE, J. D., 1951, Phys. Rev., 84, 1219. 

0) Woopsury, H. H., Tortestrup, A. V., and Day, R. B., 1954, Phys. Rev., 93, 1311. 

p) Hennie, G. R., 1954, Phys. Rev., 95, 92. 

q) SHacrave, J. D., 1952, Phys. Rev., 85, 197. 

r) Woopsury, H. H., Day, R. B., and Tottesrrup, A. V., 1953, Phys. Rev., 92, 1199. 

s)  Czece, A. B., and Wizxrnson, D. H., 1953, Phil. Mag., 44, 1269. ; 

t)  Hrrp, B., Warrrnean, C., Burter, J., and Coxiim, C. H., 1954, Phys. Rev., 96, 702. 

u) Kinsey, B. B., BartHotomew, G. A., and WALKER, W. H., 1951, Canad. J. Phys., 
29, 1. 

(v) BarrHotomew, G. A., Brown, F., Govz, H. E., Lirgzrianp, A. E., and Pau, 
E. B., 1954, Phys. Rev., 96, 1154. 

(w) Duncan, D. B., and Perry, J. E., 1951, Phys. Rev., 82, 809; see also Li, C. W., 
1953, Phys. Rev., 92, 1084. 

(x) Drvons, 8., MANNING, G., and Bunsury, D. St. P., 1955, Proc. Phys. Soc. A, 68, 18. 

(y) ScHarpt, A., Fowier, W. A., and Lauritsen, C. C., 1952, Phys. Rev., 86, 527. 

(z GoLpRING, G., 1954, Proc. Phys. Soc. A, 67, 930. 

(a!) Prnroup, A. 8., and Spicer, B. M., 1955 (in course of publication). 

(b!) Lauspenstetn, R. A., LAUBENSTEIN, M. J. W., Koxuster, L. J., and Mostey, R. C., 
1951, Phys. Rev., 84, 12. 

(c1) DuBrrpex, L. A., Barnzs, 8. W., Buck, J. H., and Strain, C. V., 1938, Phys. Rev., 
58, 447. 

(d+) Turtrion, J., Barnes, C. A., and Lauritsen, C. C., 1954, Phys. Rev., 94, 1076 (un- 
corrected for internal conversion). 

(e1) Wrixinson, D. H., and Cuzae, A. B., 1953, Phil. Mag., 44, 1322. 

(f1) Drvons, S. (unpublished)—see AgzENBERG, F., and LauriTsEen, T., 1955; also 

Exuiort, L. G., and Betz, R. E., 1949, Phys. Rev., 76, 168 for Li’*. 

1) Fow.er, W. A., and Lauritsen, C. C., 1949, Phys. Rev., 76, 314. 

) Watxer, R. L., and McDantet, B. D., 1948, Phys. Rev., 74, 315. 

it) Inari, E. K., and Boyuz, A. J. F., 1953, Phil. Mag., 44, 1081. 

j1) Kraus, A. A., 1954, Phys. Rev., 93, 1308. 

kt) Lrsrrman, D. A., 1955—see § XIV of AszENBERG, F., and Lauritsen, T. (1955). 

14) Mackry, R. J., 1954, Phys. Rev., 94, 648, and unpublished work. 

m1) Watxer, R. L., 1950, Phys. Rev., 79, 172. 

nt) Brown, A. B., SvypER, C. W., FowLer, W. A., and Lauritsen, C. C., 1951, Phys. 
Rev., 82, 159. 

(01) Seren, L., FRreEpLANDER, H. N., and TurKEL, 8S. H., 1947, Phys. Rev., 72, 888. 

(p?) KINSEY, B. B., BartHotomew, G. A., and WaLKER, W. H., 1951, Phys. Rev., 83, 519. 

a Cieaa, A. B., and Wiixrnson, D. H. (unpublished)—see also (e!). 

( 

( 


r+) Carver, J. H., and Witxrnson, D. H., 1951, Proc. Phys. Soc. A, 64, 199. 

s') Circe, A. B., Jonns, G. A., and Witkrinson, D. H., 1955, Proc. Phys. Soc. A, 68, 538. 

t1) poe corer Givin, H. H., Kern, B. D., and Haun, T. M., 1955, Phys. Rev., 
5 TPA, 

(u') Turion, J., and TeLeeptr, V. L., 1953, Phys. Rev., 92, 1253. 

(v') Horsrapter, R., and RaveNHALL, D. G. (private communication). 

(w1) Spicer, B. M., 1955, Phys. Rev., 99, 33. 

(xt) Wixinson, D. H., 1955, Phys. Rev., 99, 1347. 

(y') Jonus, G. A., Puitiies, W. R., Jonnson, C. M. P., and Witxinson, D. H., 1954 

Phys. Rev., 98, 547 (uncorrected for internal conversion) ; r 


Radiative Transitions in Light Elements: II 141 


are anomalously strong and probably represent a very high admixture of 
lin the predominantly 70 state at 6-89 Mev, although they only 
have Status B and we cannot yet be sure of them. Inthe aise discussion 
we set aside such isotopic spin forbidden transitions and consider only the 
‘ordinary’ ones. Figure 2 shows the ‘ ordinary ’ E1 transitions grouped 
according to their |M|? value. The ordinate is the actual number of 
transitions in the given range of |M|?; 46 transitions are represented ; 
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Distribution of E1 matrix elements. Transitions forbidden by the isotopic spin 
rule are omitted and transitions for whose strength we have only limits 
are not shown. 


those for which only limits are available are not shown. This distribution 
is a rather tidy one ; it centres on |M |?=0-032, and a variation by a factor 
of seven in speed on either side of this value covers 85° of the transitions.* 
This indicates the degree of reliability of an offhand estimate of the width 
of an E1 transition in a light nucleus and suggests that the best guess will 
be I,~0-032F,,,. At the same time we must bear in mind that rather 
weak transitions do occur, sometimes from the same level as a strongone 

of the same type: a good example is given by the transitions from the 
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* We must ask whether the rapid fall of the distribution on the side of low 
| 12|? is perhaps artificial and really due to the missing of weak transitions, 
not to their absence. While it is felt that this is largely a genuine fall we must 
certainly be ready to find many more weak transitions, and again it should 
be stressed that there is a lot of interest in looking for and comparing all 
transitions of a given type that can occur from a given radiating level, and at 
any rate placing and publishing limits on their strengths when the actual 
transition cannot be detected. This is particularly important for the omnibus 
comparison of the experimental distributions with the expectations of shell 
theory—see later. 
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1— state of 160 at 13-05 Mev to the 0+ ground state and to the 0+ state 
at 6-05 Mev; the former has a ||? value of 0-098, while the latter has 
| MP<8-4x 10-4. There seems to be a tendency for the very weak E1 
transitions to involve states that are rather unusual in other respects such 
as the low-lying odd parity state in 19F (}— at 110 kev) and the pair state of 
16 (0+ at 6-05 Mev). It may be that these weak transitions merely 
emphasize that these states are indeed ‘ curious’ and not very closely 
related to their neighbours so that they sometimes do not combine very 
well with them. A similar impression is gained from the few low-lying E1 


Fig. 3 
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Matrix elements for M1 transitions (see also caption of fig? 1): 
transitions established and measured in heavy nuclei (Goldhaber and 
Sunyar 1955) which all have very low| M |? values. In all these cases in 
heavy nuclei the low-lying state with the change of parity is a ‘ curious’ 
state that has interloped into the expected level structure. 

Figure 3 shows the | M |? values for M1 transitions in the same Way as 
does fig. 1 for E1 transitions. Again we see a considerable spread in 
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| iF values, and this is displayed in the histogram of fig. 4. This distri- 
bution is to be compared with that of fig. 2 for E1 transitions : it contains 
34 transitions. The M1 distribution differs in two ways from the EL: 
it is a wider, more rambling distribution—a factor 20 in speed is now 
needed on either side of the mean before we include 85° ,, of the transitions : 
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Distribution of M1 matrix elements. Limits are not shown. 


secondly it centres on |M |?=0-15, which is a factor of 4-7 greater than 
the corresponding centring point for the E1 distribution. We may now 
recall the remark of § 1 that, for shell model transitions, we might expect 
| M |? values about a factor of four greater for M1 than for E1 transitions 
because of the diminution in the E1 Weisskopf unit due to the effective 
charge. This expectation seems to be realized. 

The greater spread in the M1 distribution can also perhaps be understood 
qualitatively in shell model terms. In the first place M1 transitions are 
somewhat more complicated than E1 since they contain both space and 
spin contributions and this gives the possibility of greater fluctuations. 
Secondly M1 transitions take place largely by rearrangements within a 
given configuration—the Ip shell in the majority of our examples—and 
this means both that much larger values of the quantity n? times the 
parentage overlap (Lane and Wilkinson 1955) are on occasion possible 
for them than for E1 transitions, and also that, because of the multiplicity 
of common parents often associated with transitions within a configuration, 
small values of the parentage overlap due to cancellation of contributions 
from the various parents will also tend to occur rather frequently. We 
contrast this situation with that obtaining for El transitions where. 
because of the change of configuration necessarily involved for the radiating 
particle, the largest value that can be taken by the square of the product 
of the fractional parentage coefficients of the initial and final states for 
any common parent is n+. This fact must be taken together with the 
observation that 2s and 1d nucleon interactions with Ip-shell nuclei tend 
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to display reduced widths of the order of the single-particle value (Lane 
1954) which shows that there will often be a unique or almost unique 
parentage for the initial state. These two remarks imply that we cannot 
get such big enhancements for El as for M1 transitions and also, since 
the parentage is often almost unique, that we do not have so much 
possibility of weak transitions due to cancellation of contributions from 
several parents. All these effects then suggest that the larger spread 
found for the M1 strengths than for the El strengths is in qualitative 
accord with shell model expectations. 

We may now ask whether in fact the shell model can give an account 
of the actual numerical strengths and spread in strengths that we have 
displayed in figs. 2 and 4. We have discussed in § 1 how we should in 
fact expect, within the pure shell model, a rather considerable straggle in 
the theoretical ||? distribution. The situation in respect of M1 
transitions is really too complicated to do anything about statistically 
because of the great multiplicity of parents but, as we have just noted, 
for the El transitions we can assume as a very crude approximation, and 
neglecting all configuration interaction, that the parentage is unique. 
This effects an enormous simplification. Consider now an E1 transition 
from an initial state of quantum numbers L, S, J, 7’ consisting of a 2s or 
1d nucleon (the /, orbit) plus n— 1 1p nucleons into a final state of quantum 
numbers L’, S’, J’, 7’ with all n nucleons in the Ip shell. (All but two of 
our experimental E1 transitions are in lp-shell nuclei.) We now make 
our assumption that the initial state has some unique parent in Lp”~!; 
the final state has fractional parentage coefficient <|) for this parent 
state whose quantum numbers are L,, etc. Then (Lane and Radicati 
1954): 


> on LT Ta 
M\?= —sU(11L0,;1,L')UOL' IS; LS’ U(14s es he J 
| rt ( p? "0 ) ( )<}> (1577 eee Pik ins 
where 


U(abed ; ef)=V (2e+1)(2f+1)W (abcd ; ef). 


Note that [7 |? in this theoretical formula is not defined in quite the 
same way as we use throughout this survey in connection with eaperimental 
transitions. For the experimental transitions ||? always means the 
width in the appropriate Weisskopf units ; in the above formula it means 
the width in terms of that of the appropriate single-particle transition 
(2s+1p or ld-+1p). There is little difference between the two definitions 
of |M/|? for these 1p-shell transitions, and it is not here worthwhile to 
apply a correction factor before comparing the experimental results with 
the theoretical. In an attempted one-to-one correspondence between a 
particular transition and theory the correction factor to allow for the 
actual orbits involved is of course essential as is the use of more realistic 
wave functions than the uniform wave functions of the Weisskopf unit. 

As mentioned earlier we do not now attempt the dangerous task of 
specifying the quantum numbers and parentage relationships of individual 
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experimental transitions, but rather simply ask how the gross assembly 
of experimental E11 transitions summarized in fig. 2 compares with the 
gross assembly of theoretical E1 shell-model transitions computed for 
the same region of the periodic table and between states whose excitations 
it is reasonable to suppose might correspond roughly with those of the 
states actually concerned in the experimental transitions (the review of 
Inglis (1953) is taken as a guide here). In other words we ask in general 
terms if the shell model predicts the same sort of ||? values and spread 
in those values that obtain in practice. In this spirit of general enquiry 
we have therefore used the above expression to compute the theoretical 
| |? values of a great number of E1 transitions (about 250 of them) 
between states satisfying the above criteria.* These ||? were reduced 
to a histogram as were the experimental values ; the smoothed theoretical 
distribution is shown in fig. 5. It bears a most remarkable resemblance 


Fig. 5 
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The expected distribution of El matrix elements computed in LS coupling 
for lp-shell nuclei and states of ‘reasonable’ excitation under the 
assumption of unique parentage. 


in absolute |J/ |? values and in spread to fig. 2. While we are very willing 
to admit that this detailed resemblance is due largely to chance, we feel 
that we should insist that this computation, or jew d’esprit, does show 
that the pure shell model predicts that the run-of-the-mill El transition. 
in the Ip shell should have a strength much less than a single-particle 
unit, and that a large straggle in strengths should obtain, just as 1s found 
_ in practice. 


* The largest | Mf |? value found was 5/3— for the transition 
bps id (D-s1, laa —, T=1) > 1p?? (L'=0, S’=0, J’=0+, 7” =0). 
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We have used LS coupling in drawing up fig. 5 because it seems from 
other data (e.g. Christy 1953) to be the better approximation of the extreme 
schemes in the Ip-shell ; jj coupling gives a very closely similar distribution 
—and so, we may feel sure, would any intermediate coupling. — 

We may notice in this omnibus comparison that the rapid fall of the 
theoretical shell-model distribution on the side of low [M |? is reproduced 
by experiment. However, as we have stressed earlier, the possible 
importance for the experimental distribution of weak missed E11 transitions, 
although probably not great, is not really known at the moment, and will 
not be known until a conscious effort is made to look for and set limits on 
all transitions of given type that can occur from given levels. If transitions 
of low |M|? do not appear in relatively greater abundance than at the 
moment we may probably construe the agreement between figs. 2 and 5 
on the side of low ||? as some evidence for the absence of much remote 
configuration interaction and similar effects in light nuclei at these 
excitations. 

It is of considerable practical interest to see how far we may base 
determinations of parity change in light nuclei on the strength of gamma- 
ray transitions. That is, we ask whether we can tell apart El from M1 
transitions on the basis of their radiative widths. Now although the El 
Weisskopf unit is some 20 to 30 times greater than the M1 unit in our 
light nuclei, we have seen that |M/?| values are in the mean some four or 
five times lower for El than for M1 transitions and also that each 
distribution shows a considerable spread. We therefore expect a great 
deal of overlapping of strengths for transitions of the two types. This is 
shown in fig. 6 where we display El and M1 transitions together, merely 
taking out the E,? factor to compare ‘reduced’ strengths. The open 
circles are the E1 transitions and the closed circles are the M1 transitions. 
It is seen that there is indeed a great deal of overlap between the two 
distributions and that in general we cannot base a distinction on the 
strength of the transition. It is seen, however, that the very strong 
transitions tend to be E1 and the very weak ones M1. This former, E1, 
tendency is made clearer by a few remarks about the strong M1 transitions. 
The two strong M1 transitions at A=7 are the low-lying mirror transitions 
in ‘Li and 7Be. These transitions take place between closely related 
simple states in an especially favourable manner, and this situation is 
unlikely to occur again ; and we may therefore regard them as a special 
case. The fast transitions at A—=15 are derived from 14N(n, y)15N ; this 
derivation of widths from thérmal neutron capture cross sections is the 
least satisfactory of the methods we have employed and, as we have 
remarked, the error in the deduced widths may be considerable ; in 
addition, these transitions have only Status B. The fast point at 4—12 
is from the low-lying resonance in 4B plus proton, some of whose cross 
sections, as Ajzenberg and Lauritsen remark, derive from work of 
considerable antiquity and uncertainty. So, excluding the special points 
at A=7, there remains only one reliable fast M1 transition above 
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I’, (ev)/E,? (mev)=0-02 as against twenty-odd E1 transitions above the 
same value. It is felt that if any new transition shows a value of 
I, (ev)/E,3 (Mev) greater than 0-02 then there is something like a 10: 1 
chance that it will be El. As mentioned above we have already used 
this criterion for classing as El transitions four gamma-rays in 14N whose 
strengths exceeded this limit by factors of three and more. 


Fig. 6 
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The comparison between El and M1 radiative widths after removal of the 
E,? factor. Isotopic spin forbidden El transitions are omitted as also 
are all limits and an M1 point at 4 =15 with I, (ev)/H,? (mev)=8 x 10° 
(see also caption to fig. 1). 


Because of the possibility of yet-to-be-discovered long, low-strength 
tails to both El and M1 distributions, we feel that it would be unsafe to 
base any present recommendation on the apparent predominance of M1 
transitions among the weakest gamma-rays. _ 

It was mentioned in § 1 that, when configuration interaction and more 
complicated effects set in strongly, we must expect a dependence of 
transition strength on level density to appear, and that such an effect is 


148 D. H. Wilkinson on the 


manifest among the heavy nuclei. It is therefore of interest to look for 
this effect in the light nuclei of the present survey. It does not seem to 
be strongly evident. D, the level spacing in Mev, varies from perhaps 
three to six for low-lying levels of the lightest nuclei down to about 0-2 
for the sequence of 1— levels found in 1%O(y, n)!°O, and yet no parallel 
diminution in ||? occurs. Even if exception is taken to this illustration 
on the ground that these 1O states are in the region of the * giant 
resonance ’ and so perhaps not very typical, yet there remain examples 
such as 13C(p, y)!4N where levels are several times as dense as in the 
lightest elements but radiate just as strongly. It may be that this 
observation of the relative insensitivity of ||? to D means that in these 
light nuclei for dipole transitions at the excitations with which we concern 
ourselves, effects such as configuration interaction and multiple-particle 
excitation have not yet become very important, and that we are indeed 
dealing chiefly with states whose shell-model description may still be made 
largely in terms of the change of orbit of one nucleon or rearrangement 
within a pure configuration. We note that this is also the conclusion at 
which we have tentatively arrived above if we are to regard the sharp 
falling-off of the |M |? distribution at low | |? values as genuine. Perhaps 
with 2°Ne we see the beginnings of a dependence of ||? on D. It will be 
of interest for this reason to extend this survey to the region above A= 20 
but at the moment this is not considered very profitable because there is 
there a definite experimental tendency to establish information about the 
strongest transitions, and so a distorted picture would be obtained. 

There is also no clear dependence of | |? on the evenness or oddness 
of NV, Z and A for either the E1 or the M1 transitions. Such a dependence 
is not especially to be expected from the shell-model computation for E1 
transitions or from more general considerations. 


§ 4. Discussion or HigHER MuLtrreoLe TRANSITIONS 


As we may see from the table of E2 and E38 transitions, our knowledge 
of the higher multipole transitions is very restricted. But already it can 
be seen that they tend to be fast, some of them very fast. It might be 
argued that those higher multipole transitions which have inherently low 
width will be essentially chosen by the experimenter for their large |M [? 
value because he will be able to see only the strongest ones. This is a valid 
criticism at this stage, and would preclude any reliable attempt to 
systematize, say, K2 transitions in the way we have done for El and M1 
transitions. It applies to the “B(p, y)!2C and 1O(y, p)°N entries (and 
in any case we should strongly discount the strength of this transition 
in 17C as based on the old and unreliable work referred to in §3; the true 
width may well be an order of magnitude less). But it does not apply 
to the other transitions which we shall now discuss, all of which are based 
on lifetimes measurement which would have been made more easily had 
the widths been less rather than more. So although we may not treat 
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the present sample as representative, we may examine the transitions as 
individuals, many of which are very instructive.* 

The transitions in !°B, from the 4-4 Mev state of Cand from the 6-91 Mev 
state of 'SO are fast but not very remarkably so and need not detain us. 
The rest, however, contain noteworthy examples which seem to indicate 
quite clearly that the shell model is not adequate to explain all aspects of 
radiative transitions in light nuclei; it must be admitted that cooperative 
model shave a place even below A=20. 

Consider first the 14-7 Mev transition in 1°O. If it turns out that this 
is indeed an E2 transition as we suppose it to be, we have a very clear 
contradiction with the shell model which is not capable of finding a |M (2 
value greater than 7/5 for this transition—which it gives for LS coupling 
when we take the upper state as Ip-!1f (L=2, S=0, J=2+,7=0); it has 
been suggested (reference (x!) above) that this may be a one-phonon 
transition in the Bohr—Mottelson sense, the state at 14:7 Mev being the first 
state of surface oscillation of 1*O. The hydrodynamical parameters needed 
to describe this experimental transition agree to within a factor of about 
two with the ‘ classical’ parameters of the Bohr—Mottelson scheme, the 
senses of departure being those naively expected for a double-closed-shell 
nucleus such as 1*O. Because of the high excitation involved and the 
attendant breakdown of the adiabatic approximation it may be that this 
state of surface oscillation is in fact shared among several ‘ ordinary ’ 
states, but most densely around 14:7 Mev. The existence of such a 
powerful cooperative E2 state as this, although we cannot yet regard it 
as definitely established, is of great use in discussing the next of these 
interesting fast E2 transitions, namely that in !7O. This appears to have 
about the right strength of one Weisskopf unit until one notices that the 
‘radiating particle ’ is a neutron and not a proton and so should, on the 
shell model, give a strength of |M |?~(Z/A?)?=10-? rather than unity. 


* We might at this point mention a curious phenomenon in connection with 
the transition from the 7-48 mev 2— state in !B to the ground state for which 
Devons and Goldring (1954) have measured the angular correlation of internally- 
converted pairs. These workers find a pair distribution that is inconsistent 
with a simple El transition but which could be explained if roughly half the 
radiative width were made up of M2 or E38 radiation. If this half were M2 
it would have | M |? ~ 3x 10%, while if it were E3 it would have | M |? ~ 108. 
These figures are so huge and so much in excess even of Z° that we cannot 
accept them at their face value and must suppose that the apparent anomaly 
of the distribution is due to some unsuspected mechanism or to approximations 
in the as yet incompletely established theory. 

+ The parentage situation in 1°B is complex and so shell-model calculations 
are not especially trustworthy. It is, however, difficult to achieve | ba=2 
in any E2 shell-model calculation. The situation in "C for the 4-44 Mev 
transition on the other hand is rather simple in both extreme coupling schemes, 
and we may have more trust in the computation ; we find | M |zs?=14/25 
and | U |,?=1/5. ; "7 ee 

+ Actually for a pure one-particle 2s>Id(/2) transition such as this the 
shell model leads us to expect, for a radiating proton, | M [’=38 (see the 
footnote on statistical weights in § 2). But this factor is of little use against 
the discrepancy of x 10°. 
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This transition is therefore an anathema maranatha for the shell model. 
As Thirion and Telegdi (reference (u1) above) have remarked, the way 
out of this dilemma that does least violence to our hope that 170 should 
a rather good shell-model nucleus (ef. its small ground-state quadrupole 
moment) is to suppose that there exists somewhere in '%O a 2+ state 
possessed of a very fast E2 transition to ground which can act as an 
auxiliary parent for the 170 states and which is responsible for the 'O 
transition. (In this inverted sense we should rather regard the odd neutron 
as the parent of both the 17O states, the E2 emission being effected by the 
one-phonon transitions of the attached 160.) We now see that the above- 
discussed fast E2 transition from the 14-7 Mev state of 1°O is ideal for this 
role of ‘radiating core’ in 170; an admixture of about 7°% in intensity 
of this 1®O state in the 0-87 Mev state of 170 suffices to explain the observed 
transition. This does not seem quite unreasonable. The same deus ex 
machina may suffice to explain the too-fast E2 transition in !%F where 
some such help is certainly needed since addition of neutrons to a proton 
can never speed an E2 shell-model transition (contrast the situation for 
E1 transitions).* , 

The last of these fast transitions is the E3 in 1®O. Again the shell model 
cannot account for so great an enhancement; the largest shell-model 
for |M|? is found in LS coupling and is 5/7 when the upper state is 
lp-tld (L=3, S=0, J=3—, T'=0) (unless we are prepared to go into 
the g-shell where an enhancement of 9/7 can be found). The increase in 
nuclear radius of 30°, to 40°, needed to reconcile this transition with 
shell-model expectation seems rather implausible and it seems likely that 
we are again witnessing some sort of cooperative phenomenon. © This 
time it does not seem probable that this is itself a state of surface oscillation 
since the E3 one-phonon mode is expected at a considerably higher energy 


* An equivalent way of speaking of these too-fast transitions in 170 and 
19F is to say that the motion of the external particles induces a quadrupole 
distortion of the 1®O core, and, since this distortion in the Bohr—Mottelson 
weak coupling approximation is determined by the number of the external 
particles and not by their charge, the result is to add a constant effective 
additional charge to each external particle irrespective of its own proper 
charge. In this way we can include these anomalous transitions in 170 and 19 F 
within a single scheme (Elliott and Flowers 1955). An even more picturesque 
alternative view is to think of a neutron’s motion as that of a neutral volume 
in a nucleon fluid of mean charge density e/2 per nucleon volume, which it 
displaces as it moves so that its motion is equivalent to that of a body of 
charge —e/2 moving in a vacuum. Some such picture is held by many people 
who contemplate, say, the independence of E3 transition speed in heavy nuclei 
on whether the jump is that of a neutron or a proton. The present writer first 
heard of such a picture from O. R. Frisch. It is a very attractive picture but 
one which seems difficult to fit in with the idea of shell model orbits since it 
seems to depend on one nucleon rather actively pushing others out of its way. 
Perhaps the attempts of Brueckner and his co-workers to reconcile the shell 


model with strong nucleon-nucleon interactions should permit a quantitative 
formulation of this notion, 
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than the E2 which, as we have seen, is quite probably at 14:7 Mev. It is. 
still possible that the transition from this 3— state at 6-14 Mev is speeded 
by some admixture of the remote one-phonon state (just as that from the 
2-- state at 6-91 Mev may benefit by admixture of the one-phonon state 
at 14-7 Mev). However another possibility resides in the possible applica- 
bility of the alpha-particle model for 160 (Dennison 1954). According to 
this model the 3— state is one of simple rotation of the basic tetrahedron, 
and although the dynamical properties of the model have not yet been 
fully worked out, it seems rather likely that it will predict a rather fast 
H3 decay for this state. This transition will in fact provide a good test 
for a model that in other respects is not without its difficulties. 


§ 5. CONCLUSIONS 


We may summarize the findings of this survey of the radiative properties. 
of nuclei of A <20: 


(i) The [7 |? values of E1 transitions show a (logarithmic) distribution 
which is rather symmetrical about the value 0-032 and of a width such 
that 85°, of observed transitions have a strength within a factor of seven 
of this value. Isotopic spin forbidden transitions are excepted from this 
classification and tend to show much smaller widths. 


(ii) The |M|? distribution for M1 transitions centres on [M |?~0-15 
from which an excursion by a factor of 20 in strength in either direction 
is needed to cover 85°% of the observed transitions. 


(iii) These E1 and M1 distributions are consistent with expectation 
based on crude shell-model computations. 


(iv) It is in general impossible to distinguish between El and M1 
transitions on the basis of radiative width though if a transition shows 
a value of I’, (ev)/Z,? (Mev) of greater than 0-02 then it is very probably 
an E1 transition. 

(v) There is no clear dependence, in the 1p shell, of dipole radiative 
width on level density. There is some suggestion from this and from the 
shape of the ||? distribution for E1 transitions that most states so far 
met with in light nuclei are rather simple ones. In this connection it is 
urged that weak transitions are of considerable interest and merit study. 


(vi) From the few E2 and E3 transitions so far measured there seem to 
be clear signs that collective mechanisms such as surface oscillations may 
have a place in the structure of these light nuclei; the E2 or K3 transition 
rates predicted by the shell model fall well short of experiment in several 


cases. 
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ABSTRACT 
A method of measuring the dynamic tensile strengths: of liquids, 
involving the use of an electrical pressure bar, is described. Cavitation 
phenomena in water subjected to large transient tensile stresses have been 
observed using a spark shadowgraph technique. Under similar conditions 
of transient loading the cavitation thresholds of water and of glycerol are 
approximately 17 x 10° dynes/em? and 63 x 10° dynes/cm? respectively. 


§ 1. INTRODUCTION 
EXPERIMENTAL determinations of the tensile strengths of liquids have 
been made by a large number of workers, using a wide variety of tech- 
niques, during the last hundred years. Critical discussions of some of 
the methods have been given by Temperley and Chambers (1946) and by 
Blake (1949). These surveys reveal that the tensile strength of a liquid 
may vary by a factor of 100 or more, depending both on the observer and 
the method; for example, the maximum value claimed for the tensile 
strength of water is 200 atmospheres (Briggs 1950) while the lowest is 
1 atmosphere. It seems reasonable to conclude that the experimental 
knowledge of the tensile strenths of liquids is meagre and uncertain. 
Even the highest experimental value claimed for water falls far short 
of the theoretical predictions of Temperley (1947), Fisher (1948) and 
Fuerth (1941). Although these theoretical estimates are deduced using 
very different approaches to the problem, all three yield values of the 
tensile strength of water at 17°c of approximately 1000 atmospheres. 
The principal difficulty in realizing high tensions in liquids arises from 
the presence of gaseous or solid nuclei in the liquid. Just as the tensile 
strength of brittle solids is determined by the size and distribution of 
Griffith flaws in the material, so the maximum tension a liquid can with- 
stand is limited by the size of gaseous bubbles or unwetted nuclei of 
solid matter present in that liquid. In practice, the tensile strength of 
a liquid under static loading is governed by the size of the largest bubble 
or nucleus present. 
The majority of methods for the determination of the tensile strength, 
or cavitation threshold, of a liquid employ static or quasi-static loading, 
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and bubbles or voids within the sample have time to grow as the tension 
increases. If, however, a transient tensile stress of short duration is 
applied to a liquid, any bubbles which tend to become unstable may not 
have time to enlarge sufficiently, either to be observed photographically 
or to affect the further transmission of tension. A liquid may therefore 
remain non-cavitated for quite large stress amplitudes. This effect has. 
been described by Kolsky, Lewis, Sampson, Shearman and Snow (1949) 
in work on the surface phenomena attending underwater explosions. 
The present paper describes two types of experiment which are part of a 
programme of research to investigate cavitation effects in liquids caused 
by large transient tensile stresses. The experiments reported were made 
to determine the cavitation thresholds of water and of glycerol under 
dynamic loading conditions. 


§ 2. EXPERIMENTAL METHODS 


The technique for producing a large transient tensile stress in a liquid 
is based on the fact that a compression wave incident normally at a free 
surface is reflected as a wave of tension, which travels through the tail of 
the compression pulse. The resultant stress at a point within the liquid, 
at any instant during reflection, may be found by adding the amplitudes 
of the incident and reflected waves at that point. Figure 1 illustrates 
this type of reflection for a steep-fronted wave, represented at any point 
by a discontinuous rise in pressure followed by an exponential decay in 


Fig. 1 
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Reflection of a shock wave at free boundary, 
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time. ABC represents the pressure characteristics of such a wave, which 
has been partially reflected at the free liquid surface DOH , the charac- 
teristic being shown as a function of distance. The leading portion ABH 
of the shock wave is reflected into the liquid as a tension wave DEF. and 
in this particular case the stress in the liquid at any depth between O and 
F is tensile, and equal in magnitude to (DE—-HG) at the depth concerned. 
When the value of (DE-HG) exceeds the cavitation threshold of the 
liquid it ruptures and ideally, if the amplitude of the incident wave is 
sufficiently great, a layer of liquid is thrown up with part of the momentum 
of the incident wave trapped in it. Reflection of the remainder of the 
wave then occurs at the new surface created by the detached layer. 

By observing the depth at which cavitation occurs, using high-speed 
photographic methods, during the reflection of a wave of known pressure- 
time characteristics or by measuring the pressure changes during reflection, 
it is possible to deduce the cavitation threshold or dynamic tensile strength 
of the liquid. These two methods are described in this paper. In both 
cases the test liquid was contained in a vessel through the bottom of which 
protruded a steel or Perspex rod, supported vertically. The firing of an 
electric detonator in close contact with the lower end of the rod resulted 
in a large amplitude compression pulse travelling up the rod and into 
the liquid. 

The electrical pressure bar technique for measuring transient pressure 
waves in liquids, described in section §3, was developed by Professor 
R. M. Davies and his colleagues at University College of Wales, 
Aberystwyth. The method has been used both for measuring transient 
pressures in detonating gases (Davies, Edwards and Thomas 1950) and 
for the measurement of the tensile strengths of liquids under dynamic 
stressing (Davies, Lewis, Rees and Trevena 1955). 


§ 3. THe MEASUREMENT OF PRESSURE CHANGES BELOW A FREE LIQUID 
SuRFACE DuRING THE REFLECTION OF A COMPRESSION PULSE 


3.1. Apparatus 

(a) General 

The arrangement of the apparatus is shown diagrammatically in 
figure 2. The test liquid was contained in a steel cylinder, length 18 in. 
with a nominal bore of 1 in. and outer diameter 1} in. A silver-steel rod 
(hereafter called the firing rod), 1 in. diameter and length 4 ft. was sealed 
into the bottom of the cylinder. Both ends of the rod were carefully 
ground. The bore of the cylinder was such that the rod made a good 
sliding fit into it. A stress wave was produced. in the firing rod by 
exploding a Nobel’s No. 6 electric detonator in close contact with a 
cylindrical hardened steel anvil wrung on to its lower face. 

Transient pressures at a point in the liquid were measured using an 
electrical pressure bar with a parallel-plate condenser microphone (Davies 
1948). The pressure bar, of silver steel, diameter Lin. and 3 ft. long, 
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was sealed into the side of the cylinder at its mid-point. The face of the 
measuring bar in contact with the liquid was machined flush with the 
inner wall of the cylinder. It was essential that neither the initial 
compression pulse set up in the firing rod nor the pulses produced in the 
measuring bar should be distorted during passage through the container 
seals. Following numerous tests, the most suitable seals for this purpose 
were found to be Weston Oil seals, which are manufactured complete with 
metal housing. The characteristics of the stress wave in the firing rod 
could be measured using a cylindrical condenser microphone placed at any 
desired position on the bar. 


Fig. 2 
_TEST LIQUID 
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C/M CRYSTAL OSCILLATOR 


PRESSURE BAR 


U FEED UNIT 
FIRING ROD 
PHOTOCELL DELAY SWEEP UNIT 


General arrangement of apparatus. 


MPLIFIER 


The output of the electrical pressure bar feed unit was amplified and 
displayed upon one set of Y-plates of an independent-gun, high-speed 
cathode-ray oscillograph. A 100 ke/s crystal oscillator, connected to the 
second pair of Y-deflector plates, provided a time calibration for the 
traces. Both pairs of X-plates of the cathode-ray oscillograph were 
connected across a balanced, single-sweep time base unit (Kolsky 1949). 
A brightening pulse for the oscillograph was derived from the Sweep unit. 
The oscillograph traces were recorded with an f/1-0 lens camera upon 
Ilford H.P.3 plates. 

Light from the detonator explosion was picked up by a photocell and 
the resulting pulse amplified before being delayed by a variable thyratron- 
controlled RC delay circuit. A sharp-fronted pulse from the cathode 
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circuit of the final thyratron triggered the Sweep unit a few microseconds 


before the stress pulse in the measuring bar arrived at the condenser 
microphone. 


(6) The electrical pressure bar 


Descriptions of the use of an electrical pressure bar for the measure- 
ment of transient pressures have been given in papers by Davies (1948) 
and Kolsky (1949). It is shown that if a plane compression pulse is 
propagated in a cylindrical bar under the conditions that (a) the pressure 
does not exceed the elastic limit of the material, and (b) the pulse is long 
in duration compared with the time taken to travel a distance equal to 
the diameter of the bar, then the following relations hold 


c=(H/p)¥? and p=pcv 


where c=velocity of propagation of the pulse, #=Young’s modulus for 
the bar material, v—particle velocity at any point, p—pressure amplitude 
at that point, p=density of the material. 

At the free or measuring end of the bar a compression pulse is reflected 
_ as one of tension and vice versa. The velocity of the end of the bar is 
therefore twice the corresponding particle velocity associated with a 
pulse travelling in the bar remote from either extremity. Thus if the 
displacement of the end of the bar is x at time ¢ after the arrival of the 
pulse we have 


A parallel-plate condenser microphone records the variation of x. If the 
time interval during which the variation takes place is known, together 
with pc for the bar material, then the pressure—time relation of the pulse 
may be found. 

In the experiments described here the pulse duration was about 
35 microseconds, and the maximum pressure recorded with the bar 
approximately 7-7 107 dyn/cm? or 0-5 tons/in.?. The time taken for 
a compression pulse to travel 4 in., the diameter of the pressure bar, 
is about 2 microseconds, and the elastic limit of silver steel is of the order of 
30 tons/in.2; thus the pulses were propagated through the pressure bar 
with very little distortion. 

Ideally, the sensitive element of a mechanical device for measuring 
transient pressures should be dimensionless, in order that the pressure 
over the element should be uniform at all times during the passage of the 
transient across it. A 35 microsecond pulse in water has a length of 
approximately 5 cm, whilst the diameter of the bar used in the experi- 
ments was 1:27 cm, hence the ideal conditions were not realized here. 
This factor was, to a certain extent, offset by the use of a cylindrical 
condenser microphone on the firing rod, which provided a check of the 


characteristics of the pulse entering the liquid. 
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3.2. Results 

Two series of experiments were performed using first untreated tap 
water and then glycerol as the test liquids. The procedure was first to 
line up the sensitive face of the measuring bar flush with the inner face of 
the cylinder. The cylinder was then filled with the test liquid to a known 
height above the end of the measuring bar and allowed to reach tem- 
perature equilibrium with the surroundings. The separation between the 
insulated plate of the condenser microphone was adjusted with slip 
gauges. The detonator was then fired and the output of the condenser 
microphone feed unit, suitably amplified, recorded with the high-speed 
cathode-ray oscillogrpah. This was repeated for various distances of the 
measuring bar below the free surface. In order to determine the initial 
form of the compression pulse in the liquid, the cylinder was filled until 
the depth of liquid above the measuring bar was 20 cm, and the procedure 
repeated. Under these conditions there was no interference, at the 
pressure device, between the initial wave and any reflected wave. 

Figure 3(a)(Plate 9)is an enlargement of a plate taken with 20 cm of water 
above the end of the measuring bar. The upper trace is a timing wave 
with a frequency of 100 ke/s, increasing time being measured from left to 
right. The lower curve is due to the amplified transient voltage from the 
condenser microphone. It will be seen that after the initial rise in voltage 
due to the compression pulse, the voltage remains constant and is followed 
by a number of oscillations due to dispersion in the bar. In fig. 3(6). 
which is an enlargement of the record obtained with 1-0 cm of water 
above the measuring bar, the initial rise in voltage, due to the initial 
compression pulse, is followed by a sharp decrease, lasting for about 
10 microseconds, and then a more gradual and sustained decrease with 
small oscillations superimposed on it. The sharper decrease is due to the 
effect of reflection of the compression pulse at the free surface. To deter- 
mine the variation of pressure below the surface with time at a given 
point, the displacement time curves were differentiated numerically. 

Figure 4 gives curves showing the variation of maximum observed 
tension in the two test liquids with distance of the measuring bar below 
the free surface. No existence of tension was observed at distances 
greater than about 4:3 cm below the free surface. The maximum tensions 
observed were 18x 10%+5%) dyn/cm? for water and 63 x 10%(+5%) 
dyn/em? for glycerol. It was considered that possibly these figures gave 
only an indication of the tension at which these liquids cavitated at the 
face of the measuring bar. Therefore, the high-speed photographic 
method of determining cavitation thresholds, described in $4, was 
used to investigate cavitation in water. 


§ 4. PHOTOGRAPHY OF CAVITATION IN WATER 
The spark shadowgraph method was used to record the type of transient 
cavitation with which this paper is concerned. This method is useful 
because not only may very short exposure times (<1 microsecond) be 
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obtained using spark gaps in air but the auxiliary equipment is exceedingly 
simple. . Further, the density changes in a fluid caused by a pressure 
wave give rise to corresponding changes in refractive index and, owing 
to the deviation of light caused by these changes, the pressure distribution 
in the fluid may be ascertained at any instant. 


Fig. 4 
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4.1. Apparatus 

The liquid was contained in a large, rectangular metal tank with two 
parallel plate-glass sides. A Perspex rod, diameter 1 in. and 4 ft. long, 
supported vertically, passed through a Weston Oil seal into the centre of 
the base of the tank. An electric detonator, fired against a } in. hardened 
steel anvil wrung to the lower face of the Perspex rod (firing rod), resulted 
in a compression wave travelling up the rod and into the liquid. Perspex 
was chosen as the material for the firing rod since its acoustic impedance 
is only about twice that of water, giving a pressure amplitude transmission 
of approximately 70%. 

Supported parallel to and 8 cm away from one of the glass windows of 
the tank was a photographic half-plate behind a mechanically operated 
shutter. The spark shadowgraph was recorded on this plate by light 
from an enclosed point-source spark gap, positioned 98cm from the 
other glass window and level with the surface of the liquid. The spark 
gap electrodes, with an emission area of 4x 10-2 cm?, were constructed 
of aluminium and supported in a cylinder of similar material. The spark 
was produced by the discharge of a 0-1 microfarad condenser, charged to 
10 ky, across the electrodes. The condenser was specially selected for 
its low resistance supports and plates. The mechanical arrangement of the 
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spark gap and condenser was such that the inductance of the discharge 
loop was a minimum. 

Light from the detonator explosion was picked up by a photocell and 
the ensuing electrical impulse used to trigger both a microsecond counter 
and a preset RC delay circuit. The delayed pulse triggered the spark gap. 
A second photocell unit was triggered by part of the light emitted from 
the spark gap and provided a pulse which stopped the microsecond counter. 
A record was therefore obtained of events occurring in the tank and of the 
time between the firing of the detonator and the peak output of the 
spark gap. 

The characteristics of the initial pulse of compression in the firing rod 
were measured with a cylindrical condenser microphone. The surface of 
the Perspex rod was used as the earthed conductor of the condenser 
system and for this purpose was coated with a smooth, even layer of 
graphite over required regions. The determination of the shape and 
amplitude of the compression pulse in the liquid, as a function of distance 
from the upper face of the firing rod, was made in auxiliary experiments. 
An electrical pressure bar was sealed horizontally into one side of the tank 
in such a way that only its sensitive face, machined plane, was in contact 
with the liquid. The face of the bar was positioned perpendicular to, 
and directly above the upper face of the firing rod. The distance between 
the measuring bar and firing rod was varied by moving the latter 
vertically. 

4.2. Results 

With a depth of 8 cm of water above the end of the firing rod a series of 
shadowgraphs was obtained of cavitation effects at increasing times after 
the firing of the detonator. The first five photographs of fig. 5 (Plate 10) are 
of such shadowgraphs, with times after detonation given below each. 

The plate at 537 microseconds was taken before the compression wave 
had been transmitted to the water. The upper end of the firing rod is 
shown at the bottom of the picture, and the free water surface is indicated 
by the horizontal black line near the top. At 575 microseconds the 
compression wave has been transmitted to the water. In this photo- 
graph appear two sets of bands, each set consisting of a dark band followed 
by a lighter one. The upper set of bands, which are rather faint, but 
which can be clearly seen on the original plate, is very narrow and 
represents the position of the head of the compression wave. The origin 
of the broader bands will be discussed later. This identification of the 
upper bands is supported by the evidence in the sixth photograph. Here 
a large depth of water was used and a steel plate arranged obliquely in the 
tank. The photograph shows the narrow set of bands, representing the 
head of the wave, near the upper edge of the photograph and also the 
reflection of the wave at the steel plate, giving rise to a wavefront moving 
from left to right in the picture. Referring again to the photograph 
taken at 575 microseconds, a disc of cavitated water appears on the face of 
the firing rod with a few cavitation bubbles above it. 
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At 592 microseconds the head of the compression wave has advanced 
further into the water, and a multitude of secondary wavelets appear to 
originate from a point a few centimetres below the end of the rod. The 
photograph at 611 microseconds shows the cavitation about to 
develop below the free surface, since part of the initial wave has already 
been reflected. At 625 microseconds the cavitation has developed quite 
considerably. The bubbles at the lower edge of the cavitated region are 
smaller than those in the centre, the latter region having been under 
hydrostatic tension for a greater period. A region of non-cavitated 
liquid exists immediately below the free surface, indicating that here the 
liquid was able to withstand the transient tensile stress applied to it. 

The broad set of bands appearing in some of the photographs are 
assumed to indicate the tail of the compression pulse, although the reason 
that they are so well defined is uncertain. The large bubble or group 
of bubbles which always forms on the face of the firing rod grows with 
time and then collapses, the whole process lasting about 1-5 milliseconds. 
There are two possible explanations of this effect. First, it is known 
(Davies 1948) that the very sharp compression pulse produced by an 
explosive in contact with the end of a cylindrical bar is lengthened as it 
travels along the bar. For steel the pulse reaches an approximately 
stable form after it has travelled a distance equal to about twelve times 
the diameter of the bar; for Perspex this distance is much greater owing 
to the added effects of the absorption. The stable form of the pulse 
consists of a main compression pulse followed immediately by a train 
of high frequency oscillations which extend into the tensile region. The 
reflection of the tensile portion of the pulse at the Perspex-water inter- 
face results in the end of the rod moving momentarily away from the 
liquid, this movement being sufficient to cause cavitation at the end of the 
rod. The second possibility arises from the fact that here we are dealing 
with hydrostatic pressures, which act equally in all directions. The 
liquid, while under compression, tends to flow radially outward from the 
face of the rod, leaving a vapour-filled cavity at the interface. 

The spark shadowgraphs show only the central section of the cavitated 
zone, since the width of the photographic plates was only one-third that 
ofthetank. The full extent of the cavitated zone is shown in fig. 6, Plate 11. 
The photograph shows the transient cavitation in water a few micro- 
seconds after the initial compression wave has been completely reflected 

-at the free surface immediately above the rod. The event was recorded 
with a quarter-plate camera, the illumination being provided by a 
Mullard L.S.D.2 flash tube. These flash tubes give an effective exposure 
of about 1 microsecond, if adequate care is taken in the selection of 
discharge condenser, and the inductance and resistance of the discharge 

loop kept to the minimum. The main cavitated zone is lenticular in 
form. A region of non-cavitated water remains between this zone and. 
the water surface, indicating that here the water has been able to with- 
stand the applied tension. The height of the upper limit of the cavitated 
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zone decreases with distance out from the axis of the firmg rod. This 
effect might be expected from the geometry of the system. The greater 
the distance from the axis of the bar the weaker the reflected wave 
becomes, and has therefore to travel a greater distance through the tail 
of the compression pulse before the necessary negative pressure to cause 
visible cavitation is reached. The lower limit of the cavitated zone 
indicates that the front of the compression wave is approximately spherical 
in form. 

In order to estimate the cavitation threshold of water from the spark 
shadowgraphs, it was required to know the shape and amplitude of the 
wave impinging on the free surface. Using the arrangement described 
earlier, the transient pressure was measured as a function of the vertical 
distance R above the firing rod. The depth of water above the measuring 
bar and the dimensions of the tank were such that the initial wave was 
measured before any reflected waves arrived at the detector. The 
pressure amplitude at zero distance was computed from a knowledge of 
the stress amplitude in the Perspex rod, given by measurements with the 
cylindrical condenser microphone, and the relative acoustic impedance 
of Perspex and water. The pressure amplitude varies as k~°® for values 
of R between 0 and 2-5cm and as R-}? for R between 2-5 and 8 cm. 
Theoretically the propagation of a small-amplitude spherical wave from 
a point source in an ideal infinite fluid results in the pressure amplitude 
decreasing as R-! where # is the radial distance from the source. In_ 
the apparatus described here, the arrangement was equivalent to a piston- 
like source in an infinite medium for small values of R. For distances 
from the piston comparable with the piston diameter, the propagated 
wave is intermediate between a plane and a spherical form, hence the 
pressure amplitude does not decay as rapidly as R-!. At larger values 
of Rk, departure from the acoustic case becomes more pronounced for a 
different reason; the pressure amplitude falls off more rapidly than R-?, 
owing to absorption brought about by inhomogeneities in the water. 

The simplest and most direct way of estimating the cavitation threshold 
P, from the shadowgraphs is to base the measurements on the upper 
limits of the cavitated region. Figure 8 shows diagrammatically the 
stress distribution in a liquid with an infinite value of P, when a com- 
pression wave has been completely reflected normally at a free surface. 
The stress distribution is indicated by the solid lines and is shown for two 
waves of differing amplitude. If now the liquid cavitates under a 
tension P,, the upper limit for cavitation is the line AB in the case of 
wave (a) and CD for wave (b). These limits are indicated in the shadow- 
graphs. 'To convert measurements on the shadowgraphs to corresponding 
distances in the central vertical plane of the tank a spark shadowgraph 
was taken of a wire grid placed in this plane. The shadowgraph had a 
magnification of 1-2 diameters. 

Measurements made with water depths of 8 and 10 cm indicated that 
the minimum tensile stress required to cavitate water was about 
15 x 10®dyn/cm?. 
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§ 5. DISCUSSION 

It is of interest to compare the cavitation threshold of water, found 
from photographic measurements, with the maximum tensile stress of 
18 x 106 dyn/cm? recorded using the pressure bar technique. The photo- 
graphic technique gives a measure of the minimum tension necessary, 
at the particular loading employed, for nuclei in the water to become 
unstable and grow to a visible size. The electrical pressure bar measure- 
ments indicate the extent to which a cavitated liquid may continue to 
transmit tension. From these considerations it is to be expected, as 
found in the experiments described here, that the electrical pressure bar 
would yield slightly higher results. 

The equation of static equilibrium for a spherical bubble in the body 
of a liquid is i aT 
where p, and p, are the external and internal pressures respectively, and 
cis the macroscopic value of the surface tension. The condition that the 
bubble should grow is that py, be reduced below the value p;—2o0/R. 
Taking the mean value of 16x 10° dyn/em? as the cavitation threshold 
of water, the minimum radii of the bubbles which become unstable under 
tensions greater than this were of the order of 10°? cm. This is many 
times the molecular dimensions, hence the cavities which became unstable 
were almost certainly gas or vapour filled bubbles. 

The lowest resonant frequency f of a bubble in a liquid is given by 


Smith (1935) as 1. [8p (20/r)]—20/r) 12 
i =) p } 

where p is the density of the liquid, p, the hydrostatic pressure, y the 
ratio of specific heats of the gas, and o the surface tension. The lowest 
natural pulsation period of an air or vapour filled bubble in water with a 
radius of 10°-°cm is thus about 10-8 seconds. The duration of the 
transient tensile stress was about 30 microseconds, i.e. much greater than 
the pulsation period of the bubbles or nuclei, which therefore had time 
to enlarge. Strasberg (1955) has shown that for water the maximum 
radius Rmax to which a bubble grows is given by the approximate 


relation , 
Rmax=7 { Ap/p ye 


where Ap is approximately numerically equal to the peak amplitude of 
the transient pressure and 7’ its duration. In the experiments described 
in §4 the peak amplitude of the negative pressure transient was about 
3-3 x 10’dyn/cm?. Substitution in theequation above yields Rnax=0-17em. 
This value is in good agreement with the maximum radius of 0-1 em 
for bubbles observed in these experiments. 

The maximum tensile -stress recorded in glycerol was approximately 
6-3 10’ dyn/cm?, i.e. roughly four times that recorded in water. At 
17°c the surface tensions of glycerol and water are about the same, whilst 
the viscosity of glycerol is nearly fifteen hundred times that of water. 
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As no attempt was made to degas or purify these liquids it would seem 
reasonable to assume that viscous forces markedly affect the growth of 
cavities or bubbles in liquids under transient tensile loading. 
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ABSTRACT 


Three semi-circular focusing spectrometers, making use of a common 
magnetic field, have been used to measure the angular correlation of the 
internal conversion electron pairs produced in the decay of the 1:85 Mev 
state of 88Sr. The results support the assignment 2+ for this state. 
Similar equipment has been used to measure the angular correlation of 
certain conversion electron cascades. With a polarimeter having good 
energy resolution, direction—polarization measurements have been carried 
out on the 0-90, 1:85 Mev gamma-ray cascade in **Sr. These show that 
the 1-85 and 2-75 Mev levels in this nucleus have even and odd parity 
respectively. 


$1. LyTRODUCTION 


THE measurement of gamma-ray angular correlations has proved a 
powerful technique for determining nuclear spins and multipole mixture 
ratios (Kraushaar and Goldhaber 1953, Steffen 1954), particularly in 
even-even nuclei, for which the ground state spins are known to be zero, 
and transitions to the ground state must be pure multipoles. Occasionally, 
however, cascade transitions occur in which the conversion coefficients. 
are so great that only conversion-gamma or conversion—conversion 
correlation measurements will give adequate coincidence rates. Such 
measurements encounter several difficulties. If the conversion coefficient 
of the final transition is appreciable (25°) the lifetime of the inter- 
mediate state is likely to be long enough (210-° seconds) for perturbation 
effects to be serious. To minimize scattering the sources must be very 
thin, so that material of high specific activity is required, and the techniques 
successfully employed to minimize perturbations in gamma-gamma 
correlations are difficult to apply. Finally, conversion electron spectra 
are usually rather complex, so that apparatus with good energy resolution 
is required to isolate and observe particular transitions. Siegbahn (1952). 
has described a system of two thin lens ring-focusing spectrometers which 
es be at for this purpose. Equipment described later in this paper 
achieves the same end by much simpler i 

eek Diese ee pler means, though with some loss in 


ee 
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The angular correlation between internal conversion electron pairs 
(Horton 1948) can also be measured with this equipment. Since such 
a measurement involves only one transition its interpretation is fairly 
straightforward, and perturbation effects do not arise. On the other 
hand the energy of the transition must be well in excess of 1 Mev, and the 
conversion coefficients are very small (510-3). To avoid interference 
from external pair creation very thin sources are again required. The 
most serious obstacle to such measurements is the fact that the transition 
in which the pairs are produced is usually preceded by beta radiation with 
a maximum energy exceeding that of the pair spectrum. Calculations 
suggest that even under these conditions it should be possible, by using 
the shortest resolving times at present available, to obtain pair coincidence 
rates approaching ten counts per hour which account for the greater part 
of the total coincidence rate. Because of the strong dependence of this 
correlation on multipole order and parity, such rates will still be large 
enough to determine these parameters. 

For the determination of relative parities the most generally useful 
correlation is the gamma _ polarization—direction correlation. The 
realization that the standard polarimeter can be used as a low resolution 
Compton spectrometer (Bishop and Perez 1955, Wood and Jastram 1955) 
has enabled this type of measurement to be carried out on moderately 
complex decay schemes. In the present paper a technique is described 
which gives considerably better resolution without any loss in the overall 
efficiency of the polarimeter. 


§ 2. EQUIPMENT FOR ELECTRON ANGULAR CORRELATION 
MEASUREMENTS 


The apparatus to be described makes use of conventional semi-circular 
focusing in a uniform magnetic field. Figure 1 shows electron orbits, 
projected on a plane perpendicular to the applied field, for two spectrometer 
systems, with focusing radii 7, and ry, receiving electrons from a common 
source. It is clear that if one of the spectrometer systems is rotated 
about the field line through the source the angle @ between the average 
directions of emission of the electrons accepted by the two systems will 
be changed, but the energy selection conditions will be completely 
unaffected. The magnetic field can be adjusted to bring electrons of 
any desired energy into focus in one of the spectrometers, and electrons 
of some other energy, which can also be given any desired value by 
suitable choice of the ratio 7/75, will be simultaneously focused in the 
other. Thus this system provides all the essential facilities of what 
Siegbahn has called a ‘ spectro-goniometer ’. It does, however, expose 
the source to a field of up to a few hundred oersteds perpendicular to the 
plane of observation, which could produce an appreciable perturbation of 
the angular correlation to be measured if the intermediate state half life 


were > 10-8 seconds. 
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Geometrical conditions in the spectrometers can be made broadly similar 
to the ‘high source ’ specification (Geoffrion 1949), with one important 
exception. To collect electrons with equal efficiency from all points on 
the surface of the source, each exit slit must be wide enough to accept 
the whole image of the source, including all its aberrations. To reduce 
these it is desirable to reduce the height of the source to about half the 
height of the exit slit for the spectrometer of smaller radius. Failure to 
satisfy this condition leads to the coincidence line-sharpening effect noted 
by Fowler and Shreffler (1950), and might easily produce spurious 
anisotropies. It is experimentally convenient to use a fixed value for the 
greater of 7;, y, so that each measurement requires a new set of baffles to 
give the appropriate value of 7/75. 


Fig. 1 


Layout for electron—electron angular correlation measurements. 


For the experiments described below the solid angle of collection was 
fixed at 0-12 steradians, and the greater of the two spectrometer radii at 
6 cm. This allowed the whole system to be set up in an evacuated 
cylindrical chamber 11 in. in diameter and 2} in. high. Each spectrometer 
unit was made up of a set of baffles, an exit slit, a geiger counter, and a 
gamma-ray shield, and could be rotated by making the appropriate 
changes in the positions of these four components. Since there was very 
little room in the experimental chamber, small geiger counters were made 
up with side windows of 2} mg/cm? mica covering bridged slots, 2 in. 
high and 7 mm side, just large enough to admit the image of the source 
produced by the 6 cm radius spectrometer. These counters were found 
to be capable of working into coincidence units of resolving time 
27=4 10-7 seconds without appreciable loss of coincidences even in fields 
as great as | kilo-oersted. 
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$3. CASCADE CONVERSION ELECTRON CORRELATION MEASUREMENTS 


The following results are given primarily as examples of the use of 
the above equipment. All three involve L-shell conversion electrons, for 
which theoretical calculations are not yet available. The tellurium 
results were obtained with a source thick enough to produce appreciable 
distortion of the conversion lines, so that the original anisotropy is 
appreciably reduced by scattering. The intermediate state lifetime for 
‘Hg is long enough, as the measurements show, for perturbation effects 
to be serious. 

21Te.—The angular correlation of the 82:1L; 213K cascade in this 
nucleus, measured with a tellurium chloride source, was found to be 


W(180°) : W(150°) : W(120°)=(1-00-L0-05) : (0-91-05) : (1-01-05). 


™Te—The angular correlation of the 88-5L ; 159K cascade in this 
nucleus, obtained from results similar to those for 121Te, and corrected for 
finite solid angles but not for scattering, was found to be 


W(0)=1-+ (0-16-+0-07) cos? 6. 


1% Ho.—The angular correlation of the 164L; 133L cascade in this 
nucleus was found to be given by 


W (180°) : W(150°) : W(120°)=(1-00-+0-02) ; (1-02-L0-02) : (0-98-+0-02). 


The source was a mixture of HgCl, and Hg§, thin enough to give sharp 
lines in the spectrometers (fig. 2). The small observed anisotropy con- 
trasts sharply with the 20% effect found by Frauenfelder et al. (1950), 
who used absorption techniques for energy selection, and prepared their 
sources by simultaneous condensation of gold and mercury vapour. 
Though the intermediate state is metastable (t,j.—7 x 10~° sec) the direct 
coupling between the spectrometer field (~200 oersteds) and the nuclear 
moment is much too small to account for the discrepancy. Indirect 
coupling by way of a non-zero electronic moment may however cooperate 
with the effects of atomic and crystalline fields in destroying the 
correlation. 


§ 4. CONVERSION ELECTRON ParR CORRELATION IN **Sr 


For this measurement three spectrometer systems, all of the same radius, 
were arranged, one to select positrons, the second to select negatrons 
emitted in the same average direction as the selected positrons, and the 
third to select negatrons emitted at an average angle of 110° to the 
selected positrons. The spectrometer field was then scanned through the 
neighbourhood of the pair line. With this equipment the ratio 
W(0°)/W(110°) (see fig. 3) was measured for the 1-85 Mev transition in 
88S; To correct for the effects of external pair creation two sources 
were used, of which one had an average surface density twice that of the 
other. The experimental value for W(0°)/W(110°), corrected for external 
pair creation, was found to be 4:5-+1:0. Theoretical values for this ratio, 
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in the case where the positron and negatron have equal energies, were 
taken from a paper by Horton (1948) using the Born approximation 
(Horton and Phibbs 1954). When modified for the experimental angular 
resolution they become 2-2-++0:3, 7-5-+0°9, and 5:5--0-7 for El, M1, and 
E2 multipoles respectively. Thus the measurement supports the assign- 
ment 2+ for the 1-85 mev level in **Sr. 


Fig. 2 
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Conversion electron lines of 17Hg observed in electron—electron angular 
correlation equipment. 


$5. EQUIPMENT FOR POLARIZATION—DIRECTION CORRELATION 
M®ASUREMENTS 


The polarimeter is physically similar to the system described by 
Metzger and Deutsch (1950). The scatterer is a truncated cone of 
tetraphenylbutadiene in styrene 2 in. long, with a base diameter of 1? in. 
and a semi-vertical angle of 10°. The side counter is a 2 in. x 2 in. x 1} in. 
block of thallium activated sodium iodide, with its centre 2+ in. from the 
polarimeter axis. The polarimeter remains fixed, and the polarimeter 
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angle is changed by rotating the other directional counter about the 
polarimeter axis, SO that the system contains only one moveable counter. 
Figure 4 shows the experimental arrangement. The polarimeter axis is 


Fig. 3 
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MAGNET CURRENT ————> 
The **Sr internal conversion pair line observed at 6=0° and 6=110°. 


vertical. The bearings of the swinging arm carrying the unpolarized 
detector are fastened to a horizontal ring mounting, which rotates to give 
the polarimeter motion. The directional angle @ can be set at any multiple 
of 15° from 90° to 180°, and the apparatus then records data automatically 
at polarimeter angles following the sequence... —90°, 0°, 90°, 180°, 
—90° ..., and inserts a measurement of the accidental coincidence rate 
after every third run. 

The energy discrimination of the polarimeter is obtained by noticing 
that, since the energies of the scattered quanta reaching the side counter 
are rather small, many of them release all their energy in the NaI—Tl 
phosphor. Thus if pulses from the scattering counter and the side 
counter are treated so that a given energy release in either phosphor 
produces a pulse with the same height, shape, and delay, and are then 
added together, the spectrum of the summed pulses obtained in this way 
from a source of mono-energetic gamma-rays will show a peak at the high 
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energy end (see fig. 5). This peak arises from those events in which the 
energy of the scattered gamma-ray is completely absorbed in the side 
counter, and shows approximately the same resolution which the scattering 
phosphor would display in a two crystal spectrometer with a well defined 
scattering angle. Those scattered quanta which lose only part of their 
energy in the side counter produce a continuum which extends down from 


Fig. 4 
UNPOLARISED 
DETECTOR 
@ VARIATION fk 
POLARIMETER | as 
AXIS De ex" ) 
ay WA 7 
Ae apa h eS g VARIATION 
oS se po oy ~ oe Be 
Pin 2 A ee Tat 
a -_ -—- - wow K—-e — — -)-\ 
ae ire <s | a \ 
~ we %, fate - 
/ oan Gada te I ; 
/ i 1 
| Es 
/ - & 
' SOURCE ; yn 
pee... ae a 


NORMAL PLANE 
TO POLARIMETER AXIS 


POLARISED DETECTOR 
SCATTERING PHOSPHOR 
SIDE COUNTER 


Diagram illustrating the mechanical arrangement of the components for 
polarization—direction correlation measurements. 


the total absorption peak to a cut-off at the low energy end. This cut-off 
corresponds to the energy of the Compton electrons produced in the 
scattering phosphor by those quanta which are scattered back only just 
far enough to reach the side counter. To eliminate the very large back- 
ground produced in the summing unit by non-coincident pulses from one 
or other of the polarimeter detectors, only those pulses are selected which 
are accompanied by coincidences between these detectors. The low- 
energy cut-off visible in the polarimeter spectrum shown in fig. 5 is partly 
due to amplitude dependent time shifts in the auxiliary coincidence unit 
(2r7~2 10°" sec) used for this selection. This summing technique has 
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the advantage over the technique of gating the pulse distribution in the 
scattering counter (Bishop and Perez 1955, Wood and Jastram 1955) 
that the side counter can subtend a large solid angle at the scattering 
phosphor without loss of resolution. 


$6. PoLARIZATION—DIRECTION CORRELATIONS IN esr 


The polarization—direction correlations of the 0-90, 1-85 Mev cascade 
in 88Sr were measured with the high energy radiation selected first in the 
polarized and secondly in the unpolarized detectors, corresponding to 
observation of the polarization of the 1-85 and 0-90 Mev transitions 
respectively. The experimental results, shown in fig. 6, determine the 
parities of the 1-85 and 2-75 Mev levels as even and odd respectively, im 
agreement with the results of Bishop and Perez (1955), and with the 
internal pair conversion correlation measurement described earlier. 
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ABSTRACT 


A polarized beam of 100 Mev neutrons is obtained from the A.E.R.E. 
cyclotron ; this is scattered from uranium at sufficiently small angles 
that there is an appreciable scattering from the interaction of the magnetic 
moment of the neutron with the Coulomb field. The scattering and 
polarization of the Coulomb field scattering are in agreement with theory, 
and enable a value for the sign and magnitude of the polarization of the 
neutron beam to be obtained. 


§ 1. INTRODUCTION 


THE scattering of neutrons by a Coulomb field has been investigated 
theoretically by Schwinger (1948). The calculation shows that neutrons 
will be appreciably scattered at small angles and that the scattering is 
spin dependent. Thus a study of the spin dependence enables the 
process to be used as an analyser for a polarized neutron beam. This 
experiment was undertaken, partly to study the process itself, and partly 
to find the magnitude and sign of the polarization of a neutron beam 
which is used in another experiment (Hillman 1955). 


§ 2. METHOD 


A neutron beam which is partially polarized was obtained by allowing 
protons of 170 Mev energy in the A.E.R.E. cyclotron to strike a beryllium 
target. Neutrons were observed at an angle of 26° to the incident 
proton beam. Other targets were also occasionally used. 

The neutrons were collimated by a 6-5 mm steel collimator, and then 
allowed to strike a uranium scatterer. The scattered neutrons were 
detected by a large liquid scintillator (Thresher 1955 b) used as a threshold 

detector. The scattered neutrons were detected on both sides of the 
neutron beam. 

The asymmetry is given by (Cz,—Cp)/(C_+Cp)=PiP2 where P, is 
the (unknown) polarization of the neutron beam and P, is the calculated 
polarization in the magnetic scattermg of the neutrons. Hence by 
measuring C', and ('g—the counts for the counter scattering to the left 
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and right—P, is determined. In this paper the sign of the asymmetry 
is taken to be positive if the scattering is greater in the same direction 
as the initial pn scattering, and negative if in the opposite direction. 


§ 3. APPARATUS 


The counter used was designed to give a good angular resolution. 
It was 44cm longx2:54cem widex7-6cm high. It was viewed from 
the top by 4 E.M.I. 6260 photomultipliers as before (Thresher 1955 b). 

The electronic circuits were the same as used in earlier experiments 
(Thresher 1955 a), and two biases were used on the counter which gave 
two energy thresholds for detection of the neutron beam. 

The effective energies were measured as before ; for the +26° beryllium 
targets they were 104-+-2 mev for the higher bias, and 97-2 mev for the 
lower bias. For the carbon and uranium targets the effective energies 
were not measured but were almost certainly lower. The energy spread 
of the neutrons detected was large; it may be deduced from the data 
of Snowden (1952), and was about +20% at half maximum. 


§ 4. TARGETS 


The targets were placed in a single sight line through the cyclotron. 
Two beryllium targets, the 26,° and 26_°, were placed where this line 
intersected a 49 in. radius, so that the neutrons were produced respectively 
to the left and to the right of the incident protons. The maximum 
bombarding proton energy is 170 Mev, but the mean is about 150 Mev, 
due to target thickness and radial oscillations. Another beryllium 
target is placed so that neutrons from 0° to a 133 Mev proton beam are 
counted. A carbon target was used at +10° (140 Mev bombarding 
energy) and a uranium target at +20°. The main measurements were 
made with the 26,° and 26_° targets. One run was also made with 
a 112 Mev maximum (95 Mev mean) bombarding energy and a 26° Be 
target. The same sight line and collimation were used in all but the last 
experiment. 

The targets were all designed by Dr. P. Hillman to give the best possible 
symmetry of the neutron beam. 


§ 5. BACKGROUNDS 

In spite of great efforts, a large background remained of neutrons 
scattered from the collimators. This was reduced by the use of two 
‘clearing’ collimators which prevented many of these neutrons from 
reaching the counters. The total background was still about } of the 
counting rate from the target and it had to be specially considered. 
The scatterer considerably attenuated the beam, and also attenuated 
the background from the collimator scattering. Part of the background, 
from cosmic rays etc. was not so attenuated. The two could be separated 
by completely stopping the beam. This procedure has been discussed 
before (Thresher 1955, van Zyl 1955). The background was measured 
for every asymmetry cycle, 
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§ 6. ToTrat Cross Section or URANIUM 
The attenuation of the neutron beam by the scatterer was 
experimentally determined, and also yields a value for the total cross 
section for uranium which was necessary for a computation of the 
analysing efficiency. The total cross section becomes 


5-16+-0-06X10-*4 cm? at 106-2 mev effective energy, 
5:16-0:06 10-"4 cm? at 99-+2 Mev effective energy, 
4:90+0-15x10-*4cm? at 73-2 mev effective energy. 


The energy scale is directly related to that of Taylor (1953). The 
results agree well with those of Linlor (1953). 


§ 7. ALIGNMENT OF THE BEAM 
The theoretical cross section varies with angle as A+B/6. It is 
therefore important at the small angles to align the counters with respect 
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The variation of counting rate as the counter is moved through the beam. . 
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to the beam accurately. This has been done by moving the counter 
through the beam (with a reduced beam level) (fig. 1). Although the 
collimation remained fixed, the optical alignment and determination 
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of the effective centre of the targets was uncertain and different for 
each target. The three curves are not, therefore, identical. The curves 
are a fold of the counter efficiency and the beam intensity. It is seen 
that the curves are symmetrical and the midpoint may be determined 
to about +0-2mm. This inaccuracy leads to an uncertainty in the 
asymmetries at }° of 1:2%, at 4° of 0-4% and less at wide angles. 

It was necessary for the experiments to be performed at the same time 
as those of Hillman and Stafford; as a result the collimator had to be 
realigned every few days, and the measured asymmetries include several 
such counter alignments. The error in alignment becomes negligible 
compared with the other errors. 

Checks were also made to show that placing the counter skew to the 
neutron beam by +14° made less than 2% difference to the counting 
rate. Since the maximum accidental misplacement could only be 
+0-1°, no source of error arises. 


§ 8. ANGULAR RESOLUTION 


The curve of fig. 1 also defines the major part of the angular resolution 
of the counter which is about 4° full width at half maximum. Because 
of the height of the counter, there is a small further spread at small angles. 

The experimental points have all been corrected using the theoretical 
variation of cross section and polarization with energy. This correction 
is quite small for the value of the polarization, but is appreciable for the 
value of the unpolarized cross section. Appropriate errors are assigned 
to this procedure and included in the final results. 


§9. SCATTERER 


The scatterer was chosen to be uranium to give a large amount of 
scattering, and a nucleus which was as nearly opaque as possible, so that 
the analysing efficiency should be high. The scatterer was thick to give 
a high intensity and a high ratio of effect to background. 

The thick scatterer involves much multiple scattering. The correction 
for this follows on the lines of Bratenahl (1950). 

Most of the scattering cross section lies within the first diffraction 
minimum. In this region the cross section may be simply represented 
by a Gaussian to a high accuracy. With this approximation, the angular 
distribution of singly scattered neutrons becomes 


o,(0)=exp (—67/a2). 
The angular distribution of doubly scattered neutrons is 
o4(0)=4 exp (—6?/2a2) 
and of neutrons scattered n times is 
o,(8)=(1/n) exp (—6?/na?). 
The total number of doubly scattered neutrons is a times the number 
of singly scattered neutrons where « is the ratio of the scatterer thickness 
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to the mean free path for scattering. The measured cross section is then 


ao 
greater than the singly scattered cross section by the factor Y (x*-1/nn!). 
. ~ * . . tect 
The correction to the asymmetries are less easy ; it is necessary to assume 


that the neutrons entering the counter after multiple scattering are 
unpolarized. This is probable, because the multiple scattering for the 
uranium is principally from addition of two scatters of about 10° where 
the polarization is very small (Dickson 1955). Thus we must multiply 


the measured asymmetries by the factor 2 (x"-4/nn!) which is about 
n=1 
1-13 in our case. 


§10. RESULTS 
The measured asymmetries are summarized in the table. 


Table of Results 


Asymmetries (%) Polarization (°%) of beam 


Target | Angle ————_—_—— 
Low Bias High Bias Low Bias High Bias 


26° Be + —7-9+1-9 | —10-8+1-9 8-9+2-4 12-3+2-4 
26,°Be 4 —5-6+1-5 —7-:7+1-2 7-3 +2-0 10-4+1°8 
26,,°Be 3 —5-9+1-6 —3-2+1-3 —13-2+3-9 —T-4+3-2 
26,°Be | —8-6+2-0 —5:0+1-7 —28 +7-:0 | —17-2+6-0 
26_°Be 4 —7-2+1-7 —6-4+2-] —9-342-5 —8-4+3-1 
26_°Be 1 —2-5+1:8 —13+2-1 —8-2+6-0 —4:2+7-0 
0 Be | 4 | 40-7415| —2-641-5 
02 Be 1 —0-9+2-0 —3-4+3-0 
10° C 4 —0-8+0-8 —0-9+0-7 —1-06+1-0 | —1-22+0-9 
20° U 4 —3-0+2-2 —6-9+3-3 —3-8+2-6 —8-9+4-0 
Weighted average of 26° B =9-22483 | =-10-24-1-3 
95 Mev bombardment :— | 
26,°Be +1-7 | —3-2+1-6 —3-2+2-0 


0-42 | 594414 | 2-34 


Each result is the average of several asymmetries. Standard deviations 
_were calculated in the usual way, and were in agreement with the errors 
expected from the counting statistics. The errors shown are the larger 
_of the two errors. For the 20° U target, however, only 3 asymmetries 
were measured, and counting statistics only are quoted. 

The unpolarized angular distribution is shown in fig. 2, and the 
asymmetries at various angles in fig. 3. The theoretical curves are 
discussed later. The asymmetries from the 0° target are shown in fig. 4 
together with the difference in the asymmetries for the 26,° and the 
26_° targets which, if the experiment is correct, should be zero. Only 
the figures from the higher bias are plotted in the figures, and discussed 
in the text. 
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Small angle differential cross sections averaged over polarization. These are 
compared with the theoretical curve normalized to unity at 1/3°. 


Fig. 3 


5% 
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The asymmetry using a beryllium target at 26°, plotted versus angle, and 
compared with the theoretical curve, normalized to give the best fit 
to the data. 
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The tabulated polarizations are also shown in the table, and include 
the estimated errors of the calculation. The values for the 26° Be targets 
are averaged on the assumption that there are no systematic errors due 
to scattering. The beam as used by Hillman has an effective energy 
of 100 mev, and a best value for the polarization can therefore be deduced 
of 9-8-1-4%. 

Fig. 4 


5% 


nit 


° | 0 
Spurious asymmetries using a beryllium target. 


§ 11. Spurrous ASYMMETRIES 


There are two possible causes of spurious asymmetries apart from the 
counter misalignment. Firstly, the photomultiplier gain could change 
as the counter position was varied ; secondly, there could be scattering 
of neutrons off the walls of the counting area increasing the counts in 
one position. It was unfortunately not possible to arrange for the beam 
to traverse the counting area symmetrically. Both these effects would 
be most marked for the 1° asymmetries, where the counter positions were 
separated by 21 cm. 

To check that these spurious asymmetries were small and to make 
allowances for them, some runs were performed with a target at 0°; 
the asymmetry should be zero. Runs were also performed with the 
+26° targets, so that the neutron beam polarization had the opposite 
sign from each. The difference between the asymmetries so obtained 
gives a measure of the spurious effects, and the average gives the correct 
asymmetry after correction for spurious effect. 
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Figure 4 shows result of both these checks. It is suggestive that 
there might be a spurious asymmetry at 1°. If this were due to a gain 
change, there should be a larger spurious asymmetry for the higher 
bias than for the lower bias. If due to scattering, the scattered neutrons 
would enter the counter at an oblique angle, and there should be a larger 
spurious effect for the lower bias than the higher. No systematic effect 
is observed and it is possible that the results are a statistical accident. 
In any case, the averaging procedure tends to correct for any effect which 
may be present. 

The use of a 0° target and the 26° and 26_° targets does not correct 
for spurious asymmetries due to beam misalignment, which are important 
at small angles and are different for each target. 


§ 12. CALCULATIONS 


Following Schwinger we can define a scattering amplitude 
(f6)+6 .ng(6) 


where f(6@) has real and imaginary parts; and g(@) is totally imaginary. 
Then the scattering cross section for unpolarized neutrons is given 
by o(0)=f?(@)+-g?(6) and the polarization in the second scattering by 


y é)+1 0) 7] 
P(0)= ee ) Im g(6)} 


g(@) is determined both as to sign and magnitude by electrodynamics. 
To determine P,(#), Re f(@) and Im f(@) must be known. 

Over the range of angles used here, we assume that f(@)=/(0)[-/ ,(R6)/6]. 
Also we know, from a well known theorem, Im f(0)=(k/47)oyo4; o¢¢ has 
been measured for this experiment. 

If the nucleus were opaque, Re f(0) would be zero. For a transparent 
nucleus the total cross section varies with energy in such a way that there 
is a maximum, and near this maximum, Re/f(0) is close to zero. The 
maximum for uranium is near 90 Mey (Linlor 1953), so that Re f(0) is small, 
and would be zero if the nucleus had sharp edges. Since nuclei have 
diffuse edges, Re /(0) is not zero. Experimentally Bratenahl (1950) 
at these energies finds that o(0) for lead may be up to 6% higher than that 
given by a square well optical model. Here it has been assumed that the 
nucleus has a constant potential out to a radius R, tapering linearly 
to zero, either by 1:25R or 1:8R. Errors are given corresponding to 
these two extremes. 

We then find by calculation : o(0)=k®o4o¢2/1672 (1-05--.0-03). 

It is, of course, important to know that the Born approximation 
calculations of Schwinger are adequate. Calculations using the potential 
in a small angle approximation (Glauber 1953) show less than 1/10% 
deviation from the Born approximation value. There is, also, a small 
correction due to the failure of the assumptions made by Schwinger 


Scattering of High Energy Neutrons by a Coulomb Field 183 


of a point nucleus. These have been calculated by assuming a uniform 
distribution of charge over the nucleus. At 1°, the largest angle studied, 
the polarization is smaller by 2% and at 1/3° by 0:2%. These are 
negligible. At these large impact parameters no error arises from the 
assumption of a point magnetic moment for the neutron. 


§ 13. SigN or THE POLARIZATION 
The interest of the sign of the polarization is to discover whether the 
sign agrees with the prediction of Fermi (1954) that a nuclear spin-orbit 
coupling responsible for the shell model is also responsible for high energy 
polarization. 
The overall potential predicted by Fermi is 


(1+z2.)U(r)+o.rx(—i) V7 { EAS ae io Un}. 
ie y ar 

Since U(r) is negative, (d/dr) U(r) is positive, and the spin orbit potential 
is the opposite sign to the Coulomb potential. Thus the polarization 
in Coulomb field scattering should be the opposite sign to that in nuclear 
scattering. The asymmetry in this experiment is negative ; the neutron 
_ beam from a charge exchange scattering is polarized in the same direction 
as the proton beams from nuclear scattering, as shown by Bradner (1954) 
and Chamberlain (1954). Thus, the sign of the polarization in nucleon— 
nucleus scattering is indeed in the opposite direction to that of Coulomb 
field scattering, in agreement with the theory of Fermi, and the other 
experiments of Marshall (1955) and Brinkworth (1955). The use of this 
method of obtaining the sign of a polarized neutron beam has been 

suggested independently by Shapiro (1955). 


§ 14. PoLARIZATION IN NEUTRON PRODUCTION 


The polarization in the production of neutrons at 26° by 150 Mev 
protons on Be, is found to be 10-2+1:3%. Roberts et al. (1954) assumed 
that a pn reaction in C at 240 Mev was similar to a np reaction in C at 
150 mev to derive a value for the polarization of their neutron beam. 
By neglecting the energy dependence, they may have obtained too low 
a value of the polarization of their beam by 2°% ; their free np scattering 
polarizations may therefore be too high by 2%. In fig. 5 we have plotted 
the polarization in the pn reaction on beryllium at 26° as a function 
of mean bombarding energy from our data and those of Roberts and 
Bradner (1954). Although these authors did not make measurements 
at 26°, the values may be obtained by interpolation. Also plotted are 
the free np scattering polarizations at 150 Mev (Roberts 1954) and 100 Mev 
(Hillman 1955) similarly obtained. It is seen that there is no discrepancy 
except at the lowest energy. 

_ The polarization of neutrons at 10° from carbon bombarded by 130 Mev 
protons is 1-2+1°% which is lower than expected from free np scattering. 
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On the other hand Roberts and Bradner, both found polarized beams 
from carbon higher than that from beryllium. This low polarization 
at 10° could be explained as follows. 


Fig. 5 
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Polarization in the pn reaction in beryllium at 26° as a function of energy. 
Free np scattering points are included for comparison. 


When neutrons are produced in the forward direction, from a ‘ p, n’ 
collision inside a nucleus, the proton left behind has little energy. The 
resultant !2N nucleus must then be in its ground state. There must 
be a spin change from the 0+ of the ”C nucleus to the 1+ of the 12N 
nucleus. This enhancement of the spin change effect results in a lower 
polarization than in the free np scattering, following an argument by 
Tamor (1954). A similar argument applied to the deuteron has 
been used to explain a low cross section at 0° for the reaction 
N(+D)>P(+N-++N) by Dzelepov (1955). At a lower energy (110 Mev) 
Hoffman (1953) has shown that the cross section for neutron production 
by protons on carbon is nearly flat from 0° to 10° whereas those for 
production by protons on lithium or beryllium rise at 0°, in qualitative 


agreement with this argument. Quantitative agreement is not obtained 
as discussed by Hoffmann. 
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SUMMARY 


Evidence that ductile fracture requires a compressive stress on the 
dislocation arrays close in value to that required for cleavage is presented. 
It is shown that the ductile fracture stress o, is related to the grain 
diameter | by the equation o,=o )+kl- 2, where og and & are constants. 
This is the same form of relationship as alresdy found for cleavage and 
the slope k is practically the same for both types of fracture. 


§ 1. [INTRODUCTION 


Ir now seems fairly certain that cleavage fracture originates at the head 
of an array of dislocations pressed against a grain boundary or other 
obstacle and that the array produces a sufficiently large crack to allow 
this subsequently to propagate in the Griffith manner by its own 
concentration of the applied stress (Koehler 1952, Mott 1953, Petch 1953, 
1954, Stroh 1954, Greenwood and Quarrell 1954). It also seems probable 
that ductile fracture at higher temperatures originates in the same way, 
but that, in this case, the stress concentration produced by the spread 
of the crack activates neighbouring Frank—Read sources and the energy 
of the dislocation array is so dissipated that the crack is arrested. Further 
plastic extension of the specimen is then necessary to complete the 
fracture. 

The mechanism of this completion is of interest. One possibility is 
that the applied stress has to be increased until the initial crack can 
propagate in the Griffith manner against the plastic work that its 
extension entails. But the crack would then propagate as a fast one, 
since the required energy would come from the stored elastic energy, and 
a fast-moving crack is contrary to experience for the initial stages of a 
ductile fracture. It seems more probable that, once a crack fist formed 
ahead of one array, increase in the strain will produce neighbouring 
arrays and neighbouring cracks, and eventually there will be so many 
of these cracks that they will link together into the fracture. The 
subsequent extension of this fracture will then involve further plastic 
deformation ahead of itself, further crack formation and linkage of the 
new cracks with the main one. 


* Communicated by the Author. 
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On this view, the maximum stress required in the fracture process (the 
ductile fracture stress) is determined by the condition for the formation 
of a crack ahead of a dislocation array, but this will not be the first such 
crack that forms, but some subsequent one. The stress required on the 
dislocation array that produces this subsequent crack should not be much 
different from that required for the initial crack. Thus, although there 
are clear differences between cleavage and ductile fracture, the stress 
on the dislocation arrays at the stress for ductile fracture should be 
much the same as would have produced cleavage, if cleavage had been 
possible at the temperature in question. 

To complete the picture of ductile fracture, it can be said that the 
maximum stress should occur at an early stage in the fracture process, 
because of the stress concentration produced by the fracture itself. 
Subsequent to the maximum stress, further strain by the testing 
machine will still be necessary for extension of the fracture. Eventually, 
this will grow to such a size that the work required for its extension can 
come from the stored elastic energy. The fracture will then become a 
fast one, independent of the strain by the machine, and the condition 
for its propagation will resemble that for a Griffith crack, with plastic 
energy added to the surface energy. 


§ 2. EXPERIMENTAL EVIDENCE 


Some measurements have been made that bear on this discussion. 
The figure (a2) shows the dependence of the cleavage strength o, of a 
polycrystalline mild steel on the grain diameter / at —196°c. This is 
for zero reduction of area prior to fracture. The relationship is of the 


form o,=o+hl-V2, 


where o, and k are constants. As discussed previously (Petch 1953), 
this is consistent with the production of cleavage by the stress concen- 
tration ahead of an array of dislocations proportional in length to the 
grain diameter. The term oy represents the stress required simply to 
move a dislocation, so (c¢,—o,) compresses the array. The slope & 
includes the proportionality between the length of the array and the 
grain diameter and is also a measure of the stress required ahead of the 
array. ; 

The figures (b) and (c) refer to similar measurements of the relationship 
between the fracture stress and initial grain size at 18°c and—115°c. 
These are now ductile fractures and there is considerable plastic 
deformation prior to the fracture; on average, 67%, reduction of area 
occurs at 18°c and 61% at —115°c, and the variation with grain size 
is quite small. The ductile fracture stress o, is here taken to be the 
maximum true stress in the tensile test. 

It is apparent that the stress at fracture is given by 


Op=O9+ ki-Ve, 
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which is the same form of relationship as observed for cleavage. Further, 
although the actual magnitude of the fracture stress at constant grain 
size varies considerably between the three cases in the figure, the variation 
resides in the o, term, while the slope & is nearly constant. 

The common form of the relationship between fracture and stress 
suggests that ductile fracture and cleavage are closely similar and that 
the ductile fracture stress is concerned with the production of cracks 

ahead of dislocation arrays that are proportional in length to the grain 
diameter. It is also apparent that (c,—o») at the ductile fracture stress 
is nearly the same as that required for cleavage at the same grain size. 
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The relationship between fracture stress and grain size for a low-carbon steel 
at (a) —196°c, (b) 18°0, (c) —115°e. 


This is what would be expected from § 1, provided that the dislocation 
arrays at constant / have the same length and require the same stress 
ahead of them for the production of cracks at all three temperatures. 

This suggestion of the same array length is rather surprising in view of 
the differences in plastic strain prior to fracture. The dislocation arrays 
might well be considerably shorter in the ductile fractures and, indeed, 
the length might become independent of the grain size. But, it follows 
from the observation that & is constant that the arrays cannot be shorter 
than in cleavage, unless the stress required ahead of them is also less. 
The magnitude of this stress will depend on the rigidity modulus and 
the surface energy, so temperature differences will have only a small 
effect and the plastic deformation prior to ductile fracture, through the 
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surface energy will, if anything, probably cause an increase, not a 
decrease in the stress required. Thus the stress needed ahead of the 
arrays should not be any less in the ductile fractures. It follows that 
shorter arrays are improbable. 

It is concluded that the ductile fracture stress involves dislocation 
arrays of nearly the same size and compressed by nearly the same 
resultant stress as in cleavage at the same grain size. This agrees with 
the ideas in $ 1. 

The lack of an effect of plastic deformation on the size of the arrays 
may indicate that the dislocations in iron are only successfully held up 
by grain boundaries, or, alternatively, that the size appears unchanged 
because a number of small arrays produce the same effects as a large 
one. 

The variation of oy in the figure has still to be considered, but this seems 
readily understandable. Presumably, it represents the change in the 
stress required to move a dislocation as the temperature is lowered and 
the iron strain-hardens. The big difference in the strain prior to fracture 
could explain the smaller value of op at —196°c than at —115°c. 

Since oy varies with strain, there may be some difference between the 
ductile fracture stress and the applied stress necessary to form the first 
crack, although there is little difference in the resultant stress on the 
dislocation arrays. 


§ 3. EXPERIMENTAL DETAILS 
The mild steel used had the following percentage composition : 


C Si S P. Mn Ni Cr N 
0-15 0:02 0:05 0-03 0-51 0-08 0:05 0-008 


The fracture stress was determined on specimens of | in. gauge length 
and 0-125in. gauge-length diameter, polished to an 000 emery finish. 
These were machined from in. diameter cold-drawn rod. The grain 
size range of 80 to 10000 grains/mm? was principally obtained by 
holding the specimens in vacuo for various times at 900-1200°c and 
cooling at various rates. The faster rates were all followed by a short 
anneal at 650°c with subsequent slow cooling. Straining and annealing 
were used to obtain the very coarse sizes. 

In the tensile testing, special care was taken with the axiality of 
loading and the machine was driven at a cross-head velocity of 
0-34 mm/min. , 

The ductile fracture stress was identified with the value at which there 
was a sudden discontinuous drop in load, followed almost immediately 
by fracture. For the specimens that necked, the lack of uniformity of 
stress had to be considered. Since we are concerned with the stress 
compressing the dislocation arrays, the measured average tensile stress 
in the neck was converted by Bridgman’s (1943) analysis into the tensile 
stress that would have been present if uniaxial conditions had persisted. 
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The required measurements of the diameter of the fracture and the radius 
of curvature of the neck were made with a projection microscope. 

Within each set of measurements there was some variation of the 
strain at fracture. This produces a variation in oo, so the fracture stress 
figures were corrected to a constant strain at each temperature. 

For —196°c, this correction was obtained from a separate experiment 
in which the effect of strain on the fracture stress was estimated by 
pre-straining at room temperature and then fracturing at —196°c. This 
estimate was in close agreement with another made in a different way. 
Up to grain sizes such that [-2~5-5 mm~?, there is very little strain 
prior to fracture, then the strain increases, but again becomes fairly 
constant for the finer grains (Petch 1953). By plotting the two 
constant-strain portions of the curve of o, against [~1/? without any 
correction for strain, a determination of the variation of og with strain 
can be made. 

At 18°c and —115°c a slightly different approach was possible, as 
described elsewhere (Petch 1955). 

Any uncertainty in these corrections does not affect the main con- 
clusion from these measurements, namely that the slope k is nearly 
independent of temperature. The slope at —196°c is largely fixed by 
measurements that involve very little strain before fracture, so these do 
not require any correction. At the other two temperatures, there is little 
variation of the strain at fracture with grain size, so the corrections 
required are very small and could be neglected. Also the allowance for 
necking by the Bridgman analysis is practically constant throughout 
the grain size range at these temperatures, so any uncertainties in this 
allowance have little effect on the slope. 
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ABSTRACT 


The scattering of phonons by electrons is calculated, assuming the 
usual electron-phonon interaction, for a parabolic band whose degeneracy 
temperature is comparable: with the temperature of the lattice. The 
contribution to the thermal resistance is given by an exact formula, 
subject only to justifiable assumptions concerning phonon—phonon 
interactions. With rising temperature the apparent mean free path 
of the phonons at first decreases as 1/T (or, if there are very few electrons, 
as exp (a/T')), but reaches a minimum and then increases as 7%. Energy 
and momentum conservation then allow only the tail of the electron 
distribution to contribute to the scattering. The model is thought to 
apply to certain observations on p-type germanium, irradiated sapphire and 
conducting diamond. 


Tue thermal conductivity of a dielectric solid at low temperatures is 
conveniently expressed by the equation 


k=+tevl, ee a ee. LY 


where the specific heat, c, is proportional to 7° at low temperatures, 
while the phonon velocity, v, may be taken as constant. Observation 
of « thus gives the behaviour, with varying temperature, of /, the phonon 
‘mean free path’. Without going too deeply into the meaning of this 
quantity, we may see that, as the phonon distribution at temperature 7’ 
is heaped up around a fairly high peak at the modes of angular frequency 


vr3hTh, Pagettece seme a (2) 


the variation of / as a function of 7’ represents, crudely, the variation of 
the mean free path as a function of vy. For example, if « is proportional 
to T’-1 we conclude that the mean free path of a phonon is proportional 
to the fourth power of its wavelength—the result characteristic of 
Rayleigh scattering by obstacles small compared with the wavelength. 
This argument, combined with a variant of Matthiessen’s rule for adding 
thermal resistances due to distinct scattering processes, can yield useful 
qualitative information about the types of lattice imperfections present in 
a given specimen, although it is difficult to justify it in detail (Berman 


1953). 
: * Communicated by the Author. 
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If one thinks of the various possible imperfections—vacancies, inter- 
stitials, dislocations, impurities, etc.—one can see that / should always 
increase as the phonon wavelength increases. Only total scattering by 
the surfaces of the specimen, or by large-angle grain boundaries, could be 
expected to be independent of the phonon wavelength ; even then, the 
effects of specular reflection cannot be neglected (Berman, Simon and 
Ziman 1953, Berman, Foster and Ziman 1955). However, recent 
experiments by Berman. Schneidmesser and Tirmizi (1955) have shown 
that in certain materials, such as gamma-irradiated artificial sapphire 
crystals and conducting diamonds, / can be nearly independent of 
temperature and yet very much smaller than can be explained by 
boundary scattering alone. Moreover, Rosenberg (1955) has found that 
in a certain specimen of p-type germanium the thermal conductivity 
varies as 7%, which is apparently equivalent to / decreasing as the 
phonon wavelength increases. 

To explain these anomalies, it was suggested that the phonons might 
be scattered by electrons (or holes) free to move in some sort of band 
which is not full at very low temperatures. As will be shown, a simple 
calculation along well worn tracks leads to unsuspected conclusions when 
the combined effects of momentum conservation and Fermi degeneracy 
are allowed for. This calculation is reported here ; the detailed interpre- 
tation of the experimental results will be presented at a later date. 

For the sake of simplicity, we take the most elementary model, in which 
the phonons have a ‘ Debye’ spectrum with a single constant speed »v, 
independent of their direction and polarization. The electrons will be 
supposed to lie in a parabolic band with effective mass m*. By ignoring 
anisotropy and dispersion we have the advantage that the equations may 
be solved exactly and involve only these two parameters. 

The calculation is a straightforward adaptation of the standard theory 
of the phonon-electron interaction and has been given already in the 
limiting case of a typical metal (Sommerfeld and Bethe 1933, Makinson 
1938, Wilson 1953): we here consider the slightly more complex case 
where the temperature is comparable with the Fermi degeneracy tempera- 
ture of the electron distribution. In brief, we are concerned with processes 
by which a phonon, of wave vector p, interacts with an electron in the 
state of wave vector k, which is then scattered into the state k’ whilst 
the phonon is absorbed. The rate at which this process occurs is given by 


= — Fy (1— my) Q(B’ —B), ae di 3G) 


where (a) = sin (at/h)/(ajh). EH, EH’ ave the initial and final energies of 
the whole electron-phonon system ; N, is the number of phonons already 
in the mode p; i, 2,’ are the numbers of electrons in states k, k’ ; 
p is the mass density of the lattice whose total volume is V. The 
parameter C has the dimensions of energy (a few electron volts, say) and 
measures the interaction between electrons and phonons. 
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Crystal momentum will be conserved in the transition (we are far below 
the possible limits for Umklapp processes). Hence 


k’=k-+p. Secsdeve tae ae (4) 
.For conservation of energy, we have 
el ies hk? 


Putting these relations into (3), and summing over all possible electronic 
states, we find that the total rate of change of N p due to scattering by 
carriers takes the form 

; oN, me O2m* 2 ; 
at =—L Papo | m(1—n, )k dk, a ¢ : (6) 


where ky is the minimum value of the modulus of k which will allow 
both (4) and (5) to be satisfied, i.e. 


ko= | 4p—m*ofh |. ware ere reue 
If for n,, 7’, we assume the usual Fermi functions, e.g. 
n =[exp {(L’—C)/kT}4+-1})-2 a BES) 


then (6) may be integrated exactly in terms of elementary functions. 
Since (6) is a differential equation representing the exponential decay of 
the number of phonons in the pth mode, we have at once the relaxation 
time for that mode, i.e. : 


Meecm ty [ kT | fz 
ie Ah'*p dm*y2 i= 


al l+-exp {(4m*v?—€C)/kT+22kT'/8m*v2+-32z}\ , 
x 1-2 epost 
where z=fv,/kT. 

At this point we run into a difficulty which lies deep in the theory of 
phonon transport and which has not been successfully resolved : how are 
we to calculate the contribution of this additional scattering process to 
the thermal resistance of the crystal 2? One method would be to write 
down a mean free path 

iB ; (10) 
for each mode, from (9), and to sum over the phonon distribution a series 
of contributions of the form of (1). Thus we would write 


c=} cprl, dp, ewe careee a 2 ar (11) 
where ¢, represents the contribution of the pth mode to the specific heat, 
ie 


4 UTy ; « 


c= [cy dp. ee ee ae ee 2) 


Unfortunately (11) diverges. ‘To obtain a finite answer it would be 
necessary to take, in place of /,, the mean free path for the combined effect 
of boundary and electron scattering, e.g. 

LaLa me le Pactra op 


02 
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where D is the diameter of the specimen. But the boundary and 

electron-scattering resistances would then no longer be additive, and 

the application of Matthiessen’s rule would be spurious. This method 

could be made to yield an answer, but the calculation would be very 

complicated because of the clumsy form of (13). . 
The alternative procedure* is to calculate the following average 


ay baer 

==-|dp. ee ee ee 

[=</T? (14) 
‘This integral always converges, since |, never tends strictly to zero. If 
LON OES W.=[seel], vs eth tenes 


then this is a thermal resistance which is proportional to the integral in (14), 
depends only on the electron scattering process, and may legitimately be 
added to other thermal resistances to give the total resistance of the 
specimen, according to Matthiessen’s rule. 

An intuitive argument may be given to justify (14). If we use (11) we 
are, in effect, assuming that the phonon modes are independent of one 
another, and carry their separate contributions to the heat current without 
interference. This assumption is certainly false. There are, between the 
modes, anharmonic coupling terms which cause three-phonon processes, 
where a given phonon splits into two others, or two combine to give a 
third (Peierls 1955). The Umklapp processes, which give rise to the 
‘intrinsic’ thermal resistance of the solid at high temperatures, are 
certainly frozen out of this calculation ; there are still the other processes, 
in which momentum is conserved, by which phonons interact strongly, 
but which give no direct contribution to the thermal resistance. : 

Now, suppose we have a packet of energy moving through the lattice, 
concentrated into modes of nearly a single wavelength. As a result of 
the phonon—phonon interaction, the initial phonons will soon have changed 
into a variety of others of different wavelengths. There will be a continual 
exchange of the energy in the packet from one mode to another. If this 
exchange takes place in a distance short compared with the distance 
between scattering electrons, the characteristic relaxation time for the 
packet cannot depend on the wavelength with which it started. Instead, 
all different values of the electron scattering cross section will be sampled, 
over all different wavelengths, to yield an average such as (14). 

The mathematical justification derives from the general variational 
principle of transport theory, as applied to the phonon-phonon processes. 
It can be shown (Ziman, unpublished) that W, is an upper limit to the 


excess resistance induced by scattering processes with relaxation time 7 
for phonons in the pth mode. It is assumed that the steady state phonon 
distribution is dominated by the non-Umklapp processes, but the problem 
of keeping the long-wave longitudinal phonons in equilibrium (Peierls 


1955, Herring 1954) is not considered, and may lead to some error. 


* This procedure was used by Berman, Foster and Ziman (1955), on intuitive 
grounds, before a mathematical justification was found. 
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The final result of our calculation comes from (9), (10), (12) and (14): 
the extra thermal resistance due to the scattering of phonons by electrons 
should appear to be due to a process for which the mean free path at 
temperature 7' may be written 


Lele Lgl Tati’), Site oe (16) 
where 
lp —=4rhi4p/C2m*?, id reason lz) 
kT 5=4m*v?, meee ee LS) 
1_T 15 ;°exp (3) 
a3 he sa | 9 (sinh 42)8 


{ exp (7T'5/T —C/kT +227 /16T'5)-+exp (42) ; 
xin ET LTP ee dz. vart19) 


The integration can be carried out numerically with reasonable accuracy. 
To find ¢ as an explicit function of temperature one must allow for the 
variation of the Fermi energy € with 7’. This may be taken from standard 
tables (e.g. Wilson 1953) once one has decided the degeneracy temperature 
of the electron distribution, i.e. 


T p= (h?/2m*k)(302n)? ee ee tO) 
if there are n electrons per unit volume. 


The dimensionless function ¢ has been calculated in terms of 7'/7's for 
various values of 7',/7's. The curves plotted in fig. 1 thus give the expected 


Fig. 1 
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behaviour of the mean free path, suitably scaled, as a function of 
temperature relative to the scale temperature 7's, for different values of 
the free electron density. To grasp the behaviour of this function, let us 
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consider its asymptotic forms. Simple calculations from (19) give the 
following : 

(a) Toa gas lo f>0-20457',/T'. ) Saas 

This is the usual result calculated for a metal. At low temperatures 
the electron distribution is highly degenerate. Only electrons in an 
energy range of width k7' at the Fermi surface can take up the phonon 
energy and be scattered into empty states. The mean free path comes 
out inversely proportional to this number, i.e. varies as 1/7’. 

(bs) Ty<Ty; a8 T>0, $>0-0336(7'5/T) exp {((T's—T p)/T}. (22) 

The lower limit, ky, of electron wave vectors is now operating. According 
to (7), when p is small, as at very low temperatures, only electrons with 
momentum exceeding m*v can scatter the phonons. From (18), these 
have energy k7'y. But, when 7',<7's, there are so few electrons that 
when the distribution is degenerate the Fermi level lies below k7's. The 
number of electrons available for scattering is now reduced by a Boltzmann 
factor for the energy difference between the Fermi level, k7',, and the 


level k7', ; the mean free path increases exponentially as the temperature 
falls.* 


pe a Efe \ 8128 Ne 4 T's; 1/2 |—-1 
eeeeitar dc cry) acy) || 


This is the most interesting portion of the curve. If 7’>7',, then most 
of the phonons have momenta far exceeding m*v, and the value of ky 
becomes nearly equal to }p. Let us fix our attention on the ‘ dominant ’ 
phonons at temperature 7’, i.e. those of energy around 347’. At some 
temperature the cut-off limit 4, for these might lie as in fig. 2(a), so that 
most of the electrons can interact with the phonons. But, at twice the 
temperature, this limit will have moved out, as in fig. 2(b), so that only 
about half the electrons can act. If the temperature is doubled again, 
although the electron distribution broadens out and spreads to higher 
energies the limit is so raised that only the tail lies above it. It is this 
decrease in the number of electrons available for scattering that gives 
the rise in the mean free path. The asymptotic formula (23) show ¢ 
proportional to 7. 

Rosenberg’s observations could thus be explained if he were working 
at temperatures in the region (c), i.e. where the mean free path would 
appear to decrease as the temperature falls. Preliminary quantitative 
calculations show that 7's and /) might have values which would place the 
measurements on germanium on this part of the curve. In sapphire, 
and in diamond, the higher velocities of sound would make 7’, much 


* Professor Gorter has suggested, in private discussion, that the 2nd order 
perturbation terms (virtual absorption and re-emission of a phonon by an 
electron) might be important, since the momentum limit would not then apply. 
Inspection of the matrix elements shows that the probability of such a process 


is down by several orders of magnitude below the first order terms, and hence 
would not be observable. 
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greater, so that we should be in the region of the minimum in ¢ at helium 
temperatures. This might well give apparently constant mean free path 
over a restricted range of temperatures, as in Berman’s experiments. 


Fig. 2 


Ko Ka 


Electron distribution in momentum. Electrons capable of scattering the 
phonons of energy 3k7' lie in the shaded regions. 


(a) T/Ts=4; (b) T/Tg=8; (c) T/Ts=16. 


It may be remarked that in this region the mean free path depends mainly 
on the electron density, and is almost independent of the value chosen 

’ for m*. It is not to the point here to discuss how it is possible that there 
should be free electrons in these materials at such low temperatures, 
although the answer probably lies in the presence of impurity bands. 
We have seen, however, that from a grossly over-simplified model one can 
deduce precise formulae which seem to represent, at least qualitatively, 
the experimental results. 
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ABSTRACT 


The decay of °K has been reinvestigated and values of the branching 
ratio, defined as the ratio of gamma or electron capture to beta transitions, 
have been determined by two methods, giving the values 0-124-+0-002 
and 0-121+-0-004 respectively. The absolute beta and gamma decay 
rates and the half-life have been measured, the value of the latter being 
1-28+-0-02 x 10° years. A fresh determination of the gamma energy is 
described, and evidence as to the mode of decay is discussed. Finally, 
the results are given of experiments on the reflection of electrons by 
solid source supports. 


§ 1. INTRODUCTION 

THE decay of Potassium 40 has been the subject of many investigations 
in the past and it is now fairly well established that the decay is according 
to the scheme shown in fig. 1, first suggested in a slightly modified form 
by Suess, 1948. The experiments about to be described were undertaken 
to elucidate several uncertain results, in particular the branching ratio 
of the decay which will be defined here as the ratio of the gamma or 
electron capture transition rate to the beta transition rate. This ratio 
may be determined by two methods :— 

(1) By direct counting methods and the application of suitable 
corrections for solid angle, efficiency etc. 

(2) By measurement of the amount of argon or calcium, or both, in 
a potassium-bearing mineral of known geological age (the geological 
method). 

The results obtained by method (1) show a certain variation, but in 
general tend towards a value of the ratio of 0-12. (Sawyer and 
Wiedenbeck 1949, 1950, Houtermans et al. 1950, Ceccarelli et al. 1950, 
Spiers 1950, and others.) The results of some investigators are not 
quite independent since they measured only the gamma or beta activity 
and relied on the results of others to derive the ratio. The results obtained 
by method (2) show a greater variation and the values obtained range 
from 0-05 to 0:12, more recent results being in the region of 0-09 (Shillibeer 


and Russell 1954, Wasserburg and Hayden 1954) to 0-12 (Inghram et al. 


es eo a eee eee 
* Communicated by the Authors. 
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1950). ‘The fact that the geological method tends to give a lower value 
is of considerable importance since the reverse process of measuring 
the K—4°A ratio is of value in the determination of the age of minerals. 
Two reasons for arriving at a lower ratio which one might suspect are. 
firstly, loss of argon from the mineral during geological times and. 
secondly, incomplete removal of argon from the mineral sample. That 
the latter can happen has been proved by re-examination of the minerals 
used by Russell et al. (who got a value of 0-06, Russell e¢ al. 1953) by 
Wasserburg and Hayden (1954) and by Shillibeer ef a/. (1954) the results 
of the last two investigations yielding a value of about 0-090 due to more 
efficient argon extraction. In view of the discrepancies, it was felt 
well worthwhile making, as carefully as possible, another determination 
of the branching ratio by the counting method. At the same time fresh 
estimations were made of the absolute beta and gamma decay rates 
(and, hence, of the half-life) and the energy of the gamma transition 
while further experimental evidence for the decay scheme of fig. 1 was 
sought. In order to arrive at absolute values of the beta decay rates, 
experiments on the back-scattering of electrons were made, and are 
reported here. 


Fig. | 


A+° K 40 


R- 1-35 MEV 


U 1:46 MEV 


(c +) 


(o +) 


Decay scheme of #?K. 


§ 2. THe MEASUREMENT OF THE BRANCHING RatTro 
T'wo methods were employed to measure the branching ratio. In 
both methods, it was assumed that the decay scheme of fig. 1 was correct 


and that the ratio of the absolute rates of emission of gamma rays to 
beta particles yielded the branching ratio. 


First Method 


In the first method, difficult and uncertain calculation of the beta and 
gamma detector overall efficiencies was avoided, to a first order, by 
comparing the known gamma-—beta ratio of Co to that of “°K. The 


energies of the gamma transitions of %Co (1-17 and 1-33 Mev) are fairly 


close to the "K gamma energy (1:46 Mey). Had the gamma ray energies 
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the maximum beta ray energies and the beta spectrum shapes been 
identical, it is easy to show that the branching ratio for "K, R say, would 
have been given by i 


R=2 YK Boo 
Bx Yeo 
where yx, Bx are the gamma and beta counting rates of a known ratio 
of amounts of K and yg, Bg. are the corresponding rates for the same 
ratio of quantities of ®Co. As the gamma and beta spectra were not 
identical corrections had to be made. These will be described later, 

The gamma detector is shown in fig. 2. A Nal crystal, thallium 
activated, 2in. long and ? in. diameter, is surrounded coaxially bv 
a cylindrical vessel of capacity 186 cm? which was filled to a mark A 
with the solutions which were measured. The crystal, which was 
surrounded by sufficient absorber to absorb all beta particles, was mounted 


Fig. 2 
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EMI! TYPE 6262 
PHOTOMULTIPLIER 


Gamma ray detection apparatus. Lead shielding not shown. 


on an EMI type 6262 photomultiplier and the whole encased by 4 in. of 
lead. The pulses were amplified, passed through a discriminator and 
recorded by a scaler. The solutions used were, (i) water containing 
a known aliquot of a standard °CoCl, solution and (ii) an aqueous 
solution of potassium hydroxide (493 g K/litre), the latter giving 700 ¢.p.m. 
above background (87 ¢.p.m.). It was found that the background with 
water, carbon tetrachloride or aqueous solutions of calcium chloride 
or sodium chloride in the vessel was slightly higher than when it was 
empty, the actual rate being independent, within experimental error, 
of which liquid was employed. Hence the background rate was always 
checked with water in the vessel. Integral bias curves were obtained 
for the test solutions between the energies of 10-100 key. Calibration 


202 A. McNair, R. N. Glover and H. W. Wilson on the 


was effected by the 46-7kev gamma ray of RaD. The results for 
potassium with and without background subtracted are shown in fig. 3 (a) 
and fig. 3(b). Figure 3(c) is the curve obtained after correcting for 
bremsstrahlung as will be described later. Figure 4 is the corresponding 
curve for Co. The counting rates yx. yo, to be substituted in the formula 
were obtained by extrapolating the curves to zero energy and making 
the necessary corrections to the values thus obtained. 


Fig. 3 
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Discriminator bias curves for 4K gamma rays. 


The beta detector was a cylindrical, end-corrected proportional counter 
filled with argon and methane to a total pressure of one atmosphere. 
The source was spread uniformly over an area of ~1000 cm? on a thick 
stainless steel foil pressed against the wall of the counter. The large 
surface area ensured reasonable counting rates from thin sources of 
*K(~1 mg/cm?). In order to reduce the background, the counter 
was surrounded by a double layer of G—M counters in anti-coincidence 
and the whole, in turn, by 4 in. of lead on all sides. The anti-coincidence 
arrangement was not used in the measurement of the ®°Co beta emission 
rate as a reduction in counting rate of ~4°% was found due to coincident 


detection of the gammas in cascade with the betas. Integral bias curves 
were taken from 0:2—10-0 key. 


§$ 3. CORRECTIONS 
In calculating the branching ratio from these results corrections must 
be made for (i) absorption and scattering of 4K gamma rays by the 
dense and relatively high Z potassium hydroxide solution ; (ii) the slight 
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difference in the crystal efficiency for °K and Co gammas ; (iii) gamma 
detector counts due to bremsstrahlung, especially in the case of °K - 
(iv) coincident detection in the crystal of both “69('o gamma rays. A 
separate series of experiments, described later, showed that the fraction 
of electrons reflected by the stainless steel source support was the same 


an , : 
for Co and *K so that no correction for this effect was necessary in 
determining the branching ratio. 
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Discriminator bias curve for ®Co gamma rays. 


The correction for absorption and scattermg of *K gamma radiation 
was determined by adding a known small volume of ®CoCl, solution 
to the potassium hydroxide solution and comparing the increase in 
counting rate with that due to the same quantity of ®CoCl, solution 
diluted with water to fill the vessel. It was found that the effect of the 
potassium hydroxide solution was to increase the counting rate slightly 
showing that the absorption of the gamma rays was more than com- 
pensated for by their being scattered into the crystal, in the geometry 
employed. A correction to the gamma counting rate of about 6% had 
to be applied to allow for this effect. 

The correction for the slight difference in detection efficiency for ®°Co 
and “K gammas was determined by postulating an average path length 
in the crystal, X say, which the gammas have to traverse. Then it 
can be shown that the ratio of the true count rates, [g,, 7x say, is given, 
for our experimental arrangement, by the expression 


Too 1.915 P& Yeo 
Ix Poo YK 
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where px, poo are the absorption coefficients of Nal for the #K and Co 
gamma rays and yx, yc, are the measured count rates. . 

It is necessary to make a correction to the *K gamma counting rate 
for the effect of bremsstrahlung. This effect is much more important 
for the case of “°K than for ®Co since (a) the ratio of betas to gammas 
is ~16 times greater, (b) the ®Co beta spectrum is softer, (c) the Co 
betas pass through a solution of lower density and of lower atomic number. 
The effect was estimated by examining, as a function of discriminator 
voltage, the counting rate in the gamma detector due to the addition 
of a solution of 32P which emits only electrons. The strength of the 
solution was monitored by measuring the beta rays, using the pro- 
portional counter. Figure 3(c) shows the 4°K spectrum corrected for 
bremsstrahlung. 

A small correction was made to the Co gamma rate to allow for the 
probability of coincident detection of the two gammas which are emitted 
in cascade by this source. 


Result by the First Method 


The value of the branching ratio obtained by this method was 
R=0-124+0-002 


or 11:0+0-2% of the decays are to A, 89-0-+0-2% to Ca. 


Second Method 


This method eliminated many of the corrections necessary in the first. 
It made use of the fact that ?4Na possesses a gamma transition very close 
in energy to that of *°K (1-38 Mev and 1-46 Mev respectively) while the 
beta ray energies are almost identical (1-39 Mev and 1-36 Mey resepctively). 
In principle, the areas under the photo-electric peaks of the gamma rays 
as observed by a large Nal(Tl) crystal were compared, the sources being 
monitored by their beta emission. *4Na emits also a 2-76 Mev gamma ray 
which precluded its use in the first experiment. 

The gamma spectrum of *K was examined by placing a shallow 
cylindrical glass vessel (1 in. highx2in. diameter), filled with a con- 
centrated solution of potassium hydroxide, on top of a large 
(2 in. long x 2 in. diameter) Nal(Tl) crystal which was mounted, in turn, 
on a Dumont Type 6292 photomultiplier, the whole being surrounded by 
lead. The pulses were analysed by means of a 100-channel kicksorter 
of the Hutchinson-Scarrott type. A very small quantity of 24NaCl 
solution was mixed with the potassium hydroxide solution and a second 
spectrum taken. By subtraction and knowing the form of the background 
spectrum, the spectra, corrected for background, were obtained (figs. 5 
and 6). The areas under the photo-electric peaks, corrected in the case 
of *4Na for the small contribution from the Compton electrons due to the 
2-76 Mev gamma ray, gave a relative measurement of the gamma ray 
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intensities. The beta rates were found for the potassium hydroxide 
solution and for the mixed solution, by examination in the proportional 
counter used previously. The result was corrected for the decay of the 


24Na. 
Figs. 5 and 6 
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The gamma ray spectra of #°K and ®Co. 


The branching ratio is now given by the expression 


where yx, Px and yyq; Pa refer to the gamma and beta counting rates 
of potassium and *4Na respectively. 

As the beta ray energies were practically the same, no correction for 
‘reflection of electrons by the source and source support was necessary. 
(This experiment was devised before completion of the reflection experi- 
ments which showed no change in the back scattering factor over an 
energy range which included °Co and 40K beta ray energies). Since 
the 24NaCl was added to the potassium hydroxide solution, no correction 
for absorption and scattering of “°K was necessary. The use of a mixed 
source in the proportional counter meant that no correction for source 
thickness was necessary since it was effectively the same for both the 
4K and the 24Na measurements. Bremsstrahlung from the beta rays 
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did not affect the magnitude of the peaks since the maximum beta ray 
energies are close to the gamma ray energies. 

Corrections were made for (a) the probability of coincident detection 
of the 1-38 and the 2-76 Mev gamma rays in the 74Na measurement, 
(b) the slight difference in the scintillation detector efficiency for 1-38 Mev 
and 1:46 Mev gamma rays and (c) the dead time of the kicksorter, which 
was 0:75 milliseconds. The first was calculated by estimating the actual 
disintegration rate from a knowledge of the beta counting rate and 
assuming a backscattering coefficient of 0-45 for the stainless steel support. 
The efficiency of the crystal for 1-38 and for 2:76 Mev gamma rays was 
calculated and the error introduced due to coincident detection found 
to be ~7%. Note that the assumption of a particular value of a back 
scattering factor here does not introduce any appreciable error, since the 
value of the correction is small. The second correction was calculated 
as described earlier for the first method. 


Result by the Second Method 
The value of the branching ratio is 
R=0-121-40-004 
or 10-8+0-4% of the decays are to *°A, 89-2+0-4% to *°Ca, in good 


/O 
agreement with the value obtained by the first method. 


The Absolute Decay Rates and the Half-life of *°IK 


Assuming a backscattering factor of 0-450-+0-015 for °K beta rays, 
we calculate the disintegration rate for the decay to #°Ca to be 


1650-15 betas/gram/minute 
and for the decay to 4°A to be 
204-5-++-3°8 gammas/gram/minute. 


This figure, of course, does not include argon x-rays resulting from 
electron capture. These values give a total half-life of 
1-28+ 0-02 x 108 years 


assuming Nier’s value of 0-0119-+-0-001% for the abundance of the °K 
isotope in potassium (Nier 1950). 


§ 4. Tor ENERGY oF THE GamMa Ray 


The gamma ray spectrum emitted by ~330 g of pure potassium chloride 
was examined using the crystal and kicksorter which were employed 
in the measurement of the branching ratio by the second method. 
Comparison with the gamma radiation of 22Na (1-277 Mev) and 24Na 
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(1-38 Mev) gave an energy for the K gamma ray of 1-46--0-01 mev in 


agreement with values obtained by other investigators (Bell and Cassidy 
1950 b, Good 1951 a). The gamma spectrum is shown in fig. 7. 
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§ 5. Tor EvIpENCE FOR THE MODE or Decay or ®K 


There is a large body of evidence favouring the decay scheme of fig. 1, 
although much of it is not altogether convincing, due to the difficulty 
of carrying out coincidence experiments with the very low specific 
activities involved. Perhaps the most satisfying proof is that of Johnson 
(1952) who showed by measurement of the masses of A, *K and Ca 
that the energies available for the decay to *A and to Ca are 
1-30+0-07 Mev and 1-49-+0-07 Mev respectively. The first value is in 
good agreement with the maximum beta energy (1:35 Mev) and examina- 
tion of the spectrum shape fails to reveal the presence of softer beta 
spectra. Also the energy available is less than the energy of the gamma 
ray which must therefore be associated with the decay to A. Sawyer 
and Wiedenbeck (1950) report that the electron capture disintegration 
rate approximately equals the gamma emission rate. Thus, the probab- 
ility of electron capture to the ground state of 4°A is fairly small, a result 
supported theoretically by Morrison (1951). The absence of 0-511 Mev 
gamma rays (Bell and Cassidy 1950a, 1950b, Good 1951 b) and of 
180° coincidences (Colgate 1951) show that the probability of decay 
by positron emission is very small, <0-06% according to Colgate and 
<0-:10% according to Good. Morrison’s theoretical investigation 
confirms these results. 
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An examination of fig. 7 shows no evidence of a lower energy gamma, 
transition which shows that decay to an excited state of 4#°A below 1-46 Mev 
is not probable ; nor has any such level been found by particle excitation 
experiments (Endt and Kluyver 1954). 

Johnson’s results make it most unlikely that there should be two 
gamma transitions in cascade unless one were very soft. Nevertheless 
direct experimental evidence is desirable. Figure 8 (a) shows the *®K 


Fig. 8 
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(a) *K gamma spectrum above 1-46 Mev. 
(b) ®Co gamma spectrum above 1-33 Mey. 


gamma spectrum above 1:46 Mey while fig. 8 (b) shows the corresponding 
region for ®°Co. In this case there is a fairly well defined peak A at 
2-50 Mev due to the coincident detection of the 1-17 and 1-33 Mev gamma 
rays. B represents the Compton edge of one gamma ray carried forward 
by the coincident photo-electric capture of the other. The region C 
is due to the integration of the two Compton distributions. No such 
effect is observed in fig. 8 (a) and from this graph we believe that fewer 
than 10° of the gamma transitions are in coincidence. A separate 
gamma-gamima coincidence experiment was carried out to investigate 
this point further. A Nal detector was placed against the side of the 
vessel shown in fig. 2 and the coincidence rate observed as the vessel 
was filled in turn with water, potassium hydroxide solution and CoC, 
solution. However, the overall efficiency of the system was very low 
‘and, on the basis of this experiment, we can say only that fewer than one 
in five gamma rays were in coincidence. 

Although Johnson’s figures do not permit of the possibility of a gamma 
of energy equal to, or greater than, 1-35 Mey being involved in the decay 
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to “Ca, and Endt and Kluyver report no energy level in *°Ca lower than 
3-35 Mev, several searches for beta-gamma coincidences have been made 
(Meyer et al. 1947, Houtermans et al. 1950). No beta-gamma coincidences 
were detected. We also carried out a search for beta-gamma coincidences 
using a source of enriched potassium. Enriched potassium chloride 
(60 mg) was placed on a 1 mil aluminium foil and covered with a thin 
plastic sheet and a 50 mil sheet of steel to absorb the beta rays. On top 
of this was placed a lin. cube of Nal to detect the gammas while on 
the other side the beta detector consisted of a 23 in. diameter plastic 
phosphor (1:5% 1:1:4:4 tetraphenyl 1:3 butadiene dissolved in 
polystyrene). The whole system was surrounded by 4in. of lead. The 
pulses were amplified, passed through discriminators and then 
into a coincidence circuit, the resolving time of the circuit being 
6-2 1077 seconds. The beta, gamma and beta-gamma_ coincidence 
rates were recorded as a function of beta discriminator voltage. The 
enriched potassium chloride was found to be only ~0-57°% #°K whereas 
a degree of enrichment of approximately 150 had been anticipated. 
Consequently the counting rates were lower than had been expected 
when the experiment was planned. With the source in position the 
coincidence rate increased by 1-43-+-0-20 10-3 coincidences per beta 
detected. Since some of this could have been due to bremsstrahlung, 
32P was substituted and gave an increase of 1-040--0-016x 10-3 
coincidences per beta, i.e. a somewhat smaller increase than in the case 
of °K. However, in view of the low counting rates involved the authors 
do not attach too much importance to this difference. It is hoped that, 
at a future date, the use of larger quantities of more highly enriched 
potassium and bigger crystals will permit more significant measurements 
to be made. 


§ 6. THE BACKSCATTERING OF ELECTRONS 


Two methods of measuring backscattermg were employed. In the 
first, a thin beta emitting source was deposited on a thin nylon film, 
which effectively formed part of the cathode of a cylindrical proportional 
counter. That the film caused negligible backscattermg was shown 
by reversing it, this giving no detectable decrease in count rate. A 
saturation thickness of steel or aluminium was placed behind and in 
contact with the film. The backscattering factor was then given by 
the ratio of the increase in count rate with the backing to the rate without 
backing. 

A second method was evolved to meet the criticism of Jaffe (1950) 
that a source having a film of nylon between it and the backing is not 
necessarily equivalent to a source deposited directly on the backing. 
Sources emitting both gamma and beta radiation were employed and 
the count rates of sources deposited on nylon and directly on the backing 
material were compared, the relative source strengths being obtained 
from the gamma emission rate. Essentially the same results were obtained 
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by both methods and there was no suggestion of a systematic variation 
of backscattering factor with the method employed. 

For the first method, ®Co, 22Na, 24Na, ®°Sr, and **P were employed, 
but 8°Sr and 8P were unsuitable for the second method. The results 
are shown in fig. 9 where the points indicate, in most cases, the average 


Fig. 9 
0:60 
24 89 32 
Coos Nine Na Sr Pp 

o 4 L 
Oo ae STEEL oe t 
g at F t 
x ! 

O:40F E 
0 
72 ALUMINIUM ra = 
a a 
w HE T 
- 4 
0:30 
oO 
uw) 
x 
ee 
% 0-20 [ 

ORO 

4 — 
xe) 2:0 


MAXIMUM BETA ENERGY IN MEV 


Backscattering of electrons by steel and aluminium source supports. Circles 
and squares indicate results for steel and aluminium respectively. 


of several determinations by both methods. The errors take into account 
the scatter of the individual results and are larger than those due to the 
statistics of counting alone. In the case of steel, there appears to be no 
significant energy dependence, in the region investigated, though there 
may be a slow increase with energy for aluminium. There appears to 
be some slight justification for believing that the backscattering factors 
for positrons (?2Na) are somewhat lower than for negatrons. Assuming 
no energy dependence for steel and aluminium the average values obtained 
are 

Steel ..  .. 0-449-+0-009 

Aluminium ..  0-348-+ 0-075. 


In calculating the absolute decay rates, a value for steel of 0-450--0-015 
was taken. 
§ 7. CONCLUSION 
The experiments described here yielded values for the branching ratio, 
absolute decay rates and half-life of #K. The gamma and beta ray 
energies appear to be well established and the form of the beta spectrum 
has been carefully studied by several investigators. At present, however, 
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the evidence for the mode of the decay is somewhat indirect or lacking 
in precision and it is felt that it would be worthwhile carrying out further 
experiments to verify the presently accepted decay scheme. 
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§ 1. INTRODUCTION 


THE work described below was performed in the Physics Laboratories 
at University College, in order to add to the available experimental 
information on the multiple scattering of electrons and positrons in nuclear 
emulsions in the range of energies between 1 and 3 Mev. Where it was 
found possible, plates were exposed under nearly identical conditions to 
electrons and positrons of the same energy, in order to enable close 
comparison between the behaviour of the two types of particle to be made. 


§ 2. EXPOSURE OF PLATES 


The observations were made on six plates chosen from a series which 
were exposed to nearly mono-energetic electrons and positrons obtained 
from radio active sources. 

Plates No. 1 and No. 2, (2 in. x 2in., Kodak NT4), were irradiated in 
a flat magnetic spectrometer. The magnet was one which had been 
designed for a microtron (Henderson et al. 1953) and was utilized for the 
present purpose prior to completion of this machine. The magnetic field 
distribution of this magnet is uniform to better than +0-05° over the 
region used for the present experiments, and its absolute calibration at 
the time of use is known to better than +0-5°%%. 

The vacuum chamber of the spectrometer was designed such that the 
radioactive source and the centre of the photographic plate respectively 
were at the 180° points of particle trajectories with a mean radius of 
10-5 cm. 

Baffles were inserted between the source and the plate-holder such as 
to limit the direction of entry of electrons to --6° from the mean direction. 

The plate holder was designed to allow the electrons to enter the emulsion 
at a glancing angle of incidence of 8°. 

Because of the finite size of the plates, the incident energy of the 
electrons varied by -+-6% over the area scanned. The tracks were 
found to be distributed uniformly enough over the scanned area to enable 
a simple correction for this energy spread to be made. 

Plates No.3, No.4, No.5 and No.6, (3 in. x 1 in. Ilford G5), were exposed 
to radioactive sources in a lens spectrometer which has been fully described 
elsewhere (Rundle et al. 1954). This spectrometer had been equipped with 
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a helical baffle system such that either electrons or positrons could be 
focused on to the photographic plate by a reversal of the magnetic field 
direction. 

The absolute calibration of this spectrometer is known to better than 
+0-5%, and the scanned area of the plates was chosen such as to keep the 
total spread in energy to less than +3%. 

The plates were exposed as shown in table 1. 


Table 1 
Plate No. 1 » 3 4 5 6 
Plate type NT4 NT4 G5 Gd G5 G5 
Particles Electron | Electron | Positron | Positron | Electron | Positron 
Approx. energy 
in Mev 1-0 1-5 1-0 1-5 2°5 2°5 


§ 3. OBSERVATIONS 

The plates were scanned on Cooke M 4005 microscopes, using 45 x oil- 
imniersion ‘objectives and 15 x eyepieces. 

In all cases an area less than 1 in. <}in. near the centre of the plate 
‘was scanned, in order to exclude tracks in the possibly distorted emulsion 
at the edges of the plates as well as to minimize the energy spread of the 
particles. - 

Only tracks with an initial direction within --6° of the expected mean 


- direction were taken into account. A subsidiary experiment indicated 


that this angular range would exclude the bulk of the background tracks. 

At the relatively low energies involved in these experiments, the 
tracks deviate so much from linearity that the coordinate method of 
making multiple scattering measurements due to Fowler. is not feasible 
unless inordinately small cell lengths (< 5 microns) are used. Accordingly, 
the observations were recorded in terms of the projected angles between 
successive chords measured by means of an ocular protractor. Again, 
because of the large scattering angles observed, no correction for micro- 
scope noise was necessary.. A chord length of 50 microns was used in 
all six cases. 

The number of observations made on each plate varied from 672 to 
1102. In the plates exposed to 1-5 Mev and 2-5 Mev particles, measure- 
‘ments were made on no more than four 50-micron chords on each track. 
No correction for energy loss was made. On the two plates exposed to 
1:0 Mev particles, because of the considerable difficulty in finding tracks 
not scattered out of the emulsion, readings were made (where it was 
possible) on up to ten 50-micron chords ; in consequence a correction for 


energy loss was essential. This was done by calculating the theoretical 
tate of loss of energy for electrons at 1-0 Mev (Ross and Zajac 1949) 


and assuming the rate to vary little over the range, calculating the mean 


energy over the total length of track under observation. 
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$4. RESULTS 


The angular distribution on each plate was examined and found to be 
approximately Gaussian. The number of readings on each plate (approx. 
800) was insufficient to enable more exact comparison with the theoretical 
shape of the distribution curve to be made. However a composite distribu- 
tion of scattering on all plates was obtained by summing the distribution 
in intervals of one-quarter of their separate r.m.s. scattering angles. This 
is shown in fig. 1 compared with a Gaussian of the same r.m.s. angle and. 


Fig. 1 
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Composite histogram of all electrons and positrons. 
(a) Moliére distribution. (b) Gaussian distribution. 


with a Moliére distribution with the same mean scattering angle. Similar 
composite histograms are shown for electrons and positrons separately in 
figs. 2 and 3 respectively. 
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The values of the scattering constants with and without cut-off have 
been calculated for the readings on each plate and are shown compared 
with the values predicted by Williams and by Moliére in table 2. Ky, 
was calculated directly from information supplied by the makers OE the 
respective emulsions, by the method described by Moliére (1948), and the 
method of smoothing in converting from tangent to chord angles Aeceribe 
by Menon et al. (1951). | 


Fig. 2 
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Composite histogram of all electrons. 
(a) Moliére distribution. (b) Gaussian distribution. 


All the values shown, both theoretical and experimental are for 50-micron- 
chord lengths. At both 1 and 1-5 Mev, the experimental positron- 
electron differences are largely caused by the different emulsions used ; 


216 F. F. Heymann and W. F. Williams on the M ultiple 
Fig. 3 


\ 
100 q 


Number of scatters 
oOo 
pe 


2 S 4 5 
Angle in units of R.M.S. scattering angle. 


Composite histogram of all positrons. 
(a) Moliére distribution. (b) Gaussian distribution. 


Table 2. Experimental and Theoretical Values of the Scattering Constant 


Energy 
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NT4 and G5 respectively. No calculation of the effect of different charge 
has been included in the theoretical calculations of K ypojaro. 

These results show two features of interest. The experimental 
scattering constants appear to decrease with decreasing energy to values 
appreciably lower than those predicted theoretically. This has been 
earlier remarked by Hisdal (1952) and by Overds (1954) and the graph of 
the results of many workers (fig. 4) shows a possibility that this may be 


Fig. 4 
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Results of many workers. 
Ringed points, electrons. Crossed points, positrons. 


1. Hisdal (1952) 5. Bosley and Muirhead (1953) 
2. Present work 6. Gottstein et al. (1951) 

3. Husain (1955) 7. Menon et al. (1951) 

4. Bosley and Hughes (1955) 8. Corson (1951) 


a real effect. However this impression given by fig. 4 rests largely on 
the results of the present workers and of Hisdal and Overos in all of which 
the low values may be caused by a combination of the * thin pharm bers 
effect discussed by Moliére and the difficulty of following the more boavaly 
scattered tracks biassing the results towards low values of the scattering 
constant. These considerations should not however affect the readings 
at 2-5 mev where the absolute mean scattering angle was only of the order 
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of 4°. The second feature of interest arises from the comparison of the 
scattering of positrons and electrons shown in table 3. 


Table 3. Electron: Positron Ratio of Scattering Constants 


Energy in Mev 1-0 1-5 2-5 


ONE Ee 1-08+0-04 | 1-:05+0-04 | 1-17+0-04 
Pd GS ard Ge 1-08 0-04 | 1-:04+40-04 | 1-14+0-04 


| 


The apparent difference at 1 and 1-5 Mev is largely attributable to the 
different emulsions used. This and other differences in conditions at 
these energies caused by availability of equipment made accurate com- 
parisons invalid at these energies. At 2-5 Mev however considerable 
care was taken to obtain similar conditions and G5 emulsions were used 
for both particles. The difference found here is considerable and of 
similar order to those reported by other workers (Bosley and Hughes 
1955, Corson 1951, Groetzinger, Humphrey and Ribe 1952). Furthermore 
as is shown by the ratio of cut-off scattering constants and by the separate 
histograms of the positron and electron distributions—fig. 5—this 
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Histograms for 2-5 Mev particles. 
Full line, electrons. Dotted line, positrons. 
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difference is not caused primarily by the increased frequency of large angle 
scatters in the case of electrons. At these energies single scatters of the 
order of 20° or more are easily visually observed and if these are ignored 
in calculating the scattering constants the same difference remains. 

In addition to the readings recorded and discussed above a reading was 
taken of the angle between the tangent to each track and a 50-micron 
chord. The absolute mean value of this angle was found to be one-half 
of the absolute mean value of the angle between successive 50-micron 
chords within the statistical error. The ratio between the absolute mean 
values of this angle for electrons and positrons at 2-5 Mev was found to 
be 1-15-L0-08. 


§ 5. CONCLUSIONS 


The present work shows a larger difference between the multiple 
scattering of electrons and that of positrons at 2-5 Mev than is predicted 
by theory (Mohr 1954). Both the electron—positron difference at lower 
energies and the falling off of the scattering constant at low energies 
(for which 6? <1) require further investigation and it is hoped to examine 
both these points and to bridge the gap between the present work and 
that of Hisdal (1952) and of Overés (1954) by making measurements on 
both positron and electron plates exposed and developed under similar 
conditions to 0-75 Mev and 1-5 Mev particles. 
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XX. REVIEWS OF BOOKS 


Proceedings of the International Conference of Theoretical Physics, Kyoto and 


She 1953. (Tokyo : Science Council of Japan, 1954.) [Pp. xxviti+942.] 
$10. 


THis volume records the first purely scientific international conference to be 
held in Japan. Over one hundred speakers contributed with original papers 
or informal reviews of current fields of interest. The following topics were 
discussed :— 

Field theory; cosmic rays; pions; intermediate and strong coupling 
theory ; nuclear forces; nuclear reactions; shell structure and beta-decay ; 
polymers ; liquids; transport phenomena ; irreversible processes : general 
methods in statistical mechanics ; crystal dislocation; molecules; metals : 
electron theory of intrinsic magnetization ; antiferro- and ferromagnetism : 
magnetic resonance; dielectrics; colour centres; liquid helium; super- 
conductivity. J.M.Z. 


Determination of Organic Structures by Physical Methods. Edited by E. A. 
BraupeE and F. C. Nacuop. 22 Contributors. (New York: Academic 
Press Inc.; London: Academic Books Ltd.) [Pp. xiii+810.] Price 
£5 7s. Od. : 


Tuts book should be warmly welcomed by organic chemists and by those 
physicists and physical chemists who specialize in the application of physical 
methods to the determination of the structures of molecules. The majority of 
such methods were originally developed by physicists for their own interest, 
and were only subsequently applied to organic compounds, and although often 
described individually in text-books, there was, until the appearance of this 
volume, no collective review of the principles involved in their application. 
Part I contains two chapters dealing with methods based on the phase properties. 
of molecules, and which are of importance in the determination of molecular size. 
Part II consists of five chapters describing the application of optical and magnetic 
measurements to the determination of general structural features. The third 
and longest part contains nine chapters dealing with methods which have 
contributed information about the fine structure of molecules, especially details 
of stereochemistry and electron distribution. The sixteen topics treated thus 
range from ‘ Phase Properties of Small Molecules ’ through ‘ Optical Rotation ’ 
and ‘ X-ray Diffraction ’ to (somewhat surprisingly) “ Wave-Mechanical Theory ’. 
A glance at the list of contributors shows the authors’ remark that “ each is 
an active and acknowledged leader of research in the field he surveys ”* to be 
justified. The book contains nearly 2000 references to publications, and is well 
produced but expensive. WC. 


Quantum Mechanics. Second edition. By L. I. Scuirr. (McGraw-Hill Book 
Company, Inc.) [Pp. xii+417.] Price 47s. 


Tu second edition of Professor Schiff’s well-known book on quantum mechanics 
is welcome as this is one of the most useful books for introducing a serious student 
to the subject. The book is substantially the same as the first edition though 
several minor changes have been made improving the presentation and in places’ 
bringing it up to date. 
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The Laws of Nature. By R. E. PEIERLS. (George Allen and Unwin.) 
[Pp. 284.] Price 21$. 


PROFESSOR PEIERLS was led to write this book from recent experience in 
lecturing to various audiences on atomic energy. Among his listeners he 
found some who were infected with the same scientific curiosity that scientists 
possess, but who were without the background of knowledge that a quite 
general scientific training would have given them. 

In this book, under its self explanatory title the author sets out to provide 
information for intelligent non specialists. He is right in thinking that the 
reading of it will come more easily to those who possess some elementary 
knowledge of Physics. But he is to be congratulated on the success of his 
labours. Indeed professional physicists will find in the book modes of presenta- 
tion of well-known topics which may be useful also to them if they instruct 
the young in the class-room or find themselves like Professor Peierls facing 
a wider public. A.M. T. 


BOOK NOTICES 
Principles and Applications of Physics. By O. Burin and J. D. Exper. 
(Harvard: University Press.) [Pp. 866.] Price 45s. 
Progress in Biophysics and Biophysical Chemistry. Vol.5. By J. A. V. BuTLER 
and J. T. Ranpaty. (Pergamon Press.) [Pp. 231.] Price 55s. 
Classical Electricity and Magnetism. By W. K. H. Panorsky and MELBA 


Puintures. (Cambridge, Mass.: Addison-Wesley Publishing Co.) [Pp. 400.] 
Price $8.50. 


Linear Feedback Analysis. By J.G.THomMason. (Pergamon Press.) [Pp. 355.] 
Price 55s. 


Introduction to Electronic Analogue Computers. By C. A. A. Wass. (Pergamon 
Press.) [Pp. 231.] Price 22s. 

The Unified System Concept of Nature. By 8. T. Bornemisza. (New York: 
Vantage Press Inc.) [Pp. viii+137.] Price $3. 

The Elements of Probability Theory and some of its Applications. By Haraup 
CRAMER. (London: Chapman and Hall.) [Pp.281.] Price 56s. 


Theories of Perception and the Concept of Structure. By F. H. Aiport. 
(London ; Chapman and Hall.) [Pp. xxii+709.] Price 64s. 


University Physics. By F. W. Sears and M. W. Zemansky. (Cambridge, 
Mass. : Addison-Wesley Publishing Co., Inc.) [Pp. viii+1031.] Price $10. 
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N. THOMPSON et al. Phil. Mag. Ser. 8, Vol. I, Pl. 3. 


Figs. 4-7: Development of fatigue crack. 


(4), after 74% of life; (5), as (4), electropolished; (6), after 42% of life, 
electropolished ; (7), after 77% of life, electropolished. 


N. THOMPSON et al. Phil. Mag. Ser. 8, Vol. I, Pl. 4. 


Figs. 8-10: Another part of the same specimen. 
(8), after 42°, electropolished ; (9), after 77°%, electropolished ; 
(10), at conclusion of test (3-6 x 10® evcles), not polished. 


Figs. 11 and 12: Partially fatigued specimen, electropolished, before (11) and 
after (12) 5° tensile strain. 


_ N. THOMPSON et al. Phil. Mag. Ser. 8, Vol. 1, Pl. 5. 
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Figs. 13-15: Development of fatigue crack in single crystal. 
(13), after 18%, electropolished ; (14), after 86% ; (15) as (14), but electropolished. 


N. THOMPSON et al. Phil. Mag. Ser. 8, Vol. |, Pl. 6. 


Figs. 16-18: Development of fatigue crack in single erystal. 


(16), after 54%, electropolished ; (17), after 90°%, electropolished ; 
(18), at conclusion of test (108 cycles). 


N. THOMPSON et al. Phils Magaserac, Vol. |; Pl.7. 


Fig. 19 


Fig. 23 


(a) (0) 
Fig. 19: Slip on three planes : single crystal. 


Figs. 20-22: Effect of deep electropolish on persistent slip bands. 
(20), after removing 2 p. from surface; (21), after removing 10; (22), after removing 18 p. 


Figs. 23 (a) and (b): Material extruded from slip bands of fatigued single crystal. 


N. THOMPSON et al. Phil. Mag. Ser. 8, Vol. I, Pl. 8. 


Fig. 24 Fig. 25 


Figs. 24-27: Development of fatigue crack in specimen tested in nitrogen. 


(24), after about 1°, of life in nitrogen : electropolished ; (25), after 74°, electro- 
polished; (26), after 30%, electropolished; (27), after 73%, electropolished. 


lice BULLI Phil. Mag. Ser. 8, Vol. 1, PI. 9 


Fig. 3 


(a) 
Oscillogram of a compression pulse in water. Upper trace: 100 kce/s timing 
wave. Lower trace: amplified voltage output from condenser micro- 
phone. 


(0) 
Oscillogram of pressure changes 1 cm below a free water surface during re- 
flection of a compression pulse. Upper trace: 100 ke/s timing wave. 
Lower trace: amplified voltage output from condenser microphone. 


Je Fe BULE Phil. Mag. Ser. 8, Vol. I, Pl. 10. 
Fig. 5 
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Spark shadowgraph of cavitation in water. The times shown below each 
photograph indicate the interval between the firing of the detonator 
and the emission of light from the spark gap. 
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Fig. 6 


Cavitated zone produced in water by reflection of a compression wave at the free 
surface. 
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